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ABSTRACT 


Dynamic positioning of a moving platform is analyzed using data from the 
Global Positioning System (GPS) acquired in Phase II of the Seafloor Benchmark 
Experiment on R/V Point Sur in August 1986. GPS position determinations are 
compared to simultaneous Mini-Ranger fixes. 

Ihe GPS positions computed using only broadcast ephemeris data were within 
20 m from the Mini- Ranger fixes when data from four satellites were used and within 
30 m when data from three satellites and a geoidal height constraint were used. It was 
found that the position accuracy is degraded when data from a satellite reaching 


culmination is used. 
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THESIS DISCLAIMER 


The reader is cautioned that computer programs developed in this research may 
not have been tested for all cases of interest. While every effort has been made, within 
the time available, to ensure that the programs are free of computational and logical 
errors, they cannot be considered validated. Any application of these programs 


without additional verification is at the risk of the user. 
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I. INTRODUCTION 


The Seafloor Benchmark Experiment, a project of the Hydrographic Sciences 
Group of the Oceanography Department at the Naval Postgraduate School (NPS), was 
initiated in 1984 with a goal to establish stations on the seafloor in real-time with 
geodetic accuracy (Saxena, 1982). The first experiment (Phase I), conducted in May 
1985, included a configuration of four benchmarks on the seafloor (Spielvogel et al., 
1987). Further study indicated a sizeable improvement in the positional accuracy of the 
inner stations as compared to the accuracy of the outer stations of a larger net 
configuration (Kumar and Saxena, 1985). For this reason, during the Phase II 
experiment in August 1986, a nine-station net configuration was used (Kumar, 1986). 
Three types of GPS receivers and a Mini-Ranger (MR) Falcon system were used to 
determine the position of the R/V Point Sur. 

Keeping the goal of establishing and continuing GPS research at NPS, Brown 
(1986) established the Texas Instruments 4100 GEOSTAR GPS software, written by 
the Naval Surface Weapons Center for their CDC Cyber 865 computer, on the NPS 
IBM-3033 computer. Brown was able to validate the modified software using static 
position data, but not the KALMN2 program which uses dynamic position data. 

The main objective of this thesis 1s to determine the positions of a dynamic 
platform using GPS and to compare them with MR Falcon-determined positions. 
Other objectives were to install software for processing the raw data from cassettes to 
FIC (Floating-Integer-Character) format, modify the KALMN2 program and validate 
it using Phase II data. Programs to plot and analyze the continuous ship positions bv 
GPS and by MR Falcon were also written. 

Ihe data used in this thesis were acquired during Phase II of the Experiment. 
Ihe GPS data are limited to those collected by the TI-4100 receiver from NSWC 
installed on board of the R/V Point Sur. The day chosen for analysis was 16 August 


1986, since adequate pitch and roll data are available for that day. 
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IL. INSTRUMENTATION 


A. POSITIONING SYSTEMS 
l. GPS receivers 
During the experiment, the following GPS receivers belonging to several 
agencies and companies were installed on the R/V Point Sur: 
- TI-4100 receiver from Naval Surface Weapons Centers (NSWC) 
- T1-4100 receiver from National Ocean Survey (NOS) 
- Eagle Mini-Ranger receiver from Motorola, Inc. 
- Trimble 4000A from Trimble Navigation 
There were shore based GPS receivers at the following stations: 
Beach Lab 
- TI-4100 from Pacific Missile Test Center (PMTC) 
- Trimble 4000A from Trimble Navigation 
Ferrier 
- Eagle Mini-Ranger from Motorola Inc. 
Dome 


- Magnavox Manpack from GPS Joint Program Office. 


2. Other positioning System 
For navigation and positioning on board were a MR Falcon svstem from 


Motorola, [nc., and a Loran-C receiver. 


B. OTHER SYSTEMS 

An acoustic positioning system from Oceano Instruments to acquire acoustic 
ranges (Kuo,1985) was installed on the R/V Point Sur; it also gives the position of the 
ship in relation to the bottom benchmarks (acoustic transponders) and can be 
integrated with a satellite receiver for computation of their positions using the GPS 
coordinates of the ship's acoustic transducer. 

In order to allow for corrections due to the pitch and roll of the ship, there was a 
data acquisition system, based on an HP9826 computer, to collect pitch, roll and 
heading data every second. Due to interfacing problems there was no possibilitv of 
logging heading data. The ship's Data Acquisition System (SDAS) collected, among 


other data, the ship's heading with a sampling period of about 19 s. 
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С. ANTENNAS AND TRANSDUCER POSITIONS 

The antenna $f the MR Falcon was mounted on the mast of the ship in order to 
increase its range. The antenna of the TI-4100 from `NOS was installed close to the 
MR Falcon antenna. All the other GPS antennas were mounted on an elevated 
wooden table above the ship's laboratories and close to the transducer position (Figure 
24b: 


MR Falcon antenna 


Platform with GPS antennas 


Acoustic transducer 


Not. to scale 





Figure 2.1 Relative positions of antennas and transducer. 
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ІП. DATA PROCESSING 


A. COORDINATE SYSTEM 
|. Parameters 

In order to compare the positional accuracies of the TI-4100 GPS receiver 
with the Mini-Ranger (MR) Falcon, computations were made in a single coordinate 
system. 

Although there was a Doppler survey to find the WGS 72 coordinates of all 
control stations used, due to the impossibility of setting a station over Dome Ecc, the 
site for one of the MR Falcon shore stations, the WGS 72 coordinates of nearby 
station Dome, 68 m away, were obtained. During the Phase II experiment, GPS 
provided coordinates in the WGS 72 system (Bomford, 1980, Sepplin, 1974), although 
it converted to WGS 84 on | January 1987 (Decker, 1986). 

The coordinates of Dome were computed in NAD 27 and then transferred to 
Dome Ecc. The coordinates of Dome Ecc were converted from NAD 27 to WGS 72 in 
order to compute directly the ship positions in this system, avoiding the problem of 
converting each position individually. 

The subroutines used for converting geodetic coordinates from WGS 72 to 
Universal Transverse Mercator (UTM) and vice-versa wer written in Fortran and are 
based on Basic subroutines used by DMA for an HP9826 desktop computer. A cross 
check was made which revealed a maximum conversion error of 0.01 m, which is well 
within the accuracy limits of the system and, hence, does not affect the analysis of the 
data processed. 


The following parameters were used: 


WGS72 
Semi major axis : 6 575155 m 
Flattening : 1/298.26 


UTM projection 
Central meridian: 123° 00’ 00” W 
If this package of subroutines is used for any location other than the area of the 
Seafloor Benchmark Phase II, the UTM parameters must be changed in order to get 


the correct conversions. 
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Figure 3.1 Area of the experiment. 


2. Coordinates of stations 


The two shore stations were located on high points along the coast (Figure 
3.1). Their locations and coordinates in WGS 72 are as follows: 
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Ferrier 
Γαι; ΘΠ 
Longitude 121° 33' 527.939 W 
Dome Ecc 
Latitude 36? 18' 207.897 N 
Longitude 121? 54' 007.594 W 
The elevations of the different antennas were set in the Mini- Ranger receiver and the 
measured slant distances were converted 


to appropriate horizontal distances by the receiver system during data acquisition. 


B. REDUCTIONS OF POSITIONS TO THE TRANSDUCER 
l. Offsets of transducer relative to antennas 


Since the two antennas, namely the MR Falcon antenna near the top of the 
mast and the TI-4100 antenna on the elevated table above the lab on the R,V Point 
Sur, were at different locations, a common point had to be chosen to compare both 
systems. The acoustic transducer mounted beneath the hull was chosen as a common 
point, which was used for further integration with the acoustic data, to compute the 
position of the Seafloor Benchmarks. 

For data reduction there was a need for pitch, roll and ship's heading data. 
While the existing pitch and roll data were sampled with a period of 1 5, the ship’s 
heading was sampled with a period of about 19 s. Computed headings were tried to 
improve the rate but tests indicated that the real data, even with a low sampling rate 
was better than the computed. The method used to compute the heading 1s described in 
Section Сэ. 

The geometric offsets between the MR Falcon, TI-4100 antenna and the 
acoustic underwater transducer were computed using a right-handed coordinate system 
centered in each antenna, and having the Y axis oriented towards the bow, the X axis 
to starboard and the Z axis towards the zenith. Using this convention, all offsets of 
the transducer are negative. Using the formulas developed in the next section, the 
corrections for pitch, roll and heading data were computed and applied. 

For the MR Falcon positions only horizontal corrections were applied, since 
the positions are computed in a two-dimensional system. For the GPS positions the 


corrections were applied to all three coordinates. 


For both types of positions the coordinates are converted to UTM, corrected 
and converted back to WGS 72. For the GPS positions the Z correction is applied to 
the geoidal height. 

The following offsets were measured and computed by Prof. Tucker and Mr 
James Cherry, and changed to the above format. 

Transducer offsets relative to the GPS antenna: 

XOFF + -4.268 га 
ΥΟΓΕ δα πι 
ZOFF = -9.505 m 
Transducer offsets relative to the MR FALCON antenna : 
XOFF = -5.345 m 
YOFF = -9.453 m 
ZOFF = -16.152 m 


И 


The course was computed in a range 0° to 360°, and the pitch and roll data 
were stored in a range -90? to +90°; following this convention the roll is positive when 
the port is up and the pitch positive when the bow is down. The sign of roll was 


changed to conform with the mathematical model described next. 


2. Math model 
The method to compute the correction to the coordinates 1s based on a seven- 
parameter transformation described in Moffitt and Mikhail (1980), so successive 
rotations of the vector defined by XOFF, YOFF and ZOFF, give the corrections to the 
coordinates illustrated in Figure 3.2. 
The rotation matrices are defined as follows: 


Mp - 1 0 0 
0 cos(PITCH) sin(PITCH) 
0 -sin( PITCH) cos(PITCH) 
M cos( ROLL) 0 -Sin( ROLL) 
0 1 0 


sin( ROLL) 0 cos(ROLL) 
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p - pitch r -roll h - heading 


A = Antenna T=Transponder 


Figure 3.2 Rotation of axis. 


Np cos( HEADING) sin( IEADING) 0 
-sin( HEADING) cos(« HEADING) 0 
0 0 l 


The orientation matrix is defined by the following equation: 


ше = M, M, Mj (eqn 31) 


and the corrections are computed by multiplying the orientation matrix by the offset 


vector, to get the following equations: 


DX = XOFF cos(ROLL) cos(HHEADING) + (ео 2) 
YOFF (cos(PITCH) sin( HEADING) + 
sin(PITCH) sin(ROLL) cos(HEADING)) + 
ZOFF (sin(PITCH) sin( HEADINGO]) - 
cos(PITCH) sin(ROLL) cos(IIEADING)} 


DY = XOFF cos(ROLE) sin( HEADING (i=. (3422) 
YOFF {tos(PITCH) cos( HEADING) - 
sin(PITCH) sin ROLL) sin( HEADING)) + 
ZOFF (sin ΡΙ ΟΠ ο ΓΕ ο ο) ` 
cos(PITCH) sin ROLL) sin( HEADING)) 


DZ = XOFF sin ROLL) - (едп 3.4) 
YOFF sin(PITCH) cos(ROLL) + 
ZOFF cos(PITCH) cos( ROLL) 


C. FALCON MINI-RANGER DATA PROCESSING 
1. Error of system 
The MR Falcon was used during the cruise for navigation and positioning. In 
post-processing it is used to evaluate GPS data. Since the Mini-Ranger is a system 
frequently tested and used, its error is known to be = 3 m (Laurila, 1976: Munson, 
[ОЛ 
It would be betier to have had more than two lines of position (LOP) to allow 
a minimization of the error in the computation of the Mini-Ranger position and also 
to allow a smoothing of the data in places where the data gets noisy as was seen in 
post processing. In Figure 3.3 it is possible to see the areas where the data is noisy and 
to see the small jumps characteristic of the svstem. 
2. Antenna positions 
For computing the antenna positions the coordinates of the shore stations 
were converted to WGS 72 UTM, and then each position was computed in the WGS 
72 UTM projection. The subroutine used to compute the positions was picked from 
program UCOMPS, a utility package used in the Hvdrographic Sciences Group of the 


Oceanography Department. 


3. Heading Calculation 
As mentioned previously, the ship's heading was needed to compute the 
transducer's position, but there was no data available with the l-s periodicity needed, 
so headings were calculated from Mini-Ranger data. The calculation of the ship's 
headings was made during the computation of the two antenna positions and kept on a 
separate file. These computed headings were affected by the pitch and roll of the ship 


and also bv the noisy data of the Mini- Ranger. 


18 


e 
= 


1910 





Figure 3.3 Plot of ship's track using MR Falcon data. 
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Some trials were made using different running averages for smoothing the 
computed heading in order to get a course close to reality. It was found that the 
courses obtained, even if filtered by a running average, oscillated and hence were not 


suitable for this analysis. 


4. Transducer positions 
The transducer positions were computed using the model described in 
Section B. The pitch and roll data used were sampled with a periodicity of I s most of 
the time when there was positioning data. The ship’s heading was sampled with a 
periodicity of from 19 s to 1 minute. Although tests were made in order to compute the 
heading, as stated in the previous section it was decided to use recorded data, 
constraining the analysis to periods of constant course. 
5. Programs. Inputs and outputs 
During the development of the programs many more routines were written 
than mentioned here. [n order to make the use of the programs easy, the total number 
of programs to process the MR Falcon data was reduced to four. 
A unique format for the different position files was defined, allowing only one 
plotting program in which the user must change the title according to the type of 


positions plotted. 


Program: FALCON 
Input: Month, day, vear, hour, minutes, seconds (date time tag), 
codel, rangel, strengthl, codel, range2, strength2 
Output: Lisung with date time tag and geographic positions 
File with date time tag and geographic positions 
File with date time tag and rough course 


Source: See Appendix A 


Program: COURSE SMOOTH 

Input: File with date time tag and rough course 

Output: File with date time tag and smooth course 
Listing with date time tag and smooth course 


Source: See Appendix B 
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Program: TRANSDUC FALCON 
Input: -File date time tag and geographic positions of antenna 
Output: File with date time tag and geographic positions of transducer 
Optional listing with date time tag and geographic positions 
of transducer 


Source: See Appendix C 


Program: PLOT 

Input: File with date time tag and geographic positions 
Output: Plot with track of the ship in UTM projection 
Source: See Appendix D 


All files are in free format, and, if for by any reason a fixed format is needed, the user 
should change the programs to get the output desired. The MR Falcon data were 
originally on floppy disks in HP9816 format. A program was written in Basic in order 
to convert and write the data onto a 9-track magnetic tape. The transfer of data to 
mass storage on the NPS IBM mainframe computer and data manipulation were made 
using existing procedures for the mainframe (Favorite, 1986; Mar, 1984). 

The scale and the area of the plot are fixed for the particular area of the 
Benchmark Experiment but may be changed in the program; thus any proportional size 
can be obtained by changing the control card as is explained in the manual (W.R. 
hur h Computer Center, 19815. 

As the plot program can be used for plotting the outputs of different 
programs, the title should be changed accordinglv, leaving the appropriate lines 


uncommented. An hour-minute tag is plotted every 5 minutes. 


D. GPS DATA PROCESSING 
1. Reading the TI-4100 cassettes 
a. System configuration 

The data from the TI-4100 receiver is recorded on digital cassettes in a 
format not directly compatible with the existing programs at NPS. To read and decode 
the data, it was necessary to purchase a MEMTEC cassette terminal (Model 5450XL) 
and obtain the appropriate software from the Applied Research Laboratories (ARL) of 
the University of Texas at Austin. The MEMTEC cassette reader was connected to an 


IBM PC having a hard disk of 20 Mb (Megabytes), a modem and a math co-processor. 
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Preliminary tests made to verify the software and hardware connections of 
the tape deck wit the IBM PC indicated good results. During the tests it was found 
that the format of the final output of this package was not the same as defined 
previously by Scott and Peters (1983) but was in a new format, FIC, where each block 
is organized in a structured way with control records followed in order by all floating 
type data, integer type data and character type data. The first program, CON9TR, а 
part of the library installed by Brown (1986) in the SEF (Standard Exchange Format), 


was replaced by a new program using the FIC format. 


b. Reading cassettes to disk 
Each cassette is dumped to disk using the program MFERD. The operation 
takes about 20 minutes, and the file created fits on one floppy disk which can be kept 


for raw data backup. 


c. Converting the data 

The binary image of the cassette is then converted to a FIC format file 
using the program GS2FIC, which takes about 20 minutes. This format is suitable for 
processing in the IBM-PC with appropriate software to be developed. 

In order to transfer these files to the main frame or other computer, it is 
necessary to run a program FICFICA that converts the FIC files to an ASCII, format 
where the data 1s organized in 80-column records. This operation takes about 30 
minutes, and the file created by this program can be dumped to a magnetic tape or 
sent to a nost computer. 

The files created by these two programs cannot be stored on a normal 
floppy disk and were erased as soon as the data was on the mass storage of the 
mainframe and tested. 

d. Transfer to the mainframe 

The transfer of the files to the mainframe was done using a micro-computer 
connected via modem to the main frame. A terminal emulator distributed by the NPS 
Computer Center, SIMPC, simplifies the transfer of files (Simware Inc, 1984). 

Each file, corresponding to one cassette, was sent to a disk. This operation 
took about 3.5 hours, and sometimes the transfer was stopped in the middle. When this 
happened the transfer had to be started from the beginning. 

Each file has around 15000 lines and occupies about 65% of an A-disk. 


Using an extra disk every file was converted to a job format in order to transfer it to 


mass storage. The data were sent to a member of a partition data set where all 
cassettes were stored. 
2. Conversion of formats 

After a study of the FIC format, the input of KALMN2 program, and the 
program CON9TR, a program was written called CVFICA, which had, after some 
debuging and improvements, a final form that fulfilled the needs for processing the 
data for this thesis. 

CVFICA decodes each FICA block, producing two files and a listing of 
warnings and general information. One file contains the data to be processed by the 
KALMN2 program and the other has the positions computed by the GESAR software 
installed in the receiver. 

In the present version of the program the pseudo-ranges and other tracking 
information are stored only after the navigational data for the desired number of 
satellites 15 in the output file, since this data is needed to compute the positions in the 
Navigational Mode. [t is very easy to change this program to process the data using 
the precise ephemeris. The program flags as bad the data that belong to a Space 
Vehicle (SV) whose data does not correspond to the existing SV identification in that 
tracer. The data are also flagged as bad when the status vector has any value other 
than zero. 

The program is able to convert the following data block types: 

101 - GESAR Versions 1.0+ Input data, 
3 - GESAR solutions, 
6 - Tracking data, 
8 - Tracking configuration, 
109 - Navigation Message record as transmited, 
9 - De-blocked subframes 1-3 from block 109 
11 - Receiver error block, 


13 - Tape header/trailor. 


According to the CVFICA program any other block is read and dumped to 
the listing with a warning saying that the program was not able to handle it. No such 
warning was found after using the program with data from 10 cassettes. 

The listing also indicates when the navigational data for each satellite is 


received in order to give to the user the status of the constellation of satellites. 
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As this program outputs the GESAR position solutions and the data to the 
KALMN2 program into data set files, these must be created prior to running the 
program. 

A control file is used to pass to the program information such as two lines for 
the title of the listing and the starting date of the GPS week. The starting date allows 
the conversion of the time tag in seconds to a date time tag. The program does not 


work with data that crosses the one-week limit. 


3. Program information 

Program: CVFICA 

Input: FICA files without the first two comment lines 
Control file with the title and starting date of GPS week 

(Appendix J) 

Output: File with GESAR solutions 
File in NSWC format (input to KALMN program) 
Listing with warnings and general information 


Source: See Appendix E 


4. Computation of positions 
a. Problems found 

When starting to use the Kalman filter program (Brown, 1986) for 
processing the GPS data some difficulties were encountered, either due to the 
replacement of the CON9TR program by an equivalent one or by errors existing in the 
program KALMN2. The first problems were caused by the different units used in the 
ARL and NSWC software which were a consequence of the lack of documentation 
during the development of the program CVFICA; these problems were easily found 
and solved. The last problems took a long time to find, not only due to the complexity 
of the program, but also due to the nonexistence of processed data to check the results. 
A big effort was made in order to correct all errors, but only when two cassettes were 
sent to NSWC and processed there were all errors removed from the processing. It was 
found that the last two major errors were due to causes external to the KALMN2 
program, one an error in the ARL software that affected the ambient temperature by a 
factor of 10 and caused an excessive tropospheric correction and the other a bad 


constant in the control file of the KALMN2 program. 
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b. Changes made 

In addition to the corrections, some changes were made to improve the 
outputs. Now the program prints a date time tag just prior to all outputs labeled with 
time tag in seconds of GPS week. A utle is printed as a header of the output listing. 
Other small changes were made to improve the code. All changes were made carefully, 
and the programs were throughly tested. 

c. Data processing 

After inserung data from ten TI-4100 cassettes into the mass storage, a 
control file for the KALMN2 program was made (Brown, 1986). For an initial 
approximate position, the position of the receiver based on the GESAR solutions with 
the time tag closest to the starting time was used for the data processing. For the 
receiver time bias, an iterative procedure was used; thus, the data were processed using 
a value, starting with zero, that was replaced bv one listed in the outputs of the 
KALMN2 program when the computed positions were close to the predicted ones. 

Due to the good initial results of the test runs with certain satellite 
configurations, further tests were investigated for the different types of solutions. 

The value of -37 m was used as geoidal height in order to constrain the 
positions to the geoid when using data from less than four satellites or when forced by 
the control file. This value was obtained by interpolation on a chart (Blaha et al., 1986) 
and corrected for the antenna height above sea level, giving a corrected value of -51 m. 
The control file was set up in such a way as to save the positions on another file. The 
broadcast ephemeris was used, as the precise ephemeris was not available when the 
data were processed. 

d. Math model 
To compute position, the program uses the pseudo-ranges corrected for the 


various factors affecung them and an eight-state Kalman filter (Brown, 1986). 


e. Program information 
Program: “KALMN2 
Input: File in NSWC format produced by the program CVFICA 
Control file as defined by Brown (1986)(Appendix I) 
Control file with the title and starting date of GPS week 
Output: Listing with positions, errors and other information 
as selected by the control file 
File with time tag in seconds of GPS week, X, Y, Z 
and Lat, Lon, HT, both in WGS 72 
and the number of satellites used in the solution 
Optional file with time tag 


and Geometric Dilution of Precision (GDOP) every minute 


5. Transducer positions. 


a. Procedure 


The positions of the transducer were computed using an algorithm similar 
to the one used with the MR Falcon data. The only difference 15 that the GPS 
positions are three-dimensional, so the corrections for all dimensions were computed 
and applied. 

The offsets between the transducer and the antenna and also the format of 
the input files were different. The program uses the same control file as the program 
CVFICA in order to compute the date time tag of each position. 


b. Program information 


Program: | TRANSDUC GPS 


Input: File with time tag in seconds and WGS 72 coordinates 
Control file with the title and starting date of GPS week 

Output: File with date, time tag, latitude, longitude and ellipsoid height 
Optional listing with above information 

Source: See Appendix F 
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6. Plots 
I[t is possible to plot the track of the ship using either the GESAR solutions or 
the transducer positions as described before in Section C.5. Track plotting was done in 
order to give the GPS data coverage and data status. The program to plot is the same 


one used to plot the positions computed from the MR Falcon data but with the 
appropriate title (Figure 3.4). 
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Figure 3.4. Plot of ship's track using TI-4100 data. 
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IV. DATA ANALYSIS 


A. COMPUTATION OF DIFFERENCES 
1. Approach 

To find the accuracy of the GPS positions in a two-dimensional coordinate 
system, the observed positions in WGS 72 coordinates are transformed to UTM 
coordinates and compared with the UTM coordinates of the corresponding MR Falcon 
positions interpolated for the same time tag. The differences in the coordinates of 
Mini- Ranger and GPS positions are expressed in meters. 

Due to the low sampling rate of the existing course data, only data from 
steady, straight courses were used for analvsis. The data were divided into eleven 
periods of constant headings, while the data collected during the turns of the ship were 
deleted. This procedure made it possible to analyze the differences between GPS and 
Mini- Ranger data during the different courses. 

Due to unexpected small differences found for positions computed using data 
from three satellites, several tests were made using various combinations of solutions. 
Although for data analvsis many combinations were studied, three representative cases 
are discussed. 

Different starting times during processing of the GPS data using the 
KALMN2 program were used in order to see the effects of the propagation of errors 
due to the initial noisy data, but the discrepancies were inconsistent. The term “noisy 
data” is used to refer to the oscillatory behavior resulting from the Kalman filter when 
the number of satellites changes from three to four or vice versa. The effect of ship’s 
motion on the transfer of antenna positions to two chosen common points, 1.е. the 
acoustic transducer and the GPS antenna, was found to be 0.2 m, which falls within 
the accuracy limits of GPS and hence is considered negligible for this study. 

2. Case studies 

The following cases were analyzed: 

Case A : Data from four satellites were used or data from only three available 
satellites and geoidal height constraint for position computation were 


used. 
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Case Be Same as case A but using the geoidal height constraint when both 
three or four satellites were used. 


Case C Same as B but not using data from SVII. 


3. Differences obtained 
a. Period 1. 19:14 to 19:20 (Cases À and B) 


TABLE 1 


PERIOD 1, CASES A AND B. 
DIFFERENCES BETWEEN GPS AND MINI-RANGER POSITIONS 


Ranges(m 
AVg i 9) Мах 

ров 5У5 А В А В 
179 6,8,11 20 20 24-2) 24-32 
80 6,8,9,11 24 30 του 29-31 

55 277 30 Ме? 30-30 * 

| 6,8.9.11 34 40 9 124 ДУО 

1 6.8.11 34 42 34-34 42-42 
23 6,8,9.11 27 31 25-29 30-33 | 

| ος 1 28 31 28-28 SUCEDE: 
74 6,8,9,11 29 30 20-30 28-32 


* noisy data 

Case A - using four satellites, or. | | 
four satellites with, geoidal height constraint 

Case B - geoidal height constraint in case 


pos - number of consecutive positions computed every second 


During the first part of this period, all positions were computed using data 
from three satellites. In both cases (Table 1) the differences found are in the same 
range, 24 to 32 m. For the next 80 positions, computed using data from four satellites, 
Case A lead to values lower than in Case B. The upper limit of the ranges in cases, 31 
m, 15 due to the initial oscillation from the Kalman filter. 

The next part of the period is one of noisy data, i.e. changes from four 
satellites to three and vice-versa. In both cases the positional accuracy is degraded, 
giving larger differences between GPS and Mini-Ranger positions. 

At the end of the period, with continuous data from four satellites and with 
the exception of one position, the differences are low, reaching 20 m in Case A and 28 


m in Case B. 


It was found that the use of geoidal height as a constraint for position 
computations with data from four satellites degrades the computed position. The 
reason for this is explained in Section 4.B. Noisy data causes oscillations in the 


observed differences which are due to the Kalman filter. 





TABLE 2 
SATELLITE VISIBILITY 


| 16 AUGUST 1986 
Time ως, ЕП Az В” AZ Ви A B iv 
19:00 39 ντ... τ... 57 156 -- -- 
19:20 63 21 44 44 32 321 £67 147 6 198 
19:40 66 45 36 38 39 328 75 124 15 200 
20:00 66 72 28 34 48 335 77 78 24 203 


19 57 342 711 74 
11 67 350 62 36 





The offset of the GPS positions in relation to the MR Falcon positions is 
in a southeasterly direction (Figure 4.1 on page 40). Also the ship’s course is close to 
the azimuth of SV9 (Table 2), which must be the reason for the successive tracking 
losses. The GPS antenna is in the shadow of the ship’s mast when the satellite is low. 

b. Period 2. 19:20 to 19:25 (Cases A and B) 

All positions were computed using data from four satellites. There is a nine- 
second interruption of data in the beginning of the period corresponding to the time 
needed to change tapes in the receiver during logging. 

In Case A the differences range from 17 to 26 m (Table 3). The maximum 
value was found in the first position after the lack of data. Most of the differences are 
around 19 m. In Case B the differences range from 23 to 38 m. The minimum value is 
found after the gap in data. Most of the differences are around 29 m. Visual 
comparison is possible by referring to Figures 4.2 and 4.3. As before the offset of the 


GPS positions is in the southeasterly direction. 


TABLE 5 


| PERIOD 2, CASES A AND B. 
| DIFFERENCES BETWEEN GPS AND MINI-RANGER POSITIONS 


Ranges(m 
vg : 9) маҳ 
роѕ SVs A B A B 
295 6,8,9,11 19 29 17-26 23-38 


our Satellites with geoidal height constraint 


Case A - using four satellites, or 
Case B - geoidal height constraint in Case 
} 


pos - number of consecutive positions computed every second 


c. Period 3. 19:25 to 19:30 ( Cases A and B ) 


TABLE 4 


PERIOD 3, CASES A AND B. 
DIFFERENCES BETWEEN GPS AND MINI-RANGER POSITIONS 


ο. 
Avg Min-Max 
pos Svs A B B 
[52 6,5,9,11 IS 30 17-19 29-3] 
1 6,8.11 18 32 18-18 55 
l 6,8,9,11 13. 2S8 13-13 28-28 * 
9 6,8,11 30 sil 17-37 2 2 
118 6,8,9.11 26 33 24-37 32-36 
11 6,8,1 46 36 31-49 ΙΤ 
4 6,8,9.11 39 57 26-43 21-28 7 
3 ӨЗІ о πῃ. 
* noisy data 


Case A - using four satellites, or | | | 
three satellites with geoidal height constraint 
Case B - geoidal height constraint in Case A 


pos - number of consecutive positions computed every second 


This period is characterized by having noise in different parts of the data 


and shows great variation in the differences. In Case A, at the beginning of the period 


the average differences are 18 m, after the first period of noisv data thev reach values 
οι ο τα, and at the end of the period the differences are up to 46 m. [n Case B at the 
begining of the period the average differences are 30 m. This values changes to 33 m 
after the noisy data, and at end of the period the differences are up to 38 m. It is 
possible to distinguish the oscillation in the positions due to noisy data, much more 
emphasised in Case A than in Case B (Figures 4.4 and 4.5). The offset of the GPS 
positions 15 as before in the southeasterly direction. 
The successive periods lacking data from one satellite are as before due to 
SV9, which is still in an azimuth close to the ship’s course. The fact that use of geoidal 
height as a constraint with data from four satellites for computation of positions does 
not improve the solution is evident in Table 4. 
d. Period 4. 19:30 to 19:35 (Cases A and B) 


TABLE 5 


PERIOD 4, CASES A AND B. 
DIFFERENCES BETWEEN GPS AND MINI-RANGER POSITIONS 





* noisy data 


;Θ 
Avg Min-Max 
pos SVs A B А В 
7 6,9,11 38 38 38-40 38-39 
40 6,3,9.11 3] 35 29-40 34-39 
7 6811 1 38 3243 3739 * 
| 689.1 35 57 38-38 37-37 * 
1 6811 38 58 32838 38-38 * 
49 689.11 30 35 27-358 33-39 
10 68,11 40 37 31-44 35-38 * 
150 ου 11] 3] 3% 25-41 37-40 
| 68,11 34 38 3434 3858” 
Ἱ ο 33 39 2938 36-45 * 
9 6,8,11 41 38 33-46 36-30 * 
2 ο H 39 37 38-40 33S ^ 
I 6,911 38 38 38-38 38.39. κ 
3] 6,8,9,11 38 41 35-44 40-45 
| 


Case A - using four satellites, or | | | 
three satellites with geoidal height constraint 
Case B - geoidal height constraint in Case A 


pos - number of consecutive positions computed every second 





This period is characterized by successive changes in the number of 
satellites used in the computation of positions; the differences are larger than observed 
in the previous periods.Successive losses of lock on SV9 were the cause for this. In 
Case A (Figure 4.6) it is possible to see the oscillations of the computed positions when 
the data are noisy. Improvement in the solutions is seen (Table 5) when we change 
from three to four satellites; the average differences go from 41 m to 30 m, becoming 
higher values when only three satellites are available. 

In Case B the differences have higher values than in Case A, but the ranges 
are smaller due to the smaller, noisy oscillations. As before the noisy Is caused by the 
loss of lock on SV9, and the offset in the GPS positions is in the southeasterly 
direction from the Mini- Ranger positions. 

e. Period 5. 19:35 to 1940 (Cases A and В) 

This period has the noisiest data processed, so the differences have greater 
spread. Although the range between differences is wide, 20 to 68 m in Case A and 29 
to 67 m in Case B, the differences are mostly smaller in Case A than in Case B. 

In Figure 4.7 we see successive jumps from three satellites to four and vice 
versa. The GPS positions are south and east of the Mini-Ranger positions. The 
differences in the northings are almost double of the ones observed in the eastings. 
The reason for the frequent changes in the number of satellites used is the relatively 
low altitude of SV9 and an azimuth close to ship’s heading, which together with the 
pitch and roll led frequently to shadowing of the TI-4100 antenna. 

f. Period 6. 19:47 to 19:53 (Cases A, B and C) 

Ihis period is characterized also by noisy data, but it is for a different 
reason. Here data from SV8 1s rejected during the phase of smoothing it with Doppler 
data. The observed roll is higher than in previous periods, caused perhaps by the 
almost 90° change in course. As Table 2 shows, SV8 is setting, having a lower 
elevation than any other satellite used during this period. 

There are no significant differences between Cases A and B, and average 
differences range from 26 to 30 m. In Case C we see that even lowering the number of 
satellites to two or three and using the geoidal height as a constraint, smaller 
differences are obtained than the ones computed with data from SVI1 (Table 6). The 
offset in the positions has no relation to the course of the ship, re. it is still in the 
southeasterly direction, and a small or no time lag exists between the data of the two 


systems (Figure 4.8). 
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n TABLE 6 


PERIOD 6, CASES A, B AND C. 
DIFFERENCES BETWEEN GPS AND MINI-RANGER POSITIONS 


Ranges(m) 
Avg Min- Max 

pos SVs A B C A B C 

16 89.11.13 ο ο. 09 725220 29-30 28-30 
1 911,13 DMOJNES 07297 29.29 28.28 * 
E 891113 2529 О 82570 29-30 RD e 
EIL PESOS 6 92707 28-28 26-26 * 
4 89.11.13 09 20 "27705 27-30 ολη. 
1 911,13 νο ο 552$ 28-28 26-26 * 

79 591113 27 28 26 26-29 26-10 24-28 
"11; ЖЕ 22 2797 DOE 26-26 * 

ОГИ 29 30 30 27-32 28-33 021-35 
Ж j В ІІ 220 29 27525 29-29 28-28 3 

Ius o FINIS 29 29 29 28.30 28-30 D730 
| 911 13 ЖА 20: 27 28-25 29-29 2e E 

PONES OST IIIS 29 20 28 28-31 2831 26-2 
ο P280 27 δε 29.29 oS 

162 89.11.13 2809. 27. 26-3] ТЕЗІ 24-31 
| І 911,13 29 30 27 2524 30-30 ιτ.” 
| 1 89.11.13 29 30 28 29-29 30-30 28-28 * 
| |. SE ЖӘ 7 772 29-29 27-271 
| ME IIS 29 30 28 29-29 30-30 29.24 
| l 911,13 ὄπ 777200 28-28 pops 
5 89.11.13 МЕ 27:27 76-27 26-27 20-22 ық 
ο ае 27 1... 27.27 22 21-21 # 
ο ο]: ο ο ο). 00260 29.3] 20-2] * 
1 91113 5 00-23: 2076 327 17 2ο... 

6 8911.13 ος δα οἱ 26205 VIDE 20-24 


8 


* noisy data 


Case A - using four satellites, or | : . 
three satellites with geoidal height constraint 
Case B - geoidal height constraint in Case A 
Case C - ignoring SV11 but constraining 
with geoidal height. . 
pos - number of consecutive positions computed every second 





g. Period 7. 19:57 to 20:05 (Cases A, B and C) 

In Case B there is an improvement in the positions in relation to Case A. If 
we refer to previous periods, we see that what is happening is the reverse of what was 
happening there, i.e. the use of the geoidal height in positions computed using data 
from four satellites improves the solution. As before, we see that noisy data causes an 
oscillation in the observed differences but with small amplitudes. Noteworthy are the 


highly improved positions indicated by a 13-m difference (Table 7). 
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- TABLE 7 


PERIOD 7, CASES A, B AND C 
DIFFERENCES BETWEEN GPS AND MINI-RANGER POSITIONS 


Ranges(m) : 
Avg Міп- Мах 

pos SVs А ВСС А В С 
188 οι. 28823 701 2020 19-25 [7-25 
E UNIS 2482535 τοῦ ED 4-24 23228 20-20 3 
22 91l 209850 Ws EDO-24 |8520 16-20 
1 811,13 2pm AMD 7222 21-24 18-18 * 
ΠΠ “οἱ 228801 NS 20-25 19-23 | 

JP S DHEN P E, 22897 19-19 * 
12 38.9 (ie ο. 10 MES οι ο 19-21 17-2] 

М ү [9 О Чо 19-[9 19-19 16-16 3 
e 20 18 16 18-22 17-19 13-18 
Les its 19 19 16 19-19 19-19 16-16 * 
13 οἱ) [oes Ji [7-20 LZ 13-16 
ІШІ ӘРІ? 17 Б 12-085 14 * 
З ов 21 19 16 19-26 17-24 14-21 
ОО | РГ ЛГ, К-ү АЛО 14-14 * 
35 ИЕ 20 18 16 19-22 16-20 13-18 
(2 ӘНГЕ Кас δε 16-16 13-13 * 
2 2801] I3 29 23 23 28-29 25.26 22.2) * 
І 911,13 27 24 210 δη ο λα. ο... 
27 οτι 25 23 21] 18-36 17:34 15-31 


* noisy data 
Case A - using four satellites, or | | | 
three satellites with geoidal height constraint 
Case B - geoidal height constraint in Case 
ase C - SERE SV11 but constraining 
with geoidal height. . 
pos - number of consecutive positions computed everv second 


The GPS positions are also south and east of the MR Falcon positions, but 

the fact of ignoring SV11 causes a decrease in the differences in easting (Figure 4.8). 
h. Period 8. 20:07 to 20:13 (Cases A, B and C) 

During this period all differences in all cases are smaller than 21 m (Table 
8) The higher values are obtained after the computation of a position using fewer 
than four satellites in Cases A and B, and three satellites in Case C, and are due to the 
oscillation caused by the Kalman filter. Table 8 shows that constraining the positions 
computed using data from four satellites improves the solution; the minimum difference 
reaches only 10 m. Also, ignoring SV11 causes a reduction of the differences, mainly in 
the easting, and we see that the differences are smaller than in the other cases (Table 
8). Figures 4.10 and 4.11 show the improvement in the positions when SVII is 


Ignored. 
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TABLE Š 


PERIOD 8, CASES A, B AND C. 
DIFFERENCES BETWEEN GPS AND MINI-RANGER POSITIONS 







Ranges(m) | 
Avg Min- Max 
DOS SVs A B C A B C 
37 S,9.11,13 19 19 14 19-21 16-18 2515 
ПЕ И 18 17 14 18-18 17-17 14-14 * 
153 89.11.13 ο μή. 17-21 16-19 ο] 
] ΠΠ; ШӘЛІ 1217 |7417] μι. 
28 89.11.13 ӘШІШТІЗ 191 16-19 13-16 
| Ш g RDIS ИО 17-17 19517 1, 
142 89,11,13 ο — 17-20 15-17 10-14 
* noisy data 


Case A - using four satellites, or | | : 
three satellites with geoidal height constraint 
Case B - geoidal height constraint in Case 
Case C - ignoring SV11 but constraining 
with geoidal height. . 
pos - number of consecutive positions computed every second 






i. Period 9. 20:13 to 20:20 (Cases A, B and С) 


ΙΕ» 


PERIOD 9 CASESA B AND C. 
DIFFERENCES BETWEEN GPS AND MINI-RANGER POSITIONS 


Ranges(m) 


Avg Min-Max 
pos SVs B C A B e 
ГЭЕ 8:9. 11.13 20 17 13 19-22 15-18 11-14 
1 no data 18 15 11 18-15 15-15 пеш 5 
Sp 691113 ΤΙ η 15-16 10-11 


* noisy data 


Case A - using four satellites, or | | 
three satellites with geoidal height constraint 
Case B - geoidal height constraint in Case A 
Case C - ignoring SV11 but constraining 
with geoidal height , κ 
pos - number of consecutive positions computed every second 





3 


The differences found during this period are similar to the ones obtained in 
period 8. Table 9 Indicates that the differences in all cases are correlated as before. The 
oscillations are due to the noise in the data from the MR Falcon and not to the GPS 
positions, which define a smooth and continuous course (Figure 4.12). 

j. Period 10. 20:20 to 20:25 (Cases A, B and O) 


TABĽCE 1 


PERIOD 160, CASES A, B AND C. 
DIFFERENCES BETWEEN GPS AND MINI-RANGER POSITIONS 


Ranges(m) | 
АУЕ Min-Max 
pos SVs A B C A B C 
Ә-Ә ШЕ 20 ο η ДІ 14-16 9-11 
QS 5ο. ll Пе 10-17 10-15 10-12 
206 οἱ 23 25 24 15-26 15-26 19-26 


Case À - using four satellites, or | | | 
three satellites with geoidal height constraint 
Case B - geoidal height constraint in Case 
Case C - ignoring SV11 but constraining 
with geoidal height, . 
pos - number of consecutive positions computed every second 


During the first tests, the differences were smaller when the positions where 
computed from three satellites, SV 11 excluded and the ship constrained to the geoidal 
height, than the ones using four satellites, SV11 or SV6 included. This happened 
because during data acquisition SV6 was selected to replace SV11, and while the 
receiver was not locked on the SV6, only data collected from three sateilites was good 
(Figure 4.13). 

In Case B (Figure 4.14 or Table 10) constraint of the position to the 
geoidal height improves the solution where SV11 is used. In Case C (Figure 4.15) we 
notice a reduction in the differences in the positions where data from SV11 was used. 
When SV6 replaces SV11, in all cases the differences are larger but of same 


approximate magnitude. 


k. Period 11. 20:30 to 20:37 (Cases A, B and C) 


~ 
= 


TABLE 11 


PERIOD 11, CASES A, B AND C. 
DIFFERENCES BETWEEN GPS AND MINI-RANGER POSITIONS 


Ау Min-Max 
pos SVs А B е А В С 
130 6,8,9 13 ὃν ος 35 5536 5ο 37-39 
287 8913 36 36 36 30.37 30-38 30-38 


$ 


Case A - using four satellites, or | . ' 
three satellites with geoidal height constraint 
Case B - geoidal height constraint in Case А 
Case C - ignoring SV11 but constraining 
with geoidal height, . 
pos - number of consecutive positions computed every second 


Ranges(m) 


The initial positions were computed using data from four satellites, and 
after 20h32m10.7s the receiver started to track SVI2 in place of SV8, but no broadcast 
ephemeris data were received during this period. The differences computed using data 
from three satellites are the same in all cases, but those computed using data from four 
satellites have different values. Thus in Case A the differences are smaller than in Cases 
B and C, where the positions are constrained to the geoidal height (Table 11). па! 
cases the differences in eastings are larger than in northings, contrary to what was 


happening before; this offset is similar for the other two cases (Figure 4.16). 
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Figure d.1 Period I (Case A). Using all available satellites. 
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Figure 4.2 Period 2 (Case A). Using all available satellites. 
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Figure 4.3 Period 2 (Case D). Constraining to geoidal height. 
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Figure 4.4 Period 3 (Case A). Using all available satellites. 


43 


SERFLOOR BENCHMRRK EXPERIMENT 
PHASE [I 
16 AUGUST 1986 


1930 





N 
Satellites used: 
SV6, SV8, 5У9, 5У11 1325 
1925 
+-MR FALCON TRANOUCER POSITIONS 
X GPS TRRNDUCER POSITIONS 
S | 290 2 50 


METERS 


Figure 4.5 Period 3 (Case DB). Constraining to geoidal height. 
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Figure 4.6 Period 4 (Case A). Using all avallable satellites. 
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Figure 4.7 Period 5 (Case А). Using all available satellites. 
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Figure 4.8 Period 6 (Case C). Ignoring SV11. 
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Figure 4. 9 Period 7 (Case B). Constraining to geoidal height. 
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Figure 4.10 Period 8 (Case B). Constraining to geoidal height. 
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Figure 4.11 Period 8 (Case C). Ignoring SVII. 
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Figure 4.12 Period 9 (Case C). Ignoring SV11. 
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Figure 4.13 Period 10 (Case A).. Using all available satellites. 
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Figure 4.14 Period [0 (Case B). Constraining to geoidal height. 
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Figure 4.15 Period 10 (Case C). Ignoring SV11. 
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Figure 4.16 Period 11 (Case A). Using all available satellites. 
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4. Programs. Inputs and outputs 


Two programs written for data analysis are summarized below: 


Program: COMPARE POSITION 
Input: Control file with period limits (Appendix K) 
File with MR Falcon positions 
File with GPS positions 
Output: Listing with date, time tag, and GPS geographic positions 
and differences in UTM coordinates 


Source: See Appendix G 


‘Program: COMPARE PLOT 
Input: Control file with period limits (Appendix K) 
File with MR Falcon positions 
File with GPS positions 
Output: Plots of tracks using GPS and MR Falcon data 
Source: See Appendix H 


B. CONCLUSIONS 

The results of the data analvsis using four satellites with the broadcast ephemeris 
indicate that the real-time positioning by GPS of a dynamic platform, e.g. a ship, under 
the best conditions is within 15 m but in most cases 20 m. 

Using data from three satellites and constraining the solution to the geoidal 
height, positional accuracies are within 20 m under the best conditions and 30 m in 
most of the cases. 

[t was found that the GPS positions were always southeasterly of the Mini- 
Ranger positions. Due to lack of sufficient data and due to time limitations, correlation 
between the geometry of the satellites and ship's heading could not be investigated. 

The accuracy of the positions is correlated with the relative motion of the 
satellites; this was found when SV11 was ignored as in Case C, which improved the 
solution. If the variations in azimuth and elevation of all satellites are compared for all 
periods the data was processed, SVs 6 and 11 have variations larger than 110 degrees, 


while SVs 8, 9 and 13 have variations in azimuth smaller than 35 degrees (Table 2). It 
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can be seen that both SVs 6 and 11 reach their culmination during the period the data 
is analyzed. SV6 reaches its culmination around 20:15 and SV11 around 19:50. 

When SV11 is approaching culmination, the solution is not improved when the 
position 1s constrained to the geoid, after culmination however, application of the 
constraint improves the solution (Tables 12, 13 and 14). 

The degradation in the accuracy of the GPS positions using either SV6 or SVII 
is evident in Table 15 where the differences found in these situations are bigger than for 
positions computed without data from them. It is evident that even in Case C when 
SV6, with the same characteristics as SV11, is used to compute the position, the 
differences become larger. 

There are some theoretical studies (Landau, 1986) about the selection of satellite 
configuration in order to get the best results for computed positions. However, with 
the present data it is not possible to select a best combination, since the data 
acquisition was limited to a small observational period and a small number of satellites. 
Currently a way of measuring the effect of the geometry of the satellite configuration is 
through the analysis of the geometric dilution of precision (GDOP) parameters. These 
parameters include the position dilution of precision (PDOP), reflecting the dilution of 
precision in three dimensions, the the horizontal dilution of precision (HDOP), 
reflecting the dilution of precision in two dimensions, the the vertical dilution of 
precision (VDOP), reflecting the dilution of precision in the vertical dimension, and the 
time dilution of precision (TDOP), reflecting the dilution of precision in time (Milliken 
and Zoller, 1980). Thus, a low PDOP provides a good geometric configuration. Table 
16 lists the observed satellites and their GDOP parameters. 

It was found that the position of the antenna is important, and in the case of the 
NSWC TI-4100 receiver some error was caused by shadowing due to the ship's mast 
when the azimuth of the satellites was close to the course of the ship. 

The positions used for comparison with the GPS positions are themselves 
affected by the noisy signal of the MR Falcon positioning system. This problem could 


be reduced if more than two LOPs were used in the computation of positions. 


M 
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TABLE 12 


MR FALCON POSITIONS VERSUS GPS POSITIONS. 
OBSERVED DIFFERENCES IN METERS DURING PERIOD: 


H в px DY pp рх ҚҰСЫ рх GS pP 
9550; в 15-4 52-2- 
9.2557 M ANN ж авва ....- 
190 25 10 40 DD DD -- 
955 ІІІ οἱ ы 
РІ τι... 
1025 057 с 0 M NN n 
ο 25ο 5. 2: Б 
2. 

ІІ ооо обр n 2 


19 256506 10 -16 18 13 -26 30 -- -- -- 


ОХ, DY differences in the UTM coordinates 
DP derived position difference from DX and DY 


Case A - using four satellites, or р | 
three satellites with geoidal height constraint 
Case B - geoidal height constraint in Case 
ase C - no data during this period 


x 
| 
-16 19 12 -27 30 -- -- -- 








ΕΕ» 


MR FALCONPOSITIOSS VERSUS GPS POSITIONS. 
OBSERVEDIDIFFERENCES ISN METERS DURING PERIOD 6 


Time Case A Case B Case C 
H M S DA 5 DY DP DXT DY DP DX DY DP 
| 19 52 407 17 :20 25 19 -22 29 16 -21 26 
| 19 52 417 16 -24 28 18 -23 29 15 -22 27 
| 19 52 427 16 -24 29 18 -23 30 16 -22 27 
| 19 52 437 17 -22 28 19 -23 30 16 -22 27 
| 19 52 447 15 -24 29 16 -25 30 14 -24 28 
19 52 457 16 -21 27 18 "227 72457 715 -2] 26 
19 52 46.7 10 -27 29 12 -28 30 8 -27 28 
19 52 477 14 -23 27 15 -24 25 8 20008 723 
19 52 48.7 14 -23 27 15 2235 228 9 -20 22 
19 52 49.7 15 -22 26 16 -21 26 10 -18 21 
БЛ -507- 17 -20 26 18 -19 26 12 -16 20 
| ο ο 215 ὄν. 7 19 -19 27 14 -17 2l 
БЕ 2277 19 -19 27 20 -19 27 15 | 22 
Ie οἱ 19 -19 27 20 -19 28 16 -18 24 
19 52 547 13 -27 30 13 -27 31 10 -26 28 
ID 557 |І) -26 29 13 -26 29 9 -25 26 


DX, DY differences in the UTM coordinates 
DP derived position difference from DX and DY 


Case A - using four satellites, or | | 
three satellites with geoidal height constraint 

Case B - geoidal height constraint in Case / 

Case C - ignoring SV11 but constraining with 
geoidal height 








TABLE 14 


MR FALCON POSITIONS VERSUS GPS POSITIONS. 
OBSERVED DIFFERENCES IN METERS DURING PERIOD 8 


| Time Case А Case B Case C 
ма: DX DY DP DX DY DP DX DY DP 
20 7 οὗ Ар 9. - Ἡν ου М 
20 7 207 по 0) ὃ -15 ІП i 
Ὁ τ. ΠΠ... 9 1217 2884 Б 
η." .:. . 3 7 С .αϱ ПО 
| 20 8 20.7 10 -17 20 8 С ІС 
20 8 407 БП :- ΙΙΙ 
| 20 9 07 Πα ө, 5 ІШПЕ ТЕ 
| 208 9-120 2 ΠΠ оо рота Sp 444 12 
| 208. 9 1240.7 πι ΠΠ 8 -1590 17 WS Pole 
20— 10 07 ШІ STO 2 447167820 ЕБ 
20 10 207 9 -17 19 87 5-16 7 MES. (ВЕ 
20: 10 ОТТІ ПП o2 o EP NT? 
| 21: 97 buc 3» 8 215 NUT NEA ιο Ti 
ОРОС 9 "4 NES MEI» 
20 |. ον 2 πλ 8 -4 16777 тошу 
20:12 7707 νο ШИШ: OESTE NI 


DP derived position difference from DX and DY 


Case A - using four satellites, or | | | 
three satellites with geoidal height constraint 

Case B - geoidal height constraint in Case ; 

Case C - ignoring SV11 but constraining with 
geoidal height 


x DX, DY differences in the UTM coordinates 
| 
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OBSERVED DIFFERENCES IN METERS DURING PERIOD 10 


HAs px ҰСЫ рх Ч рр px DY DP 
0 20 77 14 455 2009013 15 а 10 и 
0 20 87 14 45 20 9 4315 42 -0 1 
НИ oT lle ον 17 5 2715 а 00 1 
ІІ ο πμ 1 02,2 1, 
2 ТЕМ 10 1 (2 0 π 
ІІ КООШ у Б зз їп 
20 20 447 а JO 10 4. шю 10 4 -0 10 
20 20 467 4. 10 10 4. -ἰ 10 4 -0 10 
СИ СОС P u 


PONE] “557 B 10041 -10 ll 3 -10 10 
ДЕТІ 56.7 3 -10 ГТ -10 11 3 -10 ll 
0 МИ 37.7 8 108 15 -13 15 Š -13 15 
20 21 387 10 -14 18 10 -14 18 10 -14 18 
ΓΙ ΟΥ, 5; 12 -15 19 I» -15 720 =I? -15 20 


τ 40722215 -16 21 13 -16 20 13 -16 21 


со „э ο . 


TABLE ТӘ 
MR FALCONROSITIONS VERSUS GPS POSITIONS. 
| 


ОХ, DY differences in the UTM coordinates 
x DP derived position difference from DX and DY 


Case A - using four satellites, or 
three satellites with geoidal ΠΕΙ constraint 
Case B - geoidal height constraint in Case À 
Case C - ignoring SV11 but constraining with 
geoidal height 


61 


Time 
19:00 
19:20 
19:40 
20:20 


СУ 5ο СӨ (олу СА 


TABLE 16 


PREDICTED GDOP PARAMETERS 
l6 AUGUST 1986 


SV set GDOP PDOP HDOP YDOP TDOP 
S 6.90 55 2.47 3.19 5.82 
ΙΙ 7.61 6.30 2.05 5 4.47 


9 11 13 3:97 3.01 1.97 221 1:55 
ΙΙΙ; 3.56 555 215 2.22 22 
ο ου 3:38 2.94 1.88 2-75 1.58 
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V. RECOMMENDATIONS 


GPS data from other days of the Phase II cruise should be used with the 
modified programs discussed here to further check our conclusions. 

lests should be made in order to find the reason for offsets in a southeasterly 
direction and evaluate their correlation, if any, with the geometry of the satellites used. 

Io test the conclusion that the accuracy depends on the large variations of the 
azimuth of the satellites, SVs other than SV6 and SV11 should be used. 

The relation between satellites’ culminations during observations and the 
accuracy of the positions should be further investigated to provide a criterion for 
satellite selection. 

A better antenna site on the ship should be found to avoid shadowing and to 
minimize the effects of pitch and roll. Antennas of two similar systems should be 
installed at both high and low elevations to determine the effects of pitch and roll. 

Ihe Kalman filter should be improved to avoid the oscillation when data are 
noisy. Other algorithms should be implemented and their results compared to the 
results of the KALMN2 program. 

To provide a real-time position computation the KALMN2 program or a similar 
program should be installed in a transportable computer. A system should be designed 
to have one computer converting and logging the data, another processing and 
displaving results, saving them or sending them to still another computer or data 
recording system, and there should be a master system to control the synchronization 
of time recorded with all data and to control the flow of information. 

Processing of the same data should be done using the precise ephemeris to see 
how the solution is improved. 

Processing in differential mode should be done to see how much the solution 1s 
improved. 

Comparison of data from other receivers should be done to check which type of 
receiver provides higher accuracy. 

Next an experiment should be designed to obtain point or differential positions in 


real time for moving platforms (survey ship) with little or no shore support. 
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To allow better evaluation of the GPS data. positions 


computed with more than 
two lines of position should be used as reference. 
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APPENDIX À 
PROGRANI FALCON. SOURCE LISTING 


//EZEQUIEL JOB (0812,9999),'EZEQUIEL',;,CLASS-C 

/ /'*MAIN ORGzNPGVMI.0812P,LINES-(99),CARDS-(99) 
//*FORMAT PR,DDNAME -GO.FTO6FQO01, 

//*FORMS=SEP1 

// EXEC FORTVCLG 

//FORT.SYSIN DD * 


AUTHOR: AUGUSTO EZEQUIEL 
DATE: APRIL 10, 1987 


DESCRIPTION: 

THE PROGRAM TAKES THE RAW RANGE-RANGE DATA FROM AN INPUT FILE; 
COMPUTS THE POSITIONS AND AND THE COURSE BETWEEN TWO ANTENNA 
POSITIONS AND PRINTS THE UTM COORDINATES, THE COMPUTED COURSE 
AND THE GPS TIME OF EACH OBSERVATION. 

THE PROGRAM RUNS IN MVS 

I/O SPECIFICATIONS: SEE END OF THIS JOB. 


THIS PROGRAM IS RESTRICT AS IS TO THE SEAFLOOR BENCHMARK EXPERIMENT 
FOR OTHER USES DO NOT FORGET TO ENTER THE COORDINATES OF THE 
STATIONS AND OTHER PARAMETERS. 

COORDINATES ARE IN WGS72 DATUM. 


THE PROGRAM WILL RUN WITH ANY AMOUNT OF DATA ONLY LIMITED TO DISK 
SPACE 


ANY BLANK LINES WILL TERMINATE THE PROGRAM IN ERROR 


соо о: осоо ооо оо ооо оо осо о ооо о о) 


REAL ж8 X(2),Y(2),XPOS,YPOS,XOLD,YOLD 

INTEGER LADEG1,LADEG2 ,LAMIN1,;LAMIN2 ,LODEG1 , LODEG2 ,LOMHIN1 ,LOMIN2 
INTEGER MONTH,YEAR,DAY ;HOUR ;MIN;LADEG ; LAMIN,;LODEG LOMIN 

REAL LASEC1,LASEC2,LOSEC1,LOSEC2,LASEG ,LOSEG,SECS 


STATION 1 FERRIER 


ооо 


LADEG1=36 
LAMIN1=33 
LASEC1=53. 748 
LODEG1=121 
LOMIN1=53 
LOSEC1=52.939 


STATION 2 DOME ECC 


ооо 


LADEG2=36 
LAMIN2=18 
LASEC2=20.95 
LODEG2=121 
LOMIN2-54 
LOSEC2=00.50 


COMPUTE X-Y FOR CONTROL STATIONS 


оо о 


CALL GPUTMILADEG1,LAMIN1,LASEC1,LODEG1,LOMIN1,LOSEC1,Y(1) ,X(1)) 
CALL GPUTMILADEG2,LAMIN2,LASEC2 ,LODEG2 ,LOMIN2 ,LOSEC2 ,Y(2) ,Xt2)) 


05 


END OF INITIAL DATA 
COMPUTS THE FIRST POSITION TO HAVE A STARTING POINT FOR 
EVALUATION OF THE COURSE 
READ( 50, * ,END=100) ΜΟΝΤΗ ,ΟΑΥ YEAR, ,HOUR MIN ,)SECS, 
w KODE1,RANG1,;ISTRE1,KODE2 ,RANG2 ,ISTRE2 
CALL RRXY( X,Y »RANG1,RANGZ,XOLD,YOLD ) 


READ IN DATA FROM INPUT FILE 


ооо 


50 CONTINUE 
ВЕАО( 50, ж ΕΝΟΞΙΟΟ) MONTH,DAY,YEAR,HOUR,MIN,SECS, 
x KODELI,RANGI,ISTREl1,KODE2;,RANG2,ISTRE2 


COMPUTE X-Y AND GP OF RR POSITION 


ооо 


CALL RRXY(CX,Y ;RANG1;RANG2,XPOS,YPOS) 


COMPUTES THE COURSE 


ооо 


CALL HEAD( XOLO,YOLO,XPOS,YPOS,COURSE »>IERR ) 
IF LANE JUMPS DOES NOT SAVE THE POSITION 
IF(IERR.LT.O) WRITE(6,1)IERR 


1 FORMAT(I6) 
IF(IERR.LT.0) GOTO 5D 


O O O O oO 


COMPUTES THE GEOGRAPHIC COORDINATES 


ооо 


CALL UTMGPt( YPOS ,XPOS ;LADEG; LAMIN;,LASEG,LODEG ; LOMIN,LOSEG) 


OUTPUT THE RESULTS 


ооо 


WRITE(6,60) MONTH,DAY , YEAR; HOUR, MIN,SECS,LADEG,LAMIN,LASEG 
# » LODEG , LOMIN , LOSEG , COURSE 
60 FORMAT(2(1X,I2),1X,I4,3(1X,I3,1X,I2,1x,F7.3),1X,F5.0) 


5 
НКТТЕ( 7) МОМТН , ОДУ ,УЕДК , НОЦК „МТМ, 5ЕС5 , ГАОЕС , |АМТНМ , |АЗЕС 
+ » LODEG , LOMIN , LOSEG 
C 
C SAVES THE COURSE IN A SEPARATE FILE 
Е 
WRITE(8) MONTH,DAY,YEAR,HOUR »>MIN »SECS ,COURSE 
Е 
C SAVES POSITION FOR COURSE COMPUTATION 
C 
XOLD-XPOS 
YOLD=YPOS 
GOTO 50 
Е 
C END OF PROGRAM 
C 
100 CONTINUE 
STOP 
END 


C 
C3XGO00090900000000000000)0090)00OODXDOODODOO)0(OODOQOODO00090!01)D0O0GOOO)9GOQODO9OOOO)OO0)0O00O00]00!0€ 
SUBROUTINE RRXY(XR,YR,R1,R2,XCO,YCO) 

C3O(OOCOODO00000)00G)000Q0O0O0O000XIDDGODUODQOOOO0O00C JE JEJE HE JE JE JE IE FE HJE IE FE HE EJE JEJE FE FE EIE EH JEFE HME 
C 

IMPLICIT REAL*8 (А-Н,0-2) 

DIMENSION XR(2),YR(2) 

A = DOSQRT((XR(2)-XR(1))%%2 + (YR(2)-YR(1))%%2) 


CR = (XR(2)-XR(1)) Z À 
SR = (YR(2)-YR(1)) Z À 
ΧΡ = (R1%R1 - R2wR2 + AwA) Z (2%A) 


ARG = R1¥R1 - XP*XP 
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C 


YP - OSQRT(ARG) 


ХСО - XPXCR + YP%SR * XR(1) 
YCO = XP%SR - ҮРЖСВ * YR(1) 
RETURN 

END 


C EE 3€ X EO 99EXC EXE OE E CX XE CC CC COE OE C € EXE OXEXOOXE EXE JE EIEE JEJE JE JE % 34334 IE IESE IE IEEE HE IEE 


SUBROUTINE GPUTM(LAOEG,LAMIN,LASEG,LOOEG,LOMIN,LOSEG ,NORTH>EAST ) 


С FEE SE IE 3€3€ € E€993€ 3€ € € 3€3€3E HE JE FEE SESE CC HE JE JE JE JE JE JE JE EJE JE JE JE JE JEJE JE JE JEJE JE JE JE JE JE JE IE IE JE JE 555556565565 KE 


C 


осоо со) 


OOUBLE РВЕСТ15ТОМ А,Н›,М,АР,ВР,СР,ОР,ЕР,5,В1›Е50@›Е5$@Р›ВМ,ВР›,КО 
OOUBLE PRECISION R2,R3,R6,R5,P,P2,P3,P6,P5,P6,A6,B5,SINSEC 
DOUBLE PRECISION PHI;,DLAM,NORTH,;EAST,;B,PHIMIN,;PI,LON 

REAL LASEG,LOSEG 


THIS SUBROUTINE COMPUTS THE UTM COOROINATES OF GP IN WGS 72 


IN 


$ 


$ 


$ 


ZONE 10 CENTRAL MERIOIAN 123 00 00 W 

CM=-123.0D0 

PHI =OFLOAT( LADEG )+DFLOAT(LAMIN )/60.00+0BLE( LASEG )/3600.D0 
LON=DFLOAT( LOOEG )+OFLOAT( LOMIN )/60.00+0BLE( LOSEG )/3600.D0 
LON=-LON 

DLAM=( LON-CM )*3600.000 


A=6378135.000 
R=298. 2600 


KO=0.9996D0 
B - AX(R-1.D0)/R 


М = (А-В )/(А+В ) 


AP * AX((1.00-N)*5.00/4.00*( N92 -N3XX*3 ) 81. 00/66 . OOX (N30€& -N3*XE5 ) ) 

BP = $3.D0/2.00X*A*( (N-NXx2)«7.D0/8.D0*( Nx»x3 -Nxsxà ) 455. D0/66 . DOXNXX5 ) 
CP z 15.00/16.00XA*( NX*2-Nx*3«3.00/4 . 00X( NXx& -NXX5) ) 

2р = 35.D0/48.00*xAX( Nx*3-N39*4«11.00/16.00*XN395) 

ЕР = 315.00/512.00X*A*( Nxx«4 -NX*5 ) 


PHIMIN Ξξ PHI*60.00*2.90888208666D-4 

PI-OARCOS(-1.00) 

PHI-PHI/180.00xPI 

S - APXPHIMIN-BPXDSIN( 2. 00XPHI )*CPXOSIN(4.OOXPHI) 
- DPXxOSIN(6.00*PHI )*EPXDSIN(8.DOXPHI) 

R1 - KOx*S 


ЅІМЅЕС = (1.00/3600.00 )/180.00%РІ 
SINSEC=OSIN(SINSEC) 


ESQ < (Джжд2-Взее2 )/Джжд 
ESQP - Е50/(1.00-Е50) 
RM 3» AX(1.00-ESQ)/( DSQRT( 1.00 -ESQXOSIN(PHI )32 ) x3 


RP = RM*(1.00+ESQP*DCOS( PHI )**2 ) 

R2 = RPXOSIN( PHI )*OCOS( PHI )*SINSEC%*2/2 .00*KO*1.008 

RS s SINSECX*axRPxOSIN(PHI )*DCOS( PHI )*3/26 . 00*( 5. 00 -OTAN(PHI )*ez 
+ 9, 00%ESQP*DCOS( PHI )**2+4. OOXESQP*ESQP*DCOS( PHI )**4 )*KOX1.016 


RG = RPXOCOS( PHI )*SINSECXKO1 . 04 
R5 = SINSEC**3*RPxDCOS( PHI )**3/6.D0%(1.00-OTAN( PHT }x#*2 
* ESQPXOCOS(PHI )**2 )*KOx1.012 


P - .000100Xx0LAM 


рг = ржж2 
PZ = PxxZ 
рё = ржжё 
р5 = ржж5 
P6 - P**6 


A6 * P6XSINSEC*6XRPX*DSIN(PHI )*OCOS(PHI )x*5/720.D0 
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5 x (61.00-58.DO0XxDTAN( PHI )*3:t2 ЖОТАМ(РНІ ) еже 
$ + 270.D0*ESQP*DCOS( PHI )**2-330.D0*ESQP¥DSIN( PHI )€*2)xK0x1.D24 
BS = P5*SINSEC**5*RP*DCOS( PHI )**5/120.00*(5.00-18.DOXDTAN( PHI 122 
$ + DTAN( PHI )**44+14.D0*ESQP*DCOS( PHI }**2 
$ ~ 58.DO0xESQPXDSIN( PHI )*x«2 )xKOx1.0D20 
МОВТНЕВ 1 + К2жр2 + рЗжра+д6 
EAST=(RG*P+R5*¥P3+B5 }+500000.00 
RETURN 
END 
C 
C ARAE EEIE EIE MEIE EE JE MEIE EIEI JE IEEE IE ME JE NEEE ME MEIE EHIE ME ENEE ME ENEJ JEJE EME IE JE ME IE E JE ME EM AEAEE EE 
SUBROUTINE HEAD( XOLD,YOLD,XPOS,YPOS,COURSE,IERR } 
CE segs Oeo IE IRE ecc EO e je С зе e SCC e SE ο ο Ж RE 
C 
ΒΕΕΑΙ «8 XOLD,YOLD,XPOS»,YPOS 
IERR=0 
DXzXPOS-XOLD 
DY=YPOS-YOLD 
RADIUSZSQRT ( DX*x2 «DY3€x2 ) 
IF(RADIUS.GT.20.) IERR=-1 
IF(RADIUS.GT.0.) GO TO 10 
COURSE 20.0 
RETURN 
10 CONTINUE 
20 CCURSE= ARCOS(DY/RADIUS } 
COURSE =COURSE*180 ./3.141592654 
IF(OX.LT.O. } COURSE=360.-COURSE 
RETURN 
END 
C 
Cot ope x CL EM AGE E IE JE IE ME XC ЕЕ Ж 555 ж ЖКЖ ο. 
SUBROUTINE UTMGP(NORTH ; EAST »,LADEG,LAMIN,LASEG,LODEG,LOMIN,LOSEG } 
Ооа επ e doo mo e e e ac Mee o€ HE IE SE Se 060€ € 3€ 96 3€ 9€ 9€ 9€ 5€ X6 3€ € 0 9€ € C ХХХ ХХХ ХХ ^ 9€ € 3€ 
C 
DOUBLE PRECISION A,R,N,B, BP,CP,DP,EP,S,R1,ESQ,ESQP,RM;RP,KO 
DOUBLE PRECISION R7,R8>R9,E5,Q,Q2,Q3 ,Q4,Q5,Q6 ,06,AP,SINSEC 
DOUBLE PRECISION R10,DLAM,NORTH,EAST,RPHI,DOPHI »PHIMIN 
DOUBLE PRECISION EPRIME ,DELTA,D6NUM,CM,PHI ,LON,PI 
REAL LASEG,LOSEG 


С 
C THIS SUBROUTINE COMPUTS THE UTM COORDINATES OF GP IN WGS 72 
C IN ZONE 10 CENTRAL MERIDIAN 123 00 00 И 
С 
СМ=-125.000 
Є 
A-6378135.0D0 
R=298.26D0 
C 
KO=0.9996D0 
C 
PI-DARCOS(-1.D0) 
B - AX(R-1.D0)/R 
C 
М = (А-В)/(А+В) 
AP * A*((1.D00-N)*5.D0/4.D0*( Nx*2-Nx«3)*81.D0/64 . DOX( NYX& -N3x 5 ) ) 
BP = 3.00/2.D0*A*( (N-Nx*«2)*7.D0/8.DO0*( NXx3-N3X*& ) 55. DO/66 . DOXNXXB5 ) 
CP = 15.00/16. DOXA ( Ne 2 -N%**3+3 00/4. DO*( N¥*G-N¥5 ) } 
DP = 35.D0/48.00X*AX( Nxx3-N3*3x*4 «11.D0/16.D0XNX*5) 
EP = 315.D0/512.D0*A*( NX*& -N**5 ) 
C 
C FIRST APROXIMATION OF PHI 
C 
PHI =NORTH/30.800/3600.00 
C 
C COMPUT TRUE MERIDIONAL DIST AND APROXIMATE PHI 
E 


DO 100 I=1,8 
PHIMIN=PHI*60 .D0*2 . 90888208666D-4% 
ΒΡΗΙΞΡΗΙ/180.ΌΟΧΡΙ 
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S = APXPHIMIN-BP*DSIN( 2. DO*RPHI )*CPx*DSIN( 4. DOxRPHI ) 
$ -DP*DSIN( 6. DO«xRPHI ) *EPXDSIN(8.DOXxRPHI) 
Rl = KOxS 
DELTAzNORTH-R1 
PHI=(DELTA/30.8/3600. )+PHI 
100 CONTINUE 
RPHIZPHI/180.DO0xPI 


SINSEC 
SINSEC 


(1.00/560000?)/180.00%РІ 
DSIN( SINSEC) 


ESQ с (Дж%Ұ2-Вжже2г )/A*x2 
ESQP - ESQ/(1.D0-ESQ) 
RM * AX*(1.DO0-ESQ)/( DSQRT( 1.00-ESQXDSIN(RPHI )9*2 ) ) 3 


RP - RM*(1.DO0*ESQPXDCOS(RPHI )s*2) 

R7z(DTAN(CRPHI )/( 2. DO«RP3&2*SINSEC ) )*€( 1. DO«ESQPX«DCOS( RPHI )*x2 ) 
$ *¥1.012/KOx*2 

R8=(DTAN(RPHI )/( 24. DO*RP¥¥4G*SINSEC ) )*(5.00+3.00%DTAN( RPHI )**2 
$ +6 .DOXESQP*DCOS( RPHI )**2-6 .DOxESQPXDSIN( RPHI }**¥2-3 . DOXESQP¥*2 
$ XDCOS ( RPHI )x4&-9.DOX*ESQP3*2x*DCOS( RPHI )382xDSIN(RPHI )x2 ) 
$ #1.024/КОжжа4 

R9z1.D0/DCOS( RPHI )/( RP*XSINSEC )*1.D6/K0 

R10-1.D0/DCOS( RPHI)/(6.DOX*RP**3xSINSEC )*( 1.00*2.DOxDTAN( RPHI )x2 
$ *ESQPX*DCOS ( RPHI )3«2 )*1.018/KO**3 


EPRIME=EAST-500000.00 
Q = .00000100*EPRIME 
ажжа 
Q**3 
Q**4 
Qx*5 
Qx*6 
D6éNUM=Q6*DTAN( RPHI } 
D6z( D6NUM )/( 720. DO*RP**6*SINSEC )*( 61.00*90.DOx*DTAN(RPHI )*x2 
$ *645,D0*0TAN(RPHI )**4*107.DO0*ESQP*DCOS ( RPHI )x*2 
$ -162.DO0x*ESQP*DSIN( RPHI )**x2-45.00xESQP*DTAN(RPHI )*2 
$ XDSIN(RPHI )**2)*x1.D36/KO0**6 
E5bzZ(Q5*1.D0/DCOS(RPHI ) )7( 120. DO*RP**S5xSINSEC )*( 5.D0*28.D0 
$ ХОТАМ( КРНТ ) ее 2+24.ООЖОТАМ( КРНТ )**4+6 .DO*ESQP*DCOS( RPHI )*«2 
$ *8.DO*ESQP*DSIN(RPHI )**2)*1.D30/K03*5 


O 
£. 
uw n ou M 


DPHI-z(-R7*Q2*R8x*Q6-D6)/3600.D0 
DLAMz(R9*Q-R10*xQ3*E5)/3600.D0 
PHI=PHI+DPHI 
LON=CM+DLAM 
CALL OMS(PHI,LADEG,LAMIN,LASEG ) 
CALL OMS(LON,LODEG,LOMIN,LOSEG ) 
RETURN 
END 
C 
(3 3€3€ 9€ €3€ € 9€3€ € 3€ EOCOOCOE EE ODE XE 9€ зема IE FE IESE X жзеж жззезеезвезазезезежзезвезезезезезезезезезезезжезезезезе 
SUBROUTINE OMS(DEC,LDEG,MIN,SEC) 
QC 36(00)(009000000000)0000)00€ € X EIEEE IEE HEHE IE JEJE € 363€ X190 XC X€€ € HE JEJE JEJE € 90€ 99€ HE E JEJEJE FE 
С 
DOUBLE PRECISION DEC,;XNUM;,XMIN 
C 
XNUMSDABS(DEC ) 
LDEG=DINT ( XNUM ) 
XMIN=( XNUM-DFLOAT( LDEG ) )¥60.00 
MIN=DINT(XMIN } 
XNUMz( XMIN-DFLOAT ( MIN) )*60.D0 
SEC=SNGL ( XNUM ) 
IF (LOEG.GE.360) LDEG=LDEG-360 
RETURN 
END 
/ ж 
//GO.FT06F001 DD SYSOUT=* 
//60.ЕТ07Ғ001 00 05М-М55.50812.ҒА(СОМ.АМТЕММА.РО5,015Р-010 


09 


//60.ЕТ08Ғ001 DD DSNzMSS. 
//GO.FTSOFOO1 OD OSN=MSS. 


КА 
// 
/7 
/7 
/7 
// 
/7 


DSN-MSS. 
DSN=MSS. 
DSN=MSS. 
OSN=MSS. 
DSN=MSS. 
DSN=MSS. 


$0812. 
$0812. 
50812. 
50812. 
50812. 
50812. 
50812. 
50812. 


FALCON. 
FALCON. 
FALCON. 
FALCON. 
FALCON. 
FALCON. 
FALCON. 
FALCON. 


COURSE.ROUGH, DISP=OLD 

DATA(FALCON1),DISP=SHR 
DATA( FALCON2 ) ,DISP=SHR 
DATA( FALCON3 ),DISP=SHR 
DATA( FALCONG ),DISP=SHR 
DATA( FALCONS ) ,DISP=SHR 
DATA( FALCON6 ) ,DISP=SHR 
DATA( FALCON7 ),OISP=SHR 
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APPENDIX B 
PROGRAM COURSE SMOOTH. SOURCE LISTING 


/ / JOBSMOTH JOB (0812,9999),'EZEQUIEL',CLASS-C 
//ЖМАТМ ORG=NPGYM1.0812P ,LINES=( 99) ,CARDS=(99) 
//*FORMAT PR,DDNAME=GO.FTO6FO01, 

//*FORMS=SEP1 

/” ЕХЕС FORTVCLG 

ZZFORT.SYSIN DD ж 


AUTHOR: AUGUSTO EZEQUIEL 
DATE: JANUARY 05, 1986 


DESCRIPTION: 

THE PROGRAM MAKES THE RUNNING AVERAGE OF THE COURSE COMPUTED 
BY THE PROGRAM « COMP > IN ORDER TO TAKE OUT THE EFECTS OF THE 
PITCH AND ROLL EFFECTS. THE FIRST AND LAST THO POSITIONS ARE 
REJECTED AS NO RUNNING AVERAGE IS POSSIBLE FOR THEM. 

THE PROGRAM RUNS IN MVS 

I/O SPECIFICATIONS: SEE END OF THIS JOS. 


THIS PROGRAM IS RESTRICT AS IS TO THE SEAFLOOR BENCHMARK EXPERIMENT 


THE PROGRAM WILL RUN WITH ANY AMOUNT OF DATA ONLY LIMITED TO DISK 
SPACE 


ANY BLANK LINES WILL TERMINATE THE PROGRAM IN ERROR 


опп0попо0попопопо опопопопопопо 


REAL COURSE(5),SECS(5) 
INTEGER MONTH(5),YEAR(5) ,DAY(5),HOUR(5)},MIN(5),FLAG 
LOGICAL SAVE 


E 
© INITIALIZATION OF FLAGS 
С 
SAVE= .FALSE. 
FLAG=1 
С 
C READS THE FIRST FIVE COURSES 
С 


DO 100 J-1,5 
READ(51,END=1100) MONTH(J),DAY(J),YEAR(J),HOUR( J), 
$ MIN(J),SECS(J),COURSE( J) 
150 РОКМАТ(2(1Хх,12),1Х,14,„1<,12,1Х,12,1Х,Е4.1,1Х,ЕО07.2) 
С WRITE(6,150) MONTH(J),DAY(J),YEAR(J),HOUR( J), 
С $ MIN(J),SECS(J}),COURSE(J) 
100 CONTINUE 


E AVERAGE THE COURSE 

E 105 AVGs(COURSE ( 1)*«COURSE ( 2 ) COURSE ( 3 ) COURSE ( 4 ) eCOURSE( B ) )/5.0 
E SAVES THE AVERAGE IN THE MIDDLE POSITION 

ў COURSE ( 3 )=AVG 

: DOES NO SAVE THE FIRST TWO POSITONS OF THE FILE AS NO RUNNING 
C AVERAGE EXISTS FOR THEM 

: IF(SAVE) GOTO 50 

С 
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FLAG=FLAG+1 
IF(FLAG.EQ.3) SAVE=z.TRUE. 
GOTO 70 


WRITES THE DATA WITH THE AVERAGED COURSE INTO THE FILE 


ооо 


50 WRITE(7) MONTH(1),DAY(1),YEAR(1),HOUR(1) ,»MIN(1),SECS(1),COURSE(1) 


WRITE(6,150) MONTH(1),DAY(1),YEAR(1),HOUR(1), 
$ MIN(1),SECS(1),COURSE(1) 


MOVES FORWARD THE DATA 


оо О 


70 CONTINUE 
DO 75 J=2,5 
MONTH(J-1)=MONTH( J) 
DAY( J-1J=DAY( J) 
YEAR( J-1 JSYEAR( J) 
HOUR( J-1 J=HOUR( J ) 
МІМ( J-1 )=MIN( J) 
SECS( J-1)=SECS( J) 
COURSE ( J-1 J)=COURSE( J) 

75 CONTINUE 


READS ONE MORE COURSE 


оо О 


READ(51,END=1000) MONTH(5),DAY(5),YEAR(5),HOUR(5), 
$ MIN(5),SECS(5),COURSE(5) 


GO BACK TO COMPUT A NEW AVERAGE 
GO TO 105 


SO SAVES THE REMAINING TWO AVERAGES 
000 CONTINUE 
DO 1001 J-1,? 
WRITE( 7,150} MONTH(J),DAY(J),YEAR(J),HOUR( J), 
$ MIN(J),SECS(J) ,COURSE( J) 
WRITE(6,150) MONTH(J),DAY(J),YEAR(J),HOUR( J), 
$ MIN(J),SECS(J),COURSE( J) 
1001 CONTINUE 
STOP 


C 
C END OF FILE FOUND DURING THE INPUT OF A FITH COURSE 
E 
1 


OF THE FIRST FIVE COURSES 
100 CONTINUE 
НКТТЕ (6,140) 
140 FORMAT (' ERROR. END OF FILE DURING THE FIRST SET OF COURSES') 
STOP 
END 


C 
C NO OUTPUT IS MADE IF ΑΝ END OF FILE WAS FOUND DURING THE READING 
C 
1 


/* 

//GO.FTO6FOO1 DD SYSOUT = 

//GO.FTO7FOO1 DD DSN=MSS.S0812.FALCON.COURSE .SMOOTH ,DISP=SHR 
//GO.FT51F001 DD DSN=MSS.S0812.FALCON.COURSE . ROUGH ,DISP=SHR 
ди 


APPENDIX C 
PROGRAM TRANSDUC FALCON. SOURCE LISTING 


// JOBTRANS JOB (0812,9999), ‘EZEQUIEL’ ,CLASS=C 
//ЖМАІМ ΟΒΟΞΝΡΟΝΗΙ .0812Ρ.Ι ΙΝΕΞΞ( 99),ΟΑΒΟΞΞ( 99) 
{/ΧΕΟΒΜΑΤ ΡΕ ,ΟΡΝΑΜΕΞΟΟ.ΕΤΟέΕΕΟΟΙ, 

//*FORMS=SEP1 


// ЕХЕС FORTVCLG 
ZZFORT.SYSIN DD * 
C 


AUTHOR: AUGUSTO EZEQUIEL 
DATE: MARCH 04, 1987 


DESCRIPTION: 


THIS PROGRAMS TAKES EACH POSITION OF THE ANTENNA 

THE COMPUTED COURSE ( FILTERED OR NOT ) ,THE PITCH 

AND ROLL DATA, AND COMPUTES THE POSITIONS OF THE 
TRANSDUCER, APLYING THE SEVEN PARAMETER TRANSFORMATION 

TO THE OFFSETS OF THE TRANSDUCER IN RELATION TO THE ANTENNA 
CONSIDERED THE CENTER OF COORDINATE SYSTEM AND THEN ADDING 
THE CORRECTIONS TO THE COORDINATES. 


THE PROGRAM RUNS IN MVS 
I/O SPECIFICATIONS: SEE ENO OF THIS JOB. 


THIS PROGRAM IS RESTRICT AS IS TO THE SEAFLOOR BENCHMARK 
EXPERIMENT 


THE PROGRAM WILL RUN WITH ANY AMOUNT OF DATA ONLY LIMITED TO DISK 
SPACE 


ANY BLANK LINES WILL TERMINATE THE PROGRAM IN ERROR 


оппоппоппопоооопопопопопоопоопопопо 


DOUBLE PRECISION XPOS,YPOS,TIME1,TIME2,TIME3,SECD,SECH,RATE 
DOUBLE PRECISION TIME¢, TIMES 

REAL COURSE »PITCH( 2),ROLL(2),SEC1(2),SEC2,PIT ,ROL »SEC3(2) 
REAL 0X,DY,0Z,OFFX,OFFY ,OFFZ,HEAD(2),LASEG,LOSEG 

INTEGER MONTH1(2),DAY1(2) , YEAR1( 2) ; HOUR1( 2) ;MIN1(2) 
INTEGER MONTH2 ;DAY2 »YEAR2 »HOURZ »мтмг 

INTEGER MONTH3( 2) ,DAY3(2),YEAR3S(2),HOUR3(2),MIN3(2) 
INTEGER LADEG,LAMIN,LODEG,LOMIN 


Е 
С INITIALIZATION OF CONSTANTS 
C 
SECH=3600.00 
SECD=24.D0*SECH 
C 
OFFX2-5.365 
ΟΕΕΥΞ- 9.455 
OFFZ=-16.152 
C 
С READS THO SETS OF PITCH AND ROLL DATA 
C 
READ(51,%*,END=z100) HONTH1(1),)DAY1(1),YEAR1(1),HOUR1(1),MIN1(1), 
$ SEC1(1),PITCH(1),ROLL(1) 
TIME1=DFLOAT( DAY1(1) )*SECD+DOFLOAT(HOUR1(1 ) )*SECH 
$ * DFLOAT(MIN1(1))*60.D0*DBLE(SEC1(1)) 
C 


READ(51,*,ENDz100) MONTH1(2),DAY1(2), YEAR1( 2) ;HOUR1( 2) ; MIN1(2); 
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5 SEC1(2),PITCH(2),ROLL(2) 
TIME 2=DFLOAT( DAY1( 2} }*SECD+OFLOAT(HOURL( 2) )*SECH 
$ * DFLOAT(MIN1(2))*60.00*DBLE(SECI(2)) 
READS THO SETS OF COURSE 


ВЕАО( 52 ,ЕМО=100 ) MONTH3(1),0AY3(1),VYEAR3(1),HOUR3(1),MIN3S(1), 


5 SEC3(1),HEAD(1) 
TIMES=DFLOAT( DAY3( 1) )*SECD+DFLOAT( HOUR3( 1) )*¥SECH 

$ + DFLOAT(MIN3(1))*60.00+DBLE(SEC3(1)) 
READ(52,END=100) MONTH3(2),DAY3(2),YEAR3(2),HOURS(2),MIN3(2), 

$ SEC3(2),HEAD( 2) 
TIMES=DFLOAT( DAY3( 2) J¥SECD+DFLOAT( HOURS( 2) )*¥SECH 

$ + DFLOAT(MIN3( 2) )¥60.D0+DBLE( SEC3(2)) 


READS ONE POSITION OF THE SHIP 


20 READ(50,END=100) MONTHZ ,DAY2 ,YEARZ,HOURZ »MIN2 »SEC2, 
$ LADEG,LAMIN ,LASEG,LODEG,LOMIN,LOSEG 


TIMES-DFLOAT(DAY2 )X*SECD*DFLOAT ( HOUR2 )* SECH 
» + DFLOAT(MIN2)%60.D0+DBLE(SEC2 ) 


POSITION IN ΤΙΜΕ ΒΕΤΗΕΕΝ THE TWO SETS OF PITCH AND ROLL DATA 
AND THE THO SETS OF COURSE DATA 
25 CONTINUE 


POSITION TO EARLY FOR PITCH AND ROLL DATA 


IF(TIMES.GE.TIME1) GO TO 40 
ARITE(6,30) MONTH2;,DAY2,YEAR2;,HOUR2,MIN2,SEC2 
30 FORMAT(' PSN 3 ',I2,1X,I12,1X,I4,1X,I2,1X,12,1X,F4.1, 
$ ' REJECTED. NO PITCH AND ROLL DATA') 
GOTO 20 


POSITION TO EARLY FOR PITCH AND ROLL DATA 


40 CONTINUE 
IF(TIMES.GE.TIME4) GO TO 60 
WRITE(6,50) MONTH2,DAY2,YEAR2,HOUR2;MIN2,SEC2 
50 FORMAT(' PSN 39 ',I2,1X,I2,1X,I4,1X,I12,1X,I12,1X,F4.1, 
$ ' REJECTED. NO COURSE DATA') 
GOTO 20 


POSITION LATER THAN THE SECOND SET OF PITCH AND ROLL DATA 


60 CONTINUE 
IF(TIMES.LE.TIME2) GOTO 70 


MOVES THE PITCH AND ROLL DATA TO THE FIRST SET 


MONTH1(1)zMONTH1(2) 
DAY1(1)=DAY1( 2) 
ҮЕАК1(1 ) =ҮЕАВ1(2) 
HOUR1( 1) < НОИ 1 (2) 
МІМ1(1)-МІМ1(2) 
SEC1(2)=SEC1( 2) 
PITCH(1)=PITCH(2) 
ROLL(1)=ROLL( 2) 
ΤΙΜΕΙΞΤΙΜΕΖ 


READS А SECOND SET OF PITCH AND ROLL DATA 


READ(51;,*,ENDz2100) MONTHI1(2),DAY1(2),YEAR1( 2)  HOUR1( 2) ; MIN1( 2), 
$ SEC1(2),PITCH(2);,ROLL(2) 
TIME2-zDFLOAT(DAY1( 2) )J*SECD*DFLOAT ( HOUR1( 2 ) JX SECH 
$ * DFLOAT(MINI1(2))*60.00*DBLE(SEC1( 2)) 
GOTO 25 
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POSITION LATER THAN THE SECOND SET OF COURSE DATA 


70 CONTINUE 
IF(TIME3.LE.TIMES) GOTO 80 


MOVES THE COURSE DATA TO THE FIRST SET 


MONTH3 (1 )=MONTH3( 2) 
DAY3(1)=DAY3(2) 
YEAR3(1 J=YEAR3(2) 
HOUR3( 1 )=HOUR3( 2) 
MIN3(1)zMIN3(2) 
SEC3(2 J=SEC3( 2) 
HEAD( 1 )=HEAD( 2) 
ТІМЕФ-ТІМЕБ 


READS A SECOND SET OF PITCH AND ROLL DATA 
READ(52,END=100) MONTH3(2),DAY3(2),YEAR3(2),HOUR3(2),MIN3(2), 
$ SEC3(2),HEAD( 2) 
TIME 5=DFLOAT ( DAY3( 2 ) }*SECD+DFLOAT( HOUR3( 2 ))%SECH 
$ + DFLOAT(MIN3( 2 ) )*60.00+DBLE(SEC3(2)) 
GOTO 25 
POSITION WITHIN THE TWO SETS IN TIME 
80 RATE=( TIME3-TIME1 )/( TIME2-TIME1 ) 
PIT=PITCH( 1 )+(PITCH( 2 )-PITCH( 1) )¥RATE 
ROL=ROLL(1 )+( ROLL( 2 )-ROLL( 1) )¥RATE 
RATE=( TIME3-TIME4¢ )/( TIMES-TIMEG ) 
COURSE =HEAD( 1 }+( HEAD( 2 )-HEAD(1) )*RATE 
COMPUTS THE CORRECTIONS AND THE POSITIONS OF THE TRANSDUCER 
CALL PTHLOFFX,OFFY ,OFFZ;,COURSE ;PIT ;ROL ; DX, DY ,DZ) 
COMPUTES THE UTM COORDINATES 
CALL GPUTM(LLADEG,LAMIN;LASEG,LODEG,LOMIN,LOSEG,YPOS,XPOS) 


XPOS=XP0S+DX 
Y¥POS=YPOS+DY 


CONVERTS BACK TO GP 
CALL UTMGP( YPOS,XPOS ,LADEG,LAMIN,LASEG »,LODEG, LOMIN »LOSEG ) 
OUTPUTS THE RESULT 


WRITE(7) MONTH2,DAY2,YEAR2HOUR2 , MIN?2 ; SEC2 
5 LADEG,LAMIN,LASEG,LODEG ,LOMIN,LOSEG 


WRITE(6,90) MONTH2,DAY2,YEAR?2 ;HOUR2 ,MIN2 , SEC2, 
< LADEG,LAMIN,LASEG,LODEG,LOMIN,LOSEG 
90 FORMAT(2(1X,I2),1X,I6,1X, I2, 135 I2, 1X, F6. 1, 2C1X 13, 1X , 12, 135 F7.6)) 


READS ONE MORE POSITION 
GOTO 20 


END OF PROGRAM 
100 CONTINUE 
STOP 
END 
SUBROUTINE PTH(XOFF,YOFF,ZOFF,COURSE,PITCH,ROLL,DX,DY,DZ) 


THIS SUBROUTINE COMPUTES THE CORRECTIONS TO THE 
COORDINATES DUE TO THE OFSET OF ANTENA IN RELATION 
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ΤΟ THE TRANSDUCER 


PI=ARCOS(-1.0) 
ROLL-ROLLxPI/180. 
ΡΙΤΟΗΞΡΙΤΟΗΧΡΙΖ180. 
COURSE =COURSE*PI/180. 


CHANGE SIGNAL OF ORIGINAL DATA TO MATCH THE RIGHT HAND 
SYSTEM CONVECTION 


ROLL=-ROLL 


COSPHI= COS(ROLL ) 
SINPHI- SINCROLL ) 
COSK= COS(COURSE ) 
SINK- SIN( COURSE ) 
COSW= COS( PITCH ) 
SINW= SIN(PITCH ) 


DX = XOFF * COSPHI ж COSK + 


$ YOFF ж ( COSW * SINK + SINW * SINPHI * COSK ) + 
$ ZOFF * ( SINW * SINK - COSW * SINPHI ж COSK 
DY = XOFF * COSPHI * SINK * (-1.) + 
$ YOFF * ( COSW * COSK - SINW * SINPHI * SINK ) + 
$ COFF * ( SINW * COSK + COSW * SINPHI * SINK ) 
DZ = XOFF * SINPHI - 
$ YOFF * ( SINW * COSPHI ) + 
$ СОҒҒ * ( COSW * COSPHI ) 
RETURN 
END 


SUBROUTINE UTMGP(NORTH,EAST ,LADEG,LAMIN, LASEG ,LODEG,LOMIN,LOSEG ) 
DOUBLE PRECISION A,R>N,B, BP,CP,DP,EP,S,R1,ESQ,ESQP ,RM,>RP ,KO 
DOUBLE PRECISION R7,R8,R9,E5,0,Q2,Q03,Q6,Q5,Q6,D6 ,AP,SINSEC 
DOUBLE PRECISION R10,DLAM,NORTH;EAST;RPHI,DPHI ,PHIMIN 

DOUBLE PRECISION EPRIME,DELTA,D6NUM,CM>PHI,LON,PI 

REAL LASEG,LOSEG 


THIS SUBROUTINE COMPUTS THE UTM COORDINATES OF GP IN WGS 72 
IN ZONE 10 CENTRAL MERIDIAN 123 00 00 W 


CM--123.0D0 


Аз6378135.000 
К-298.2600 


KO=0.9996D0 


PI-DARCOS(-1.D0) 
В - Ax(R-1.D0)/R 


N * (A-B)/CA*B) 


АР = Ax((1.D00-N)*5.D0/4 . DOx( NXPF2 -N3X3 ) 81. D0/66 . DOx ( N39€6 -№өє5 ) ) 

ВР = 3.D0/2.D0*Ax( (N-Nxe2 ) *7.D0/8.DO*( N3ex 3 -N9e9x$ ) 55. D0/66 . DOXN9€5 ) 
CP = 15.00/16.00жАж( Мжзе2 -Мжз3+3.00/4.0036( Мжзеф-Межб ) ) 

DP = 35.00/48.00жАж( Мю3 -Мжж4+11.00/16.00жМжзБ ) 

ЕР = 315.00/512.00жАж( Междф-Мезеб ) 


FIRST APROXIMATION OF PHI 
PHI-NORTH/30.8D0/3600.D0 
COMPUT TRUE MERIDIONAL DIST AND APROXIMATE PHI 


DO 100 I-1,8 
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PHIMINZPHI*60.D0x2.90888208666D-4 

RPHI-PHI/180.DO0xPI 

5 < APXxPHIMIN-BP*DSIN( 2.DOx*RPHI )«CP*DSIN( 4. DOxRPHI) 
-DPx*DSIN(6.DOx*RPHI)«EPXDSIN(8.DOXRPHI ) 

R1 - KOxS 

DELTA=NORTH-R1 

PHT =(DELTA/30.8/3600. )+PHI 


100 CONTINUE 


$ 


$ 
$ 
$ 


$ 


$ 
$ 
$ 


$ 
$ 


RPHI-PHI/180.DO0xPI 


SINSEC 
SINSEC 


(1.D0/3600D0)/180.D0xPI 
DSIN(SINSEC ) 


ESQ с (Ажж2-Вжж2)/Ажж2 
ESQP = ESQ/(1.D0-ESQ) 
RM * AX*(1.DO-ESQ)/( DSQRT(1.DO-ESQxDSIN( RPHI )*x2 ) )x3 


КР = RM*(1.DO0«*ESQPXxDCOS( RPHI )x2) 
R7-2( DTANCRPHI2/(2.DOXRPX*2xSINSEC ) )*( 1. DO«ESQPXDCOS( RPHI )x2) 
x1.D12/KOxx2 
R8z( DTAN(RPHI )/( 26 , DOXRP2&xSINSEC ) )x( 5. DO«2. DOXDTAN(RPHI 1192 
+6. DOXESQP*DCOS( RPHI )x*2-6. DOXESQP*DSIN( RPHI )**2-3 . DOXESQPx*2 
*DCOS( RPHI )**4-9 , DOXESQPx*2*DCOS( RPHI )**2*DSIN( RPHI )**2) 
x1.D264/KO0»*4 
R9z1.D0/DCOS( RPHI )/( RP*«SINSEC )x1.D6/KO 
R10-1.D0/DCOS( RPHI )/( 6. DOXRPX*xZxSINSEC )*(1.D0«2.DOxDTAN(RPHI )же2 
*ESQPx*DCOS( RPHI )*«2)x1.D18/KOx*x3 


EPRIME-EAST-500000.D0 
Q - .000001D0*EPRIME 


Q2 с 0ж%Ұ2 
95 = 4 
Q4 = QxxG 
95 = д5 
Q6 - Q**6 


D6NUMzQ6 *DT AN(RPHI ) 

D6z( D6éNUM )/( 720. DOxRPxx*6xSINSEC )*( 61. D0*90.DOxDTAN( RPHI )*x2 
*$5.DOxDTAN(RPHI )4*107.DOxESQP xDCOS( RPHI )*x2 
-162.D0xESQPXxDSIN( RPHI )3&«2-645. DO*ESQP xDT AN(RPHI )sx2 
*XDSIN(GRPHI)**2)x*x1.D36/KOxx6 

E5-(Q5*1.D0/DCOS( RPHI ))/( 120. DO*RPx*5xSINSEC )*(5.D0*28.DO 
XDTAN(CRPHI )*2* 26. DO*DT ANCRPHI )4 «6 . DOXESQPxDCOS ( RPHI )*«2 
*8. DOXESQPX*DSIN( RPHI )**2)x1.D30/KO0*«5 


DPHI-(-R7*Q2*R8*Q4-D6)/3600.D0 
DLAM=( R9¥Q-R10*Q3+E5 )/3600.D0 
PHI=PHI+DPHI 

LON=CM+DLAM 

CALL DMS(PHI,LADEG,LAMIN,>LASEG ) 
CALL DMS(LON,LODEG,LOMIN,LOSEG ) 
RETURN 

END 


SUBROUTINE GPUTM( LADEG,LAMIN,LASEG,LODEG,LOMIN,LOSEG»NORTH,EAST ) 


DOUBLE PRECISION A,R,N,AP,BP,CP,DP,EP,S,R1,ESQ,ESQP ,RM,RP ,KO 
DOUBLE PRECISION R2;R3;R6,R5,P,P2,P2,P4,P5,P6,A6,85,SINSEC 
DOUBLE PRECISION PHI,DLAM;,NORTH,;EAST,;B;PHIMIN,;PI,;LON 

REAL LASEG,LOSEG 


THIS SUBROUTINE COMPUTS THE UTM COORDINATES OF GP IN WGS 72 


IN 


ZONE 10 CENTRAL MERIDIAN 123 00 00 W 


CM=-123.0D0 

PHI=DFLOAT( LADEG )+DFLOAT(LAMIN)/60.DO+DBLE( LASEG )73600.D0 
LON=DFLOAT( LODEG )+DFLOAT( LOMIN )/60.D0+D8LE( LOSEG )/3600.D0 
LON=-LON 

DLAM=( LON-CM )*3600.0D0 


TI 


$ 


$ 


$ 


$ 
$ 


$ 
$ 


A=6378135.000 
R=298.2600 


KO=0. 999600 
B = AX(R-1.00)/R 


М = (А-В )/(А+8 ) 


АР * A*((1.00-N)*5.00/64.00*( NXPE2-N93 ) 81.007664 , QO* ( N39EQ -N3XES ) ) 

BP = 32.00/2.D0*A*( (N-Nx*2 )*7.00/8. 00*( Nxx3-N»*à ) 55. 00/664 . DOXN»&5 ) 
СР = 15.00/16.00*A*( Nx*x2-N**3«3.00/4 . Q0*( NX*&4 -N3xe5 ) ) 

OP = 25.00/48.00*A*( NX*3-N**4«11.D0/16.D0xNx*5) 

EP = 315.D0/512.00*A*( N**xG -N**5 ) 


PHIMIN - PHI*60.00*2.90888208666D-4 

PI-OARCOS( -1.00) 

РНІ=РНІ/180.00#%РІ 

S  АРЖРНТМТМ-ВРЖОЗТМ( 2. ООЖРНТ ) +СР05ТМ( 4 .ООЖРНТ ) 
- DPXDSIN(6.00*PHI )+EPXOSIN(8.DOXPHT ) 

R1 - KOx*S 


SINSEC = (1.00/3600.00)/180.00*PI 
SINSEC=OSIN(SINSEC ) 


ESQ - (Axx2-B»xx2)/Ax«2 
ESQP - ESQ/(1.00-ESQ) 
RM 7 A*(1.00-ESQ)/COSQRT(1.DO-ESQXOSIN(PHI )*x2 ) )e3 


RP s RMx*(1.00«*ESQPXDCOS( PHI 9x2) 
R2 = RPXOSIN(PHI)XDCOS(PHI)XSINSEC**2/2.D0xKOX1.008 
R3 = SINSEC*x4xRPXDSIN( PHI )«0COS ( PHI 9€2/26 . 00x( 5. 00-DTAN( PHI «2 
* 9.DO0XESQPX*DCOS( PHI 3x26. DOXESQPXESQPXOCOS ( PHI 36x6 )xKOX1. 016 
Ré - RPXDCCS(PHI)XSINSECXKOx1.04 
R5 s SINSECX*3xRPX0COS( PHI 12/6 .00x( 1. 00-0TAN(PHI x2 
* ESQGPX*OCOS( PHI )xx2 )*KOx1.012 
Р = .000100*0LAM 
рг = ржж2 
P3 = P**3 
рд с ржжд 
P5 = P**5 
P6 = P%*6 
AG = P6XSINSECX*6*RPXOSIN( PHI )¥0COS( PHI )**5/720.00 


* (61.00-58.00%0TAN( PHI )x=2+OTAN( PHI )w**&% 
+ 270.00*ESQP*OCOS ( PHI ) е 2-350. ООЖЕЗаРЖОЗТМ(РНТ) 2) КО 1.024 
B5 z PS5xSINSECX*5xRPx*0COS(PHI)**x5/120.D0*(5.00-18.DO0Xx0TAN( PHI )x*2 
+ ОТАМ(РНТ ) 4 + 14. ООЖЕЗаРЖОСОЗ(РНТ ) 2 
- 58.00X*ESQPx*OSIN(PHI )**2 )*KO*1.020 
МОКТН<К1 +К2хР2 + КЪжР4 + д6 
ЕАДЗТ-( Н4ажР +КБжР3 +В5 )+500000.00 
КЕТОКМ 
ΕΝΟ 


SUBROUTINE OMS(OEC,LOEG,MIN,SEC ) 
DOUBLE PRECISION DEC ,XNUM,XMIN 


XNUM=DABS( OEC ) 

LOEG=OINT( XNUM ) 

XMIN-( XNUM-OFLOAT ( LOEG ) )*60.D0 
MIN-OINT ( XMIN) 

ХМОМ = ( XMIN-DFLOAT ( MIN) )*60.00 

SEC=SNGL( XNUM ) 

IF (LOEG.GE.360) LOEG=LOEG-360 
RETURN 

ENO 
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PROGRAM PLOT. SOURCE LISTING 
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PROGRAM PLOT 
RUNS IN FORTRAN VS 
AUTHOR: AUGUSTO EZEQUIEL 
DATE : 26 JANUARY 1987 
THIS PROGRAM MAKES A PLOT OF THE POSITIONS OF THE SHIP 
DOUBLE PRECISION XPOS,YPOS,XLEFT,YLEFT,TIMEI,TIME2,SECD,SECH 
REAL SECS,XPLT,YPLT,BLXH,BLYH,;SCALE VALUE ; LASEG, LOSEG 
INTEGER MONTH,YEAR,;,DAY ;HOUR, MIN,IPEN,LADEG,; LAMIN,;LODEG,LOMIN 
PLOTTER INITIALIZATION 
CALL PLOTS (0,0,0) 
CALL PLOT (2:52. 5-3] 
IPEN=3 
DIMENSIONS OF SHEET, LEFT CORNER AND SCALES 
BLXH=36. 
BLYH=49. 
XLEFT-566800.D0 
YLEFT=403 7000.00 
SCALE=1./10000. 
INITIALIZATION OF CONSTANTS 
TIME1=0.00 
SECH=3600.00 
SECD=24¢.D0*SECH 
GRID 
CALL PLOT (0.1,0.1,3) 
CALL PLOT (0.1,BLYH-0.1,2) 
CALL PLOT (BLXH-0.1,BLYH-0.1,2) 
CALL PLOT (BLXH-0.1,0.1,2) 
CALL PLOT (0.1,0.1,2) 
TITLE 
CALL TITLE 
SCALE 
CALL METER(SCALE ) 


READ THE DATA 


10 CONTINUE 


REAOD(51,END=50 JMONTH,DAY »YEAR,HOUR, 
$ MIN,SECS,LADEG,LAMIN,LASEG,LODEG ,LOMIN,LOSEG 


SO 
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COMPUTES THE UTM COORDINATES 


CALL GPUTM(LADEG, LAMIN,LASEG,LODEG;LOMIN,LOSEG  YPOS ,XPOS ) 


COMPUTES THE TIME IN SECS 


TIME 2=SECD*FLOAT( DAY )* SECHXFLOAT ( HOUR ) $60. DOXFLOAT(MIN)4SECS 


COMPUTES THE PLOTTER COORDINATES 


XPLT-(XPOS-XLEFTO*SCALE*100. 
YPLT-z(YPOS-YLEFTO*SCALE*100. 


TESTE IF INSIDE AREA 


IF(XPLT.LT.0.) GOTO 30 
IF(YPLT.LT.O.) GOTO 30 
IF(XPLT.GT.BLXH) GOTO 30 
IF(YPLT.GT.BLYH) GOTO 30 


TESTS IF POSITIONS ARE AWAY MORE THEN 10 SECS IN TIME 


IF(TIME2-TIME1.GT.10.D0) IPEN-3 
TIME1-TIME2 


PLOTS THE POSITION 


CALL PLOT(XPLT,YPLT,IPEN) 
IPEN=2 


PLOTS THE TIME EVERY O5 MINUTES 


IF(MOD(MIN,05).GT.O.OR.IFIX(SECS).GT.0) GO TO 10 
VALUE -FLOAT( FOUR )*100 *FLOAT(MIN) 
CALL NUMBER(XPLT*0.15,YPLT*0.15,0.25, VALUE 0. ; - 1) 


CALL SYMBOL(XPLT,YPLT,0.15,3,0. 5-1) 


GOTO 10 


PEN UP WHILE THE POSITIONS ARE OUT OF THE SHEET 


30 IPEN=3 
GOTO 10 


END OF PLOT 


50 CALL PLOT(O. ,0. , +999) 
STOP 
END 


SUBROUTINE METER( SCALE ) 


X0224.0 
YO=2.0 
CALL PLOT(X0,Y0>,-3) 


DO 10 «51,12 

XP» FLOAT(J)-1.0 

CALL PLOT(XP,0.0,3) 

CALL PLOT(XP,0.25,2) 
10 CONTINUE 


DO 20 .51,9 

XPzFLOAT(J)XO.1 

CALL PLOT(XP,0.0,3) 

CALL PLOT(XP,0.2,2) 
20 CONTINUE 


CALL PLOT(0.,0.0,3) 
CALL PLOT(11.,0.52) 


81 


G OoOo 


O 


CALL SYMBOL(00.0,1.0,0.25,29H 


CALL SYMBOL(00.0,1.0,0.25,29H 


CALL SYMBOL(00.0,0.5,0.25,29H 


CALL PLOLCO- 30-235) 
CALL PLOT(11.,.2,2) 


VALUE=1./SCALE/100. 

CALL NUMBER(-0.25,0.27,0.25,VALUE,0.0,-1) 
CALL NUMBER(O.95,0.27,0.25,0.0,0.0,-1) 
VALUE z25.0/SCALE/100. 

CALL NUMBER(5.775,0.27,0.25, VALUE, 0.0, -1) 
VALUE=10.0/SCALE/100. 

CALL NUMBER(10.7,0.27,0.25,VALUE ,0.0,-1) 
CALL SYMBOL(5.5,-0.50,.25, 6HMETERS,0.0,6) 


CALL SYMBOL(5.5,24.0,2.05,62,0.0,-1) 
CALL SYMBOL(5.7,24.5,0.4,85,0.0,-1) 
ΧΟΞ-ΧΟ 

YOz-YO 

CALL PLOT(XO,YO;-3) 


RETURN 
END 


SUBROUTINE TITLE 
XO=26.5 


YO=43.05 
CALL PLOT(XO,YO,;-3) 


CALL SYMBOL(00.0,3.5,0.25,29H NAVAL POSTGRADUATE SCHOOL 


90.929) 


CALL SYMBOL(00.0,3.0,0.25,29HSEAFLOOR BENCHMARK EXPERIMENT >0. 29) 


CALL SYMBOL(0OO.0,2.5,0.25,29H PHASE II 


CALL SYMBOL(00.0,2.0,0.25,29H 


CALL SYMBOL(OO.0,1.5,0.25,29H 


CALL SYMBOL(OO.0,1.0,0.25,29H 


CALL SYMBOL(0O.0,0.5,0.25, 29H ANTENNA POSITIONS 


CALL SYMBOL(0OO.0,0.0,0.25, 29H LCDR AUGUSTO EZEQUIEL 


CALL PLOT(-0.2,-0.2)3) 
CALL PLOT(-0.2,3.95,2) 
CALL PLOT(7.5,3.95,2) 
CALL PLOT(7.5,-0.2,2) 
CALL PLOT( -0.2,-0.252) 


XO2-XO 
ΥΟΞ-ΥΟ 
CALL PLOT(X0,Y0,-3) 


RETURN 
END 


R/V POINT SUR 


16/17 AUGUST 1986 


GESAR SOLUTIONS 


Т14100 


MINI RANGER FALCON 


TRANSDUCER POSITIONS 


;0. 


0. 


0. 


»0.,29) 
)0.»29) 
)0.,29) 
»29) 

»0.»29) 
29) 

»0. 529) 


)29) 


»0.,29) 


SUBROUTINE GPUTM( LADEG,LAMIN,LASEG,LODEG,LOMIN,LOSEG,NORTH >EAST ) 


DOUBLE PRECISION A,R,;N,AP,;BP,CP,DP,EP,S,Rl1,ESQ,ESQP,RM,RP,KO 
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DOUBLE PRECISION R2,R3,R6,R5,P,P2,P2,P4,P5,P6,A6,B5,SINSEC 
DOUBLE PRECISION PHI,0LAM,NORTH,EAST,B,PHIMIN,PI,LON 
REAL LASEG,LOSEG 


THIS SUBROUTINE COMPUTS THE UTM COORDINATES OF GP IN WGS 72 
IN ZONE 10 CENTRAL MERIDIAN 123 00 00 W 


CM=-123 .000 

PHI =DFLOAT( LADEG )+DFLOAT( LAMIN )/60.00+DBLE( LASEG )/3600.D0 
LON=D FLOAT ( LODEG ) *DFLOAT( LOMIN)/60.00«DBLE( LOSEG)/3600.D0 
LON=-LON 

DLAM=( LON-CM )*3600.0D0 


A=6378135.000 
R=298.26D0 


KO=0.9996D0 
B = A¥(R-1.D0)/R 


N = (A-B)/(A+B) 


AP = 4А5111.00-М)+5.00/4.00ж( Мее2-Мжжз )+81.00/64.003( Межд-Мее5)) 

ВР < 5.00/2.00ждж( (М- Мег )+7.00/8.00( Мее3 -Мжд )+55.00/64.ООзеНне» 5) 
CP < 15.00/16.00ждж( Мег -Ме3+5.00/4.003( Мжзед-КвееБ )) 

DP - 35.00/48.00*A*( Nxx3-Nx*x4«11.D0/16.D03N3e5 ) 

EP = 315.D0/512.D0*A*( Nx9*& -Nxx5 ) 


РНІМІМ = PHI*60.D0*2.908882086660-4 

PI=DARCOS(-1.00) 

PHI-PHI/180.00xPI 

S = AP¥PHIMIN-BPXDSIN( 2.D00*PHI )+CPxDSIN(4.D0PHI ) 
- DP%DSIN(6.D0%xPHI)+EP%DSIN(8.D0%wPHI) 

R1 = KO*S 


SINSEC = (1.00/3600.00)/180.D0xPI 
SINSEC=DSIN(SINSEC ) 


ESQ - (Asx*2-B»*2)/A*x2 
ESQP = ESQ/(1.D0-ESQ) 


RM = Ax(1.D0-ESQ)/(D0SQRT(1.D0-ESQ*DSIN( PHI ) 32 ) 33 

RP = RM¥(1.00+ESQP*DCOS( PHI )x*2) 

R2 = RP*XDSIN( PHI )*DCOS( PHI J¥SINSEC¥*2/2.D0*KOx1.0D8 

RS = SINSECH#G*RPXDSIN( PHI )*DCOS( PHI )**3/24.D0%(5.00-DTAN( PHI ) 2 


$ + 9.DO*xESQP*DCOS( PHI )**2+4.D0xESQPxESQPxDCOS( PHI )**4 )*KOx1.016 


R@ = RPX¥DCOS( PHI )*SINSEC*KO*1.DG 
R5 = SINSEC#*3*RP*XDCOS( PHI )**3/6.00%(1.00-DTAN( PHI )*2 
+ ESQPxDCOS( PHI )*x*2 )*KO*1.012 


P = .0001D0*DLAM 


P2 s P*«2 

PZ = Pez 

PG = PG 

P5 = peed 

P6 = P%%*6 

A6 3 P6X*SINSEC»**6xRP*DSIN(PHI )*x0COS( PHI ):5/720.D0 


x (61.00-58.D03x0TAN(CPHI )3€€2 «DT AN(PHI ) жә 
* 270.DO0X*XESQPX*XDCOS( PHI )*2-330.DO0XESQPxDSIN( PHI ) #2) КО» 1.024 
B5 * P5XSINSEC**5xRPxDCOS( PHI )*x*5/120.D0*(5.D0-18.00*xDTAN( PHI )**2 
+ DTANCPHI )*&«16 .00XxESQPX*O0COS ( PHI )*x2 
- 58.DOxESQPX*DSIN( PHI )**2)*xK0*1.020 
NORTHzRI*R2*P2*R3x*P&4*A6 
EASTs(RGxP«R5xP3«4B5)«4500000.DO0 
RETURN 
END 


//GO.PLOTPARM DO ж 


33 


&PLOT XMIN=0. »XMAX=38. ,YMIN=0. ,YMAX=51.,SCALE=1. ,UNITS=2.540 &END 
//60.ЕТ06Е001 DD SYSOUT=%w 


//GO.FT51F001 DD DSN=MSS.S0812.GPS.POS.ANTENNA ,DISP=SHR 
7 
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APPENDIX E 
PROGRAM CVFICA. SOURCE LISTING 


//JOBCVFI JOB (0812,9999), 'EZEQUIEL' ,CLASS=C 

/ /*XMAIN ORG=NPGVM1.0812P ,LINES=( 99) ,CARDS=( 99 ) 
//*FORMAT PR,DDNAME-GO.FTO6FO001, 

//*FORMS=SEP1 


LI 


EXEC FORTVCLG 


//FORT.SYSIN DD ж 


O O O O O O O O 


ооо ооо о ооо о ооо ооо 


ооо 


PROGRAM CVFICA 


AUTHOR: AUGUSTO EZEQUIEL 
DATE: MARCH 17, 1987 


DESCRIPTION: 
THIS PROGRAM CONVERTS THE FICA FILES INTO THE 


INPUT FORMAT OF THE KALMN PROGRAM AND GIVES 


GENERAL INFORMATION ABOUT THE DATA 


DOUBLE PRECISION FPI(500),TP,PB,HD,RLIOFF,RL2OFF 

DOUBLE PRECISION Q1,Q2,PI,C,VOS,K1,K31,K32,TLOLL,TLPLL 
DOUBLE PRECISION PDBPC ,PDBC ,CN1(4),CN2(4),CRI(G) ,CR2(4) 
DOUBLE PRECISION DOP1(4),D0P2(4),BLC,BPDC,SGR1(4),SGR2(4) 
DOUBLE PRECISION SGD1(4),SGD2(4),K2,K4,TT,M; AG,ESQ 

DOUBLE PRECISION TOCS,AS(3),ADC,CRSS,DNS,MOS,CUCS,ES,CUSS 
DOUBLE PRECISION SQAS,TOES;,CICS,OMEGS,CISS,IOS,CRCS,;MWS 
DOUBLE PRECISION DLAT,DLOG,DHT 

DOUBLE PRECISION OMEDS,ADE,;IDOTS,X,Y,Z,ROBL 

REAL LASEC,LOSEC »SEC 

INTEGER INTG(500),B8LOCK, TYPE »NCI,NII ,NFI,LAMIN,LOMIN,LADEG,LODEG 
INTEGER TRACER ,MQVEL(4),ISAT( 4) , INDEX ,SAT , IHKNOS , ТЕОАТ5 »NSAT 
INTEGER IO,STOAY ,MONTH, YEAR,HOUR ; MIN;DAY 

INTEGER STAT(4) 

CHARACTER*8 CHI(500),TITLE(10) 

CHARACTER*5 BLK 

LOGICAL SAVE ,FIRST 


SETS DEFAULT METEO DATA 
DATA TP,PB,HD /15.,980.,75./ 
SETS THE PARAMETERS OF WGS72 
DATA AG, КОВІ /6378.135000,298.26D0/ 


FLAG TO SAVE DATA ONLY WHEN THERE IS NAVIGATION DATA 
FOR ALL TRACERS 


SAVE =. FALSE. 


FLAG TO SAVE INFORMATION DATA AS TYPE 1 ONLY ON FIRST 
CASSETTE 


FIRST-.TRUE. 
ECCENTRICITY SQUARED OF REFERENCE ELLIPSOID 


ESQz(2.D0-( 1. DO/ROBL ) )/ROBL 


91=154.00 
92=120.00 


COHPUTS THE VALUE OF РІ IN DOUBLE PRECISION 
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ооо ооо 


ооо 


ооо 


ann 


G OoOo 


ооо 


PI=DARCOS(-1.00) 

SPEED OF LIGHT ( KM/SEC ) 
С:299792.45800 

NOMINAL SATELITE FREQUENCY 


VOS=10.2306 


CALCULATION OF CONSTANTS TO BE USED IN THE SGIMA COMPUTATIONS 


K1=( C ZVOS )%%2 
К51-(С/(2.00%РІжУ05ж%01) )жег 
К52-(С/(2.00ж%РІ1жУ05х02) )же2 


RL1OFF=-6000.D0 
RL2OFF=7600.D0 


М=%.00 
τιοι].Ξ. 700 
ТСРИ <. 700 
PDBPC=100.D0 
PDBC=100.D0 


SETS THE FLAGS FOR NOT HAVING THE EPHEMIS DATA 


DO 3 I=1,4 
STAT(I)20 


PRINTS OUT THE IDENTIFICATION OF THE PROGRAM 


WRITE(6,2) 
FORMAT(1H1,Z//////// 


В ,/,! зезезежезезезежеежезезезезее жзезежзезжезезезежзезежзезезезезе зе жезезежезезезезезезезежжзеж! 
в »/," зежезежзезезезезезежзезезезежзезезезежезезжезезезежзезезезезезежзезезежзезезезежезе жзезезезеж! 
в оэ“, κ жж! 
ғ ,/,! xx OUTPUT OF THE PROGRAM CVFICA THAT CONVERTS жж! 
$0 κκ THE FICA FILES TO A FORMAT READABLE BY THE жж! 
в ,/,“ xx PROGRAM KALMN2 жж! 
8 И," жж жж! 
8 ,/7/,' XX X 3€ 3900000)0000000( QPX0)900000€X X900 DODO3IOO9OQ(DXYCDDO€ X 30x ' 
8 5,/,' X X39000000€O0€O00000€!0€€IDODXDDOOOX3gM€ X J900O)900900€O9(X X 000€ ' 
B ,Z////////) 


READS THE TITLE FROM THE FILE AND CONTROL INFORMATION 


READ(31,5,END=2000 )TITLE 
WRITE(6,5) TITLE 
READ(31,5,END=2000 )TITLE 
WRITE(6,5,END=2000) TITLE 
FORMAT ( 1048 ) 


READ(31,*,END=2000 ) MONTH,STDAY , YEAR 
ARITE(6,6) MONTH,STDAY , YEAR 
FORMATt///, 
* ' DATE OF STARTING GPS WEEK (MONTH DAY YEAR) ',I3,I3,IB5) 


INPUTS THE DATA FROM THE COMMAND FILE 


READ(31,*,END=2000) NSAT, K2» K& 
READ(31,*,END=2000) X,Y,Z 
МВІТЕ(6,1) М5АТ,Қ2,Қ46,Х,Ү,2 

FORMAT (' NUMBER OF TRACERS ' ;» 1205/7» 

8 ' RANGE SIGMA BIAS FACTOR ",Е20.8,/, 
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' DOPPLER SIGMA BIAS FACTOR ',E20.8,/, 
' ESTIMATE RECEIVER POSITION',/, 

: 4 COORD',F20.6,/, 
қ Y COORD',F20.6;,/;, 
à 2 СООКО" ,Е20.6,/) 


ничия 


C 
C3O0O900000000D09)OOOO)OOOOpCODOO1ORpOOO00)1)0100)00X10900100(O)0 00(0(9000)000090090)0]01010O)0)010)Ó)!PÉ 
C%% INPUT OF A FICA BLOCK xx 


ς ΧΘΘΘΕΧΕΘΘΘΘΟΘΘΘΘΕΘΘΕΘΘΘΕΘΕΘΘΘΕΕΘΘΘΘΕΘΘΕΘΘΘΘΘΘΘΘΘΘΘΘΘΘΘΘΘΘΘΘΕΘΘΘΘΘΘΘΘΘΘΕΘΘΕ 
ς 

С READS THE CONTROL BLOCK INFORMATION 

С 

10 READ(50,20,END=2000) BLK,BLOCK ,NFI ,NII,NCI 

20 FORMAT(A5,415) 


C READS THE FLOATING PART DATA IF ANY 


IF (NFI.EQ.0) GOTO 35 
READ(50,30,END=2000) (FPI(I),I=1,NFI) 
30 FORMAT(4D20.16) 


READS THE INTEGER PART IF ANY 


5 CONTINUE 
IF (NII.EQ.0) GOTO 45 
READ(50,40,END=2000) (INTG(I),I=1,NITI) 
40 FORMAT( 6112) 


READS CHARACTER PART DATA IF ANY 


5 CONTINUE 

IF (NCI.EQ.0) GOTO 55 

READ(50,50,END=2000) (CHI(I),I=1,NCI) 
50 ҒОНКМАТ(10(А8)) 
55 CONT INUE 

IF (BLOCK.NE.6) GOTO 100 

C30000000000000000]000)00)00000000000000)00001080000O0)0)000]8O0O0OO)O0OO)3OOOOOOOO)ÓOO]|)|O)O 
C** TRACING DATA HE 
C 300090000 )9)00000000000)0]9]0909909)8O0O0)00O0O00)OOO)O0OO0O)COOQOOXOOO0DOOOO0O00O0O)O0O0!1OOO000)0OO!|BC 
С 


С USER EPOCH TIME OF PSEUDO RANGE 
C 
TT=FPI(3) 
C 
C COMPUTES THE DATE 
C 


CALL SUB( TT ,DAY »HOUR »MIN,SEC } 
DAY =DAY+STDAY 


IF NO EPHMERIS DATA AVOIOS COMPUTATIONS 


ооо 


IF( .NOT.SAVE) GOTO 10 


L1 CARRIER SIGNAL TO NOISE 


ооо 


СМ1(12-ҒР1(4) 
СМ1(2 )=РРІ(5) 
CN1(3 J=FPI(6) 
CN1(@)=FPI(7) 


L2 CARRIER SIGNAL TO NOISE 


ооо 


СМ2( 1) <ЕРТ( 8) 
CN2(2)=FPI(9) 
CN2(3)=FPI(10) 
CN2(4)=FPI(11) 


Ll PSEUDO RANGE (KM IN FICA FILES CONVERTED TO SECONDS FOR KALMN2) 


ооо 


CR1(1)-FPI(12)/C 
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CR1(2)=FPI(13)/C 
CR1(3)=FPI(14)/C 
CR1(4)=FPI(15)/C 


PSEUDO RANGE (KM IN FICA FILES CONVERTED TO SECONDS FOR KALMN2) 


CR2(1)=FPI(16)/C 
CR2(2)=FPI(17)/C 
Са2(5 )=ЕРІ( 18 )/С 
CR2(4 )-FPI(19)/C 


CARRIER DOPLER PHASE 


DOP1(1)=FPI( 20) 
DOP1(2)=FPI( 21) 
DOP1(3)=FPI( 22) 
DOP1(4 )=FPI( 23) 


CARRIER DOPLER PHASE 


DOP2(1 J=FPI( 24) 
DOP2( 2 )J=FPI( 25) 
DOP2(3 )=FPI( 26) 
DOP2(4)=FPI( 27) 


PRN OF EACH TRACER 

ISAT(1)=INTG(1) 

ISAT( 2 )=INTG( 2) 

ISAT(3 )=INTG(3 ) 

ISAT(4 J=INTG(4) 

(2 QUALITY FACTOR ( TRACER,FREQUENCY ) 
MQVEL(1)=INTG(9)-INTG(13) 
MQVEL(2)=INTG(10)-INTG( 14) 

MQVEL(3 )=INTG(11)-INTG( 15) 

MQVEL( 4 )Z INTG( 12) -INTG( 16) 

ANY ERRORS FLAGS THE DATA OF THE SATELLITE 


DO 120 INDEX=1,NSAT 


ERRORS IN SIGNAL 


IF( MQVEL(INDEX).EQ.0) GOTO 114% 
MQVEL( INDEX )=15 
I1=3+INDEX 
I2=12+INDEX 
WRITE(6,110) MONTH,DAY,YEAR,HOUR »MIN»SEC, INDEX, INTG(I1),INTG(I2) 


110 FORMAT(/;' 3e9«-30(*-** WARNING *%-%-% AT ', 
*  2(I2,1X),I4,1X,2(I2,':'),F5.25/; 


' DATA FROM TRACER ', 
12," СЕ5АН BAD STATUS ",15,1Х,15) 


MQVEL( INDEX )=15 
THE DATA IS ALREADY BAD NO NEED TO TEST THE SATELLITE 


GOTO 120 


NON EXISTENT SATELITE DATA FOR THAT TRACER 


11% CONTINUE 


C 
С Бе 
С 
€ 
C L1 
C 
C 
C Ба 
C 
C 
C SV 
G 
C 
С L1; 
C 
C 
C ТЕ 
C 
G 
C 
C 
3 
= 
C 
C 
C IF 
C 
C 
C 
C 
C 
C 


Е 


IF( STAT(INDEX).EQ.ISAT(INDEX)) GO ТО 120 


WRITE(6,115) MONTH,DAY,YEAR,HOUR ,MIN,SEC , INDEX, 
ISATC INDEX ),STAT( INDEX ) 


115 FORMAT (/, ' жж ж κ.κ WARNING жж ж жж AT ', 
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я 2Х512001Х%,12,1Х,16,2%,12,!:7;12,!:%,ББ.2,/, 
в"! TRACKING TRACER ',I2,' SV ',I2, 
8" EXISTING NAV DATA FOR SV ',12,' FLAGED BAD QUALITY DATA') 
MQVEL(INDEX)215 
120 CONTINUE 


C 
C COMPUTS THE SIGMAS 
C 


DO 15D INDEX=1,NSAT 
IF( MQVELCINDEX).NE.D) СОТО 150 


C 
C CALCULATE L1 AND L2 SIGMA RANGE 
C 
BLC-TLDLL 
BPDC-PDBPC 
SGR1(INDEX)z((BLCXM)Z(C1D.39( CN1CINDEX)Z71D. ))) 
я ж {. БЖВРОСИ( 10 . жж( СМ1 (ТМОЕХ)/10.)) J*K14K2 
SGR2(INDEX)-((BLCXM)/(C1D.39( CN2C( INDEX)Z1D. ))) 
8 ж (.БХВРОС/(10. жж( СМг2( ТМОЕХ)/10.))  К1 +К2 
С 
С CALCULATE L1 AND L2 SIGMA DOPPLER DATA 
E 
BLC=TLPLL 
BPDC=PDBC 
SGDIC INDEX )=( (BLC#M)/01D. €(CN1( INDEX )/10.))) 
я ж 11. +(.5Б#ВРОС/110.жж( СМ1(ТМОЕХ)/10.))) «31 +Ка 
SGD2( INDEX )=( (BLC*M)/(1D. **(CN2( INDEX )/1D.))) 
5 X (1.*C. 5*XBPDC/(1D.**(C CN2C INDEX)/71D. )) ) )*«K 32 4K 
15D CONTINUE 
Ë 
C OUTPUT THE DATA TO THE FILE IF THERE IS NAV DATA 
C 
TYPE=3 
WRITE(7) TYPE 
WRITE(7) TT, ISAT, CR1,» CR2, DOP1, DOPZ, 
= SGR1, SGR2, SGD1, SGD2, MQVEL 
Є 
С GO TO READ ANOTHER BLOCK 
C 
GOTO 1D 
Е 


100 СОМТІМЈЕ 
IF (BLOCK.NE.1D9) GOTO 200 
C EIEE JE HEIE JEE FE FEJE JEJE JE JE JEJE JE IE JE IE JE JE JE JE JEJE IE IE JE JEJE JE IE JEJE HEIE JE JÉ © ЭЄЗЕ ЗЕ ЭЕ ЭЕЄЭЕ-ЭЕ ЭЕ ЭЕ ЭЕ © ЭЕ ЭЄ ЭЕ JE JE JE EIE MEJEJE JEFE FE FE 
C¥* NAVIGATION DATA (AS TRANSMITED) жж 
Q OCIO 5% 555555555556 356% FE IEE Без 
С 
TRACER-INTG(1) 


SAT :INTG(2) 
С 
E GO TO READ ANOTHER BLOCK 
E 
GOTO 1D 
Е 


200 CONTINUE 
IF (BLOCK.NE.9) GOTO 500 


C 
отч | ———— XOOCOUCXOUDOIODOOOOIO I XX CK XO ODIO XOOOOC X XOOODOXOOOOp(ee eer 
Сжж NAVIGATION DATA (DECODED SUBFRAMES 1 TO 3) жж 


Схжкиххижжжжажжаежежзезезезезезезежзезезезезезежзезезезезе зе зезезезезежзезезежзежжзежзезезвезезезезезезжзезезезе ж 


Е DAY OF WEEK 


IWKNOS=DINT(FPI(6)) 
Сжж IWKNOS =IAND( IWKNOS »17778 ) 


C 
C SV HEALTH 
C 
IEDATSzDINT(FPI(9)) 
C 
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оо О 


ооо оо О ann ооо 


aan 


CLOCK EPOCH (GPS SECONDS OF WEEK) 
TOCS=FPI(13) 

CLOCK BIAS (SEC) 
AS(1)=FPI(16) 

CLOCK DRIFT ( SEC/SEC) 
AS(2 )=FPI(15) 

CLOCK DRIFT RATE (SEC/SEC*2) 
AS(3 )=FPI(14¢) 

AGE OF DATA (CLOCK) (SEC) 
ΑΌςΞΕΡΙ(10) 

RADIAL SINE CORRECTION (DIVIDED BY 1000. TO GET KM) 
CRSS=FPI(27)/1000. 

CORRECTION TO MEAN MOTION (RADIANS/SECONDS ) 
ONS=FPI( 28) 

MEAN ANOMALY AT EPOCH ( RADIANS ) 
MOS=FPI(29) 

IN TRACK COSINE AMPLITUDE (RADIANS) 
CUCS=FPI(30) 

ECCENTRICITY 
ES=FPI(31) 

IN TRACK SINE AMPLITUDE (RADIANS) 
CUSS=FPI(32) 

SQUARE ROOT OF SEMI-MAJOR AXIS CONVERTED TO SQR OF KM 
SQAS=FPI(33 )¥FPI(33) 
SQAS=SQAS/1000.00 
SQAS=DSQRT(SQAS ) 

TIME OF cPOCH (GPS SECONDS OF WEEK) 
TOES=FPI( 34) 

INCLINATION COSINE CORRECTION ( RADIANS) 
CICS=FPI(%6) 

RIGHT ASCENSION NODE ( RADIANS ) 
OMEGS=FPI( 47) 

INCLINATION SINE CORRECTION (RADIANS) 
CISS=FPI(48) 

INCLINATION (RADIANS) 


IOS=FPI(%9) 
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C RADIAL COSINE ADJUSTMENT (DIVIDED BY 1000 TO GET KM ) 


CRCS=FPI(50)/1000. 


Š ARGUMENT OF PERIGEE ( RADIANS ) 

: WS=FPI(51) 

с RIGHT ASCENSION OF ASCENDING NODE (TIME DERIVATIVE) 

C ( RADI ANS/SEC ) 
OMEDS=FPI(52) 

E AGE OF DATA ( SEC ) 

i ADE=FPI( 26) 

s INCLINATION TIME DERIVATIVE ( RADI ANS/SECONDS ) 

: IDOTS=FPI( 54) 

е SETS THE FLAG FOR EXISTANCE OF NAV DATA FOR THIS SATELITTE 

; STAT( TRACER J=SAT 

с TESTS IF THERE IS ENOUGH NAV DATA 


SAVE=.TRUE. 

DO 210 INDEX=1,NSAT 

IF(STATCINDEX).EQ.0) SAVEz.FALSE. 
210 CONTINUE 


C 
C OUTPUTS THE DATA 
C 
TYPE=2 
WRITE(7) TYPE 
WRITE(7) TRACER,SAT,AS(1),AS(2),AS(3),CICS,CISS 
© Í| CRCS, ,CRSS,);CUCS,CUSS,DNS, ES, IDOTS,IOS ,HMOS 
z »OMEDS ,OMEGS ,SQAS ,ADE , TOCS, TOES, INKNOS ,WS 
т »LEDATS,ADC 
С 
C INFORMS THAT NAVIGATION DATA WAS RECEIVED FOR THIS SATELLITE 
С 
С 


WRITE(6,220) TRACER,SAT,MONTH,DAY ,YEAR,HOUR,MIN,SEC 
220 FORMAT(///,' зезезезезезезезезезезезежзезезезежезезе жзежезезе зезе зе зе жзежзезезжежезезеж", // 


3 y X% NAVIGATIONAL жж",/, 
5 ' X* DATA FROM TRACER ',I2,' SV ',I3,' RECEIVED **',/, 
8 ж” >7Х,12,1х,12,1Хх,16,2Х,12, het >І2, Tag »F5.257X» жж! 7» 
5 ! зезезезезезезезезезезезезеж зезезезезезезезе зежзезе жзежзе же жени", /// ) 
С 
E PRINTS OUT THE WARNING WHEN THE HEALTH IS NOT GOOD 
Е 
IF( FPI(9) .NE. 0.0 ) WRITE(6,201) TRACER,SAT,FPI(9) 
201 FORMAT (AAA , | RRR IE IEEE FEE HEHE AE 9€ EJ€XE 9€9€ 9€ 3€ € 3€9€3€ 3€ 3€ 9€ HEE ' 
#,' TRACER ',I3,' SV ',I3,' HEALTH STATUS ' ,F15.6, 
%2 „ "Зезезезезе зезезезезезезезе же зезезезезезезезезезе зезезезе жзезжзезежзезееж" ) 
C 
C GO TO READ ANOTHER BLOCK 
C 
GOTO 10 
C 


200 CONTINUE 
IF (BLOCK.NE.101) GOTO 400 


C 

C30 OO COOXOCXOUOOOEOOIXDOOOOUODOOO)O0O0O)000O000000000600)0O09)1900)00O)00)0)0)9))OO|BO/B|O)!|COO!|/HOC 
Сжж BLOCK 1 PLUS INPUT DATA xx 
C3 3€ 3€ JE 3€ 36 3€ 369€ € € 9€ JE JE FE JE JE JEJE JE JEJE EXEXE CE OOOOOOXE 3OEXCCGOEOXOOXOOE(OXOOEOXOOOO)OOOO00O0)00)0€() FE X HE 
C 

C PRESSURE 

C 
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ооо aan 


ооо 


ооо 


ооо 


O00 


C 
400 


C 


PB=FPI(1) 
TEMPERATURE 

TP=FPI( 2) 
HUMIDITY 

HD=FPI(3) 


TESTE FOR BAD WEATHER DATA. CHANGE TO STANDARD VALUES BAD DATA 


IF (( PB .LT. 800. ) .OR. ( PB „СТ. 1200. )) РВ = 980. 
IF (€ TP .LT. -99. ) .OR. € TP .GT. +995 J) ТР ΠΕ, 
IF (( HD .LT. 1. ) ОВ. ( HD .GT. 100. )) HD = 75. 


OUTPUTS THE DATA FOR THE FIRST TIME 
IF(.NOT.FIRST) GOTO 10 
TYPE=1 
WRITE( 7) TYPE 
WNRITEC7) RLIOFF;,RL2OFF;PB,TP,HD,X,Y,2Z 
SETS THE FLAG 
FIRST=.FALSE. 
GO TO READ ANOTHER BLOCK 
GOTO 10 


CONTINUE 
IF (BLOCK.NE.3) GOTO 460 


C36 36 FE HE IESE FE DE FE IE IE IE DE IE IE IE IE E E IE FE IEE SEI JE JE ME JE E IEEE EIEEE IESE DE FE DE 93)909DDE EODD 3€XOOOOOOO9OO0O0O03 


Сжж 


GESAR SOLUTION жж 


C 3€3€3€ € € 33€ 3€ 3€ € 9€ X€ 3€ 3€ 9€3€ 9€ 3€ 3€ 3€ 9€39€ 9€ € 93€ 9€ 9€9€3€ CCXE EXE EXEXEXY€ 56 562 65556556565656565655565556565<5<3655565656563656 


C 
C 
C 


оо Оо 


+G O O O aan ооо 


aan 


SOLUTION USER EPOCH TIME 
TT=FPI(2) 
IF NO EPHMERIS DATA DOES NOT KEEP THE DATA 
IF(.NOT.SAVE) GOTO 10 
COORDINATES OF RECEIVER 
X=FPI(3) 
Y=FPI(4) 
ΖΞΕΡΙ(5) 
COMPUTES ΤΗΕ DATA 


CALL SUB(TT,DAY,HOUR ,MIN,SEC ) 
DAY=DAY+STDAY 


ONLY PRINTS POSITIONS IF UNCOMENTED 
WRITE(C6,450) MONTH,)DAY,YEAR,HOUR,MIN,)SEC,X,Y Z 
FORMAT (7/7, ' GESAR SOLUTION AT (MH DAY YYYY HH:MM:SS.SSS)', 
# 2X,12,1X,I12,1X,16,2X,12,': I2, ' : ' ;F5.2, 
#/,' PSN COORD X, Y, Z (KM) ',F14.3,' ! ;F16.3,' ',F14.3) 
CONVERTS TO LATITUDE / LATITUDE / HT 
CALL XYLLH(X,Y,Z,AG,ESQ,DLAT »DLOG,DHT ,PI) 


PRINTS THE RESULTS 
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о О 


HRITE(L6,455) DLAT,DLOG,DHT 
«455 FORMAT(' LAT, LON, HT (DEG, DEG, KM) ',2F17.12,2X,F12.7) 


C 
E CONVERTS LATITUDE TO DEGREES MINUTES AND SECONDS 
C 

CALL DMS(DLAT,LADEG,LAMIN,LASEC) 
C 
С CONVERTS LONGITUDE TO DEGREES MINUTES AND SECONDS 
C IN A RANGE O TO 180 

IF (OLOG.GT.180.D0) DLOG=360.00-DLOG 

CALL OMS(DLOG,LODEG,LOMIN,LOSEC ) 
C 
C WRITE(6,456) LADEG,LAMIN,LASEC ,LODEG,LOMIN,LOSEC 

456 FORMAT(' LAT, LON (DEG MIN SEC, DEG MIN SEG) ',2(214,1X,F7.4),/} 

С 
С OUTPUTS THE DATA TO A FILE 
E 

WRITE(8) MONTH;DAY , YEAR, HOUR MIN,;SEC, 

т LADEG,LAMIN,LASEC ,LODEG,LOMIN,LOSEC 

C 
C 
C GO TO READ ANOTHER BLOCK 
C 

GOTO 10 
C 


460 CONTINUE 
IF (BLOCK.NE.13) GOTO 470 


C 
C3O0O000C(OOODXOOOOODOXOOXOOOOOOOOOOIPOOUOOODOODOOOOOOOOOOOOOOODXOOOOOOOOOO[!|9POO[OO!OG 9e 
Сжж ТАРЕ HEADER/TRAILOR жж 


C EE SEHERE REE HE JEJE JEE FEFE FE H JEJE JE JE IE FE IDEJE JEE JE JEJE JE JE JE IE JE IE JE JEJE FE IE JE JE JE JE JEJE FE HEIE JE ME EIE JE JE JEE ME FE JE EJEDE JE JE JEME FE 
С 
WRITE( 6,461) 
461 FORMAT ( 7 9 5X 5 ' 336363636 36363636 266 HIE IE IE IE HE FE HE IEIENE HE HEHEHE ! ,7, 
т 5X,'** TAPE HEADER/TRAILOR  3se';/, 
z 5X,  X399000000000000000000)0)00)00)0€€ ! ,7) 
IF( INTG(1).EQ.1) WRITE (6,462) 
IF( INTG(1).EQ.2) WRITE (6,463) 
ІҒ( INTG(1).EQ.4) WRITE (6,464) 
ТЕ: INTG(1).EQ.8) WRITE (6,465) 


G62 FORMATI ' BEGINING OF DATA SET',/) 
463 FORMATI ' END OF DATA SET',/) 
G64  FORMAT(' BEGINING OF CASSETTE',/) 
465 FORMATI ' END OF CASSETTE',/) 
WRITE(6,466) INTG( 2); INTG( 3) 
466  FORMAT(' CASSETTE ',I5,' IN DATA SET ,/, 
3 ' CASSETTE SEQUENCE NUMBER ',I5) 
C 
С GO TO READ ANOTHER BLOCK 
C 
GOTO 10 
C 


470 CONTINUE 
IF (BLOCK.NE.11) GOTO 480 


С 
E IF NO EPHMERIS DATA DOES NOT KEEP THE DATA SO NO NEED TO 
С INFORM THE USER 
C 
IF(.NOT.SAVE) GOTO 10 
С 
C XOOUOGCOOCOCOOXOIOOODOOOOODOOOOOOOHIOOOOOODOOOOOOOOOOOOOOOOOOOOOO!OOOOOOOOO|OOO 
cxx RECEIVER ERROR BLOCK жж 
C EEIEIEE IE E HE IEE HEIE JE HEIE JEE JE HE JE HE IE JE FE JEJEJE JE JE JE JEJEJEJE JE JE EME IE JE JEJEJE HE JE FEJE HEFE HEIE JE JE EIEH IAAI SAAB 
C 
C 
АВТТЕ( 6,471) MONTH >,DAY,YEAR,HOUR ,»MIN,SEC 
471 FORMAT(/,;' X33 OCOCOOOOOOCXOCOOOGOOCEEOX OOOOOOOOO0O0O0O0O0O0000)060€ ! 7, 
+ : жж ERROR BLOCK RECEIVED AFTER : жж",/, 
т я 33€' ,9X,2( I2, 1X), I6, 2X, 2(I2, ' : ' );)F5.2, 8X, ЖЖ",/, 
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2 : HEE O)(O0tO)10O0O$0)I/3 oO0)QO(M(PO)XMUDOOQXOO9OQ(OQO(OOX»»€9 ' , / ) 


WRITE (6,%72)(INTG(I),I=1,NII) 

472  FORMAT(' ERROR LOG MESSAGE TI ICO FORMAT $ 1',18,/, 

г" ,18,/, 
5! ,18,/, 
4" „18, /, 
5' )18,/, 
6",18,/, 
г" :185/7, 
8',I18,/, 
ERROR LOG OVERFLOW COUNT ' »185,7) 


HHBARAH 


12g: HV #H 


GO TO READ ANOTHER BLOCK 


ооо 


сото 10 


480 CONTINUE 
IF (BLOCK.NE.8) GOTO 500 


C 
С IF NO EPHMERIS DATA DOES NOT KEEP THE DATA SO NO NEED TO 
С ΙΝΕΟΒΗ THE USER 
C 
IF(.NOT.SAVE) GOTO 10 
C 
CJO9O9OO0OO009000)00DOOOOOODOHO|OIDIM/OOOOOOOOPIOOIO]IO|IOOPOOOÓOO|OO/OOÓOOUCOOOHOOIOOHO!PHOGO|HO(POC 
οκ TRACKING CONFIGURATION жж 


С 565555565 JE E E E IE JEJE HE JE HE IE Э 63636553655 5-5-%-5 55 55555555555 56555556556 55556553655 653655% 
C 
С COMPUTES THE DATE 
С 
CALL SUB(FPI(2),DAY ;HOUR,MIN,;SEC) 
DAY-DAY «5ΤΌΑΥ 


C 
C 
МКІТЕ(6,981) MONTH,DAY YEAR HOUR ;MIN,SEC 
481 ЕОНМАТ(/»! Х5555%5555 555555 55%5<553656555 5565 55656-55 ЭКЕ" ,/, 
5 жж TRACKING CONFIGURATION AT жж",/, 
= я же" ,9х,2(12,1х),14,2х,2(Т12,":"),Е5.2,ВХ, жж",/, 
5 : 55-5 5-5 555555 COE JCXOE eX COCX EXEC XOIEOXOEXOOOO( X OOO(!O90(*099 ! , 7) 
С 
WRITE (6,482) FPI(1),(FPI(I),I=3,NFI),(INTG(I),I=1,NII) 
482 FORMAT{' PSEUDORANGE FTF OF VALIDITY SEC ',Е22.1%›/, 
* ' PREDETECTION BANDHITHS HZ FREQUENCY ,TRACER',/, 
В 2(9Е20.10,/),//, 
# ' TRACKING LOOP BANDWITHS HZ FREQUENCY , TRACER' ,⁄ 
5 9( 8Е20.19,/),/» 
# ' LOOP ROUND ТВІР СА(ТВВАСТОМ ΡΕΙΑΥ5 ТМ ЗЕСОМО5 ",/, 
в" INCLUDING PRE-ANTENNA',/,' 11 ',E20.154,' L2 ',E20.15,//, 
8 ' LOOP ROUND TRIP CALIBRACION DELAYS IN SECONDS ',>/, 
* ' (INTERNAL TO RFM)',;Z,)' Ι 1 ",Ег2О.14," |2 ",Е2О.14,//) 
+ ' SV PRN ID QF EACH TRACER ",/,4115,//, 
8 ' ANTENA INDICATOR/TRACER' >/,41I155//, 
з ' CODE INDICATOR FLAG/TRACER' »/,41I15>//) 
C 
C TESTE IF THE CONSTELATION MATCHES THE EXISTING EPHMERIS 
С DATA 


DO 495 INDEX=1,NSAT 
IF(STAT{ INDEX) .EQ.INTG{INDEX)) GO TO 495 
WRITE(6,490) INDEX,STAT{ INDEX) ,INTG( INDEX) 
490 ОРОВМАТ(////," αχ 323333 WARNING 4363606360006 жзвевеж! ,//, 
# ' ACTUAL NAVIGATIONAL DATA FOR TRACER ',I2,' IS FROM SV ',I2»/, 
# ' TRACKING CONFIGURATION SHOWS SV ',I2,//) 
495 CONTINUE 


С 
E GO TO READ ANOTHER BLOCK 
C 
GOTO 10 
C 


500 CONTINUE 
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C HEE EEIEIEE JEJEJE E JEJERE JEJE EEPE HEIE JEME JE JE JE JE JE E JE JE JE EJE JE JE JEE HE JEJE FE JE JEME JE JE JE JEJE JE JEJEJEJE JE HE X X OX» x 
Сжж BLOCKS THAT THIS PROGRAM IS NOT PREPERED TO DECODE жж 
C IE EIEH IEIR EIEEEI IEEE JE JEFE JEJE JEE JE HE JEDE E JEE EJE JE JEME XXX € EX € X E HE JE HE JE JEJE JEME JEME JEME FE FE 2E FE 
C 
WRITE(6,501) 
501 FORMAT ( // , 3HG6O)SOOOODOODOODOOPOODOODOOOOODIOOOCOOOOOOOOOOODIOHPOOOp|OO|Oeee ! , 
# Z,! κκ FOUND A BLOCK THAT THIS PROGRAM IS NOT жж", 
# /»' жж ABLE TO DECODE . BLOCK DATA IS DUMPED NEXT **', 
8 /,' ΘΘΘΘΘΘΘΚΘΕΚΚΧΘΘΘΕΘΕΟΘΘΕΕΘΟΘΘΘΘΘΘΘΘΘΟΘΘΘΕΘΘΟΘΘΘΘΘΕΚΘΚΕΙ ) 
WHRITE(6,520) BLK,BLOCK,NFI,;NII,;NCI 
IF (NFI.NE.O) WRITE( 6,530) (FPI(I),I=1,NFI) 
IF (NII.NE.0) WRITE( 6,540) (INTG(I),I=1,NII) 
IF (NCI.NE.0) HRITE(6,550) (CHI(I),I=1,NCI)} 
520 FORMAT ( 1X , A5 ,GI5) 
530 FORMAT(1X,4E20.194) 
540 FORMAT(1X,6I12) 
550 FORMAT(1X,10(A8)) 


C 
C GO TO READ ANOTHER BLOCK 
C 
GOTO 10 
C 


Сжжжзезеж жзежаеж жзежзежезезнезезезезежезезе жеззезеззе зезезезезе жаезезежзезежжзежзезезежзезезезежжзезезезезезезезе жзежзезе 
C 


C END OF PROGRAM 
C 
2000 CONTINUE 

STOP 

END 
C 
C3OOODO00000000)0)0O0O00DODDOOOOPDOODCCHODQPQOOQDOOOQQOQOIOBOODODOOUOOOIOOCOOODPOOOIHODOOODO[|P]ObOee 
E THIS SUBROUTINE WAS TAKEN FROM THE KALMN PROGRAM 


Сенна IE IE E ME IEE IE IE C Me XC ХХХ ХХХ ХХ << ХХХ Ж-Е 
С 
SUBROUTINE XYLLH 

G ( Хх › Ү › 2 › АС 5» ESQ , 

R DLAT » DLOG , DHT,PI) 
THIS SUBROUTINE CONVERTS POSITION FROM THE EARTH-CENTERED CARTESIAN 
COORDINATE SYSTEH TO THE GEODETIC COORDINATE SYSTEM. 
X,Y,Z-ESTIMATE OF RECEIVER POSITION IN KM 
AG-SEMIMAJOR AXIS OF REF. ELLIPSOID(RADIUS OF EARTH) IN KM 
ESQ=ECCENTRICITY SQUARED OF REF. ELLIPSOID=(2.-1./0BL )/OBL 
DLAT=GEODETIC LATITUDE OF RECEIVER POSITION IN DEGREES 
PLOG-GEODETIC LONGITUDE OF RECEIVER POSITION IN DEGREES 
DHT=GEODETIC HEIGHT OF RECEIVER POSITION IN METERS 
PI VALUE OF PI COMPUTED IN DOUBLE PRECISION IN THE MAIN PROGRAM 
SINCE LAT=LAT(X,Y,Z,LAT) AND HT=HT(X,Y,LAT), AN ITERATIVE PROCEDURE 
IS NECESSARY TO DETERMINE THE VALUES OF LATITUDE AND HT. 
THE WORKING EQUATIONS ARE AS FOLLOWS: 
HT=(R/COS(LAT))-AG1 WHERE AG1=AG/(SQRT(1-ESQ¥SIN( LAT )32)) 
ТАМ: 106) <У/Х 
TANI LAT J=((0Z4( ESQ¥SIN( LAT )¥AG1)/R) 


оээоооооооооооооо 


DOUBLE PRECISION LT , LT1 , LG > R › ЕР ›АС1 
DOUBLE PRECISION X ; Y » Z » AG , ESQ , DLAT , DLOG , DHT , PI 


С 

С FIND LONGITUDE 
LG » DATAN2(Y X) 
IF (LG .LT. 0.DO) LG - LG * (2. * PI) 
RzDSQRT ( x3oex2 «Y»ex2 ) 


A FIRST GUESS FOR LATITUDE WOULD BE ATAN(Z/R) 
LT1-DATAN(ZZR ) 


оо 


С SOLVE FOR LATITUDE BY ITERATIONS 
DO 10 I=1,5 
AG1=AG/( DSQRT(1.D0-ESQ*XDSIN(LT1 )#*2)) 
LT=DATAN( (Z+( ESQXDSIN(LT1 )*AG1) /R) 
EP=DABS(LT-LT1) 
IF(EP.LE. .000000001D0) GOTO 5 


J5 


LT1 = ЕТ 
10 CONTINUE 
5 AG1=AG/(DSQRT( 1D0-ESQ*DSIN(LT1)%%2)) 
C 
С FIND HEIGHT 
DHT-( (t R/DCOS(LT ) ) -AG1) 
C CONVERT DLAT AND DLOG FROM RADIANS TO DEGREES 
DLAT-LT*180.DO/PI 
0106:16#180.00/РТ 
RETURN 
END 
C 
(Lx x on MOD Mae EO Eo de sp IE IE OEE EERE IEE IE SE FE ENE IEE BEE IEE IOI KC oe J AE AE FE KAE FEFEFE IE IHR C ic Meo 
C THIS SUBROUTINE CONVERTS DECIMAL DEGREES IN DEGREES MINUTES SECONDS 
RE EHH IE EE HIE HE FEE IE A JE E DE эе ЭЁ Ж ЭЕ ЭЕ ЭЕ ЭЕ ЭЕ ЖЄ Ж JE E HEIE JE HE FE JE E IE FE FE FEJE FE IE FE PE IEEE HEH JE HE IIE IEICE 
C 
SUBROUTINE DMS(DEC,;LDEG,;MIN,SEC) 
C 
DOUBLE PRECISION DEC,XNUM,XMIN 
C 
XNUM=DABS{ DEC ) 
LDEG=DINT { XNUM ) 
XMIN=( XNUM-DFLOAT(LDEG ) )*60.D0 
MIN-DINT ( XMIN) 
XNUM* t XMIN-DFLOAT ( MIN ))*60.D0 
SEC zSNGL ( XNUM ) 
IF (LDEG.GE. 360) LDEG-LDEG-360 
RETURN 
END 
C 
Ih oe Je ο Хе 55% ο ο κκ. 
C THIS SUBROUTINE CONVERTS TIME TAG IN DAYS HOURS MINUTES SECONDS 
LR KH RRR HK KKH EHR HR HEH IEEE EHH SEI KHALI EOI OIE FOIE FE FOIE 
C 
SUBROUTINE SUB(TT, DAY, ,HOUR, MIN ,)SEC ) 
INTEGER DAY, HOUR, MIN 
DOUBLE PRECISION TIME,TT 
REAL SEC 
TIME =TT 
DAY = DINT{ TIME / {2¢.D0*3600.D0)) 


TIME = TIME - DFLOAT(DAY) * 24¢.D0 * 3600.D0 
HOUR = DINT{TIME / 5600.00) 
TIME = TIME - DFLOAT{ HOUR) ж 5600.00 


MIN = DINT(TIME / 60.D0) 

TIME = TIME - DFLOAT{MIN) ж 60.00 
SEC = SNGL( TIME ) 

RETURN 

END 


//GO.FTO6FOO1 DD SYSOUT=* 

//GO.FTO7FOOl DD DSN=MSS.S0812.CVFICA.OUTPUT ,DISP=SHR 
//GO.FTO8FOO1 DD DSN=MSS.S0812.GESAR.POS,DISP=SHR 
//GO.FT31F001 DO DSN=MSS .S0812.CVDATA,DISP=SHR 

//60.ҒТ50Ғ001 DD DSN=MSS.S0812.GPS.TAPES.DATA{DL285),DISP=SHR 


С DD DSN-MSS.S0812.GPS.TAPES.DATA(DL 286 ) DISPzZSHR 
// DD DSN=MSS.S0812.GPS.TAPES.DATA({ DL286 ),DISP=SHR 
// DD DSN-MSS.50812.GPS.TAPES.DATA(DL287) ; DISP-SHR 
// DD DSN-MSS.S0812.GPS.TAPES.DATA(DL 288 ) , DISP-SHR 
Z Z DD DSN=MSS.S0812.GPS.TAPES.DATA({D0L289 ),DISP=SHR 
// DD DSN=MSS.S0812.GPS.TAPES.DATA({DL290 ) ,DISP=SHR 
77 DD DSN:MSS.S0812.GPS.TAPES.DATA(DL291) DISP -SHR 
А DD DSN-MSS.S0812.GPS.TAPES.DATA(DL292) ,DISPzSHR 
// 
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APPENDIX F 
PROGRAM TRANSDUC GPS. SOURCE LISTING 


// JOBTRANS JOB (D812,9999),'EZEQUIEL ',CLASS-C 

/ /*MAIN ORG=NPGVM1.0812P ,LINES=( 99) ,CARDS=1 99) 
//*FORMAT PR,DDNAME -GO.FTD6FDO1, 

//7*FORMS=SEP1 


A 


EXEC FORTVCLG 


//FORT.SYSIN OD * 


OQ O O O OO O O GO O O GO O OO O O OO O GO O O OO O O O G G O O GO О 


ооо 


ооо 


PROGRAM TRANSDUC GPS 


AUTHOR: AUGUSTO EZEQUIEL 
DATE: MARCH 3D, 1987 


DESCRIPTION: 


THIS PROGRAMS TAKES EACH POSITION OF THE ANTENNA, 

THE COMPUTED COURSE ( FILTERED OR NOT ) ,THE PITCH 

AND ROLL DATA, AND COMPUTES THE POSITIONS OF THE 
TRANSDUCER, APLYING THE SEVEN PARAMETER TRANSFORMATION 

TO THE OFFSETS OF THE TRANSDUCER IN RELATION TO THE ANTENNA 
CONSIDERED THE CENTER OF COORDINATE SYSTEM AND THEN ADDING 
THE CORRECTIONS TO THE COORDINATES. 


THE PROGRAM RUNS IN MVS 


I/O SPECIFICATIONS: SEE END OF THIS JOB. 


THIS PROGRAM IS RESTRICT AS IS TO THE SEAFLOOR BENCHMARK 
EXPERIMENT 


THE PROGRAM WILL RUN WITH ANY AMOUNT OF DATA ONLY LIMITED TO DISK 
SPACE 


ANY BLANK LINES WILL TERMINATE THE PROGRAM IN ERROR 


DOUBLE PRECISION XPOS,YPOS, TIME1,TIME2,TIME3,SECD,SECH,RATE 
DOUBLE PRECISION TIME4, TIMES 

DOUBLE PRECISION TT,DLAT,DLOG,DHT>X,Y,2 

REAL COURSE ,PITCH( 2),ROLL( 2) ,SEC1(2),SEC2,PIT,ROL,SEC3( 2) 
REAL DX,DY,DZ,OFFX,OFFY ,OFFZ,HEAD( 2),LASEG, LOSEG 

INTEGER MONTH1(2),DAY1(2),YEAR1(2),HOUR1(2),MINI1(2) 
INTEGER MONTH2 »DAY2 » YEAR2 ; HOUR2 »MIN2 

INTEGER MONTH2(2),DAY3(2 ) , YEAR3( 2) HOURS ( 2) ;MINS(2) 
INTEGER LADEG,LAMIN,LODEG,LOMIN, STDAY , NSAT 

CHARACTER*8 TITLE(1D) 


INITIALIZATION OF CONSTANTS 


SECH*3600.DD 
SECD=24.DD*SECH 


OFFX=-4.268 
OPPRYS-2.955 
OFFZ=-9.502 
READS THE TITLE FROM THE FILE 
READ(3D0,5,END=10D )TITLE 


WRITE(6,5) TITLE 
READ(30,5,END=10D )TITLE 
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an 


ооо 


eo ооо ооо ооо оо О 


ооо О 


оо оО 


aan 


5 


WRITE(6,5) TITLE 
FORMAT(10A8 ) 
READ( 30 ,* ;,ENDz2100) MONTH2,;STDAY , YEAR2 
WRITE(6,6) MONTH2,STDAY , YEAR 
6 FORMAT (///; 
8 " DATE OF STARTING GPS WEEK (MONTH DAY YEAR) ‘,13,13,15) 


READS TWO SETS OF PITCH AND ROLL DATA 


READ(51,%*,END=100) MONTH1(1),DAY1(1),YEAR1(1),HOUR1(1),MIN1(1), 


$ SEC1(1),PITCH(1),ROLL(1) 
TIMEL=DFLOAT( DAY1( 1) )*¥SECD+DFLOAT( HOURI( 1) )*SECH 

$ * DFLOAT(MIN1(1))*60.DO«*DBLE(CSEC1( 1)) 
READ(51;,*,END-100) MONTH1(2),DAY1( 2), YEAR1(2) ;HOUR1( 2) ; MIN1( 2); 

$ SEC1(2),PITCH(2),ROLL( 2) 
TIME2zDFLOAT( DAY1( 2 ) JXSECD*DFLOAT ( HOUR1( 2 ) )*SECH 

$ * DFLOAT(MIN1(2))*60.DO0*DBLE( SEC1( 2)) 


READS THO SETS OF COURSE 


READ( 52, ENDz100) MONTH3(1),DAY3(1),YEARS( 1)  HOURS( 1) ; MIN3(1), 


$ SEC3(1),HEAD(1) 
TIME¢=DFLOAT (DAY3( 1) J¥SECD+DFLOAT( HOUR3( 1) )*¥SECH 
$ + DFLOAT(MIN3( 1) )*60.D0+DBLE( SEC3(1)) 
READ( 52,END=100) MONTH3(2),DAY3(2),YEAR3( 2) »HOUR3(2),MIN3(2), 
$ SEC3(2),HEAD( 2) 
TIMES5=DFLOAT ( DAY3( 2) )*¥SECD+DFLOAT( HOUR3( 2) )*SECH 
$ + DFLOAT(MINS3( 2) )¥60.D0+DBLE(SEC3(2)) 


READS ONE POSITION OF THE SHIP 
20 READ(50,ENDz100) TT , X ,Ү,2, DLAT , DLOG , DHT , NSAT 
COMPUTES THE DATE 


CALL SUB(TT ,DAY2 ;HOUR2 ;MIN2,SEC2) 
DAY2zDAY2*STDAY 


CONVERTS LATITUDE TO DEGREES MINUTES AND SECONDS 
CALL DMS(DLAT,LADEG,LAMIN,LASEG ) 
CONVERTS LONGITUDE TO DEGREES MINUTES AND SECONDS 
IN A RANGE O TO 180 
IF (DLOG.GT.180.D0) DLOG-360.D0-DLOG 
CALL DMS(DLOG,LODEG,LOMIN ,LOSEG ) 


WRITE( 6,24) LADEG,LAMIN,LASEG,LODEG ,LOMIN,LOSEG 
24 FORMAT(' LAT, LON (DEG MIN SEC, DEG MIN SEG) ',2(214,1X,F7.4),7) 


TIME3=DFLOAT( DAY 2 )*SECD+DFLOAT( HOUR2 )*SECH 
$ + DFLOAT( MIN2 )*60.DO*DBLE( SEC2) 


POSITION IN TIME BETWEEN THE TWO SETS OF PITCH AND ROLL DATA 
AND THE TWO SETS OF COURSE DATA 
25 CONTINUE 


POSITION TO EARLY FOR PITCH AND ROLL DATA 
IF(TIMES.GE.TIME1) GO TO 40 
WRITE(6,30) MONTH2 ,DAY2 ,YEAR2 »HOURZ »MIN2 »SEC2 
30 FORMAT(’ PSN 9 ",12,1Х,12,1Х,14,1Х,12,1Х,12,1Х,Ғ4.1, 
$ ' REJECTED. NO PITCH AND ROLL DATA' ) 
GOTO 20 


POSITION TO EARLY FOR PITCH AND ROLL DATA 
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&0 CONTINUE 
IF(TIMES.GE.TIME4) GO TO 60 
ARITE(6,50) MONTH2,DAY2,YEAR2 ,HOUR2 , MIN2 , SEC2 
50 FORMAT(' PSN 9 ',I2,1X,I2,1X,15,1X,12,1X, I2, 1X; F4.1, 
$ ' REJECTED. NO COURSE DATA') 
GOTO 20 


POSITION LATER THAN THE SECOND SET OF PITCH AND ROLL DATA 


60 CONTINUE 
IF(TIMES.LE.TIME2) GOTO 70 


MOVES THE PITCH AND ROLL DATA TO THE FIRST SET 


MONTH1( 1 )ΞΜΟΝΤΗΙ( 2) 
DAY1(1)=DAY1(2) 
YEARI(1)=YEARL(2) 
HOUR1(1 )=HOURL( 2) 
MIN1(1)=MIN1(2) 
SEC1( 2 )=SEC1(2) 
PITCH(1)=PITCH( 2) 
ROLL(1)=ROLL(2) 
TIMEL=TIME2 


READS A SECOND SET OF PITCH AND ROLL DATA 


READ(51,*,END=100) MONTH1(2),DAY1(2),YEAR1( 2 ),HOUR1(2),MIN1(2), 
$ SEC1(2),PITCH(2),ROLL(2) 

TIME 2=DFLOAT( DAY1( 2) )*SECD4+DFLOAT( HOURL( 2 ) )*SECH 
$ * DFLOAT(MIN1(2))*60.00*DBLE(SEC1(2)) 

GOTO 25 


POSITION LATER THAN THE SECOND SET OF COURSE DATA 


70 CONTINUE 
IF(TIMES.LE.TIMES5) GOTO 80 


MOVES THE COURSE DATA TO THE FIRST SET 


MONTH3(1)=MONTHS( 2 ) 
DAY3(1)=DAY3(2) 
YEARS(LISYEARSZ(2) 
HOURS(1)=HOURS( 2) 
MIN3( 1 J=MIN3(2) 
SEC3(2)=SEC3(2) 
HEAD( 1)-HEADt(2) 
TIMEG=TIMES 


READS A SECOND SET OF PITCH AND ROLL DATA 
READ(52,END=100) MONTH3(2),DAY3(2),YEAR3(2),HOUR3(2),MIN3(2), 
$ SEC3(2),HEAD( 2) 
TIMES=DFLOAT( DAY3(2 ) )¥SECD+DFLOAT( HOURS( 2) )¥SECH 
$ + DFLOAT(MIN3( 2 ))*60.D04DBLE(SEC3(2)) 
GOTO 25 
POSITION WITHIN THE TWO SETS IN TIME 
80 ВАТЕ=(ТІМЕЗ-ТІМЕ1 )/(ТІМЕ2-ТІМЕ1) 
PIT-PITCH(1)*( PITCH( 2 )-PITCH( 1) )*RATE 
ROL-ROLL( 1)*«( ROLL( 2 )-ROLL( 1) )*RATE 
RATES(TIMES-TIMEG )/( TIMES-TIMEA) 
COURSE :HEAD( 1)*( HEAD( 2 ) -HEAD( 1 ) )«RATE 
COMPUTS THE CORRECTIONS AND THE POSITIONS OF THE TRANSDUCER 
CALL PTH(OFFX,OFFY ,OFFZ,COURSE ; PIT ;,ROL,DX,DY ,DZ) 


COMPUTES THE UTM COORDINATES 
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CALL GPUTM( LADEG,LAMIN,LASEG,LODEG,LOMIN,LOSEG,YPOS,XPOS ) 


XPOS=XPOS +DX 
YPOS=YPOS+DY 


CONVERTS BACK TO GP 
CALL UTMGP(YPOS,XPOS ,LADEG,LAMIN,LASEG,LODEG,LOMIN,LOSEG } 
OUTPUTS THE RESULT 


WRITE( 7} MONTH2,DAYZ,YEAR2,HOUR2 ;MIN2,SEC2, 
8 LADEG,LAMIN,LASEG, LODEG, LOMIN, LOSEG »NSAT 


WARITE(6,90) MONTH2,DAY2,YEAR2,HOUR2 ,MIN2,SEC2, 


ғ LADEG,LAMIN>,LASEG,LODEG,LOMIN ,LOSEG »NSAT 


РОКМАТ(2(1Х,12),1Х,14,2(1Х,12),1Х,4.1,2(1Х,13,1Х,12,1Х,Е7.4),13) 


READS ONE MORE POSITION 
GOTO 20 


END OF PROGRAM 


100 CONTINUE 


STOP 
END 


SUBROUTINE PTH(XOFF,YOFF,ZOFF,COURSE PITCH;ROLL DX, DY ;D2Z) 


THIS SUBROUTINE COMPUTES THE CORRECTIONS TO THE 
COORDINATES DUE TO THE OFSET OF ANTENA IN RELATION 
TO THE TRANSDUCER 


PI-ARCOS(-1.0) 
ROLL=ROLL¥PI/180. 
РТТСН-РТТСНЖРТ/ 180. 
COURSE =COURSE*PI/180. 


CHANGE SIGNAL OF ORIGINAL DATA TO MATCH THE RIGHT HAND 
SYSTEM CONVECTION 


ROLL--ROLL 


COSPHI- COS(ROLL) 
SINPHI- SIN( ROLL) 
COSK= COS( COURSE } 
SINK= SIN( COURSE } 
COSW= COS(PITCH } 
ЅІМИ= SIN(PITCH) 


DX = XOFF * COSPHI αχ COSK + 


$ YOFF * ( COSM * SINK * SINW ж ЅІМРНІ κ COSK ) + 
$ ZOFF * { SINW * SINK - COSW ж SINPHI * COSK ) 
DY = XOFF * COSPHI * SINK * (-1.) + 
$ YOFF * ( COSM * COSK - SINM * SINPHI * SINK ) + 
$ ZOFF * ( SINM * COSK + COSW κ SINPHI * SINK ) 
02 = XOFF * SINPHI - 
$ YOFF * ( SINM * COSPHI ) + 
$ ZOFF * ( COSHM * COSPHI ) 
RETURN 
END 


SUBROUTINE GPUTM( LADEG,LAMIN,LASEG,LODEG,LOMIN,LOSEG»NORTH>EAST } 


DOUBLE PRECISION A,R»N>,AP,BP,CP,DP,EP»,S,R1,ESQ,ESQP ,RM,RP »KO 


100 


0000 


OQUSLE PRECISION R2,R3,R6,R5,P,P2,P3,P4,P5,P6,A6,B5,SINSEC 
DOUBLE PRECISION PHI,OLAM,NORTH,EAST,8,PHIMIN,PI,LON 
REAL LASEG,LOSEG 


THIS SUBROUTINE COMPUTS THE UTM COORDINATES OF GP IN WGS 72 
IN ZONE 10 CENTRAL MERIOIAN 123 00 00 W 


CM=-123.000 

PHI=OFLOAT( LADEG )+O0FLOAT( LAMIN)/60.00+0BLE( LASEG )/3600.00 
LON=OFLOAT( LOOEG )+OFLOAT( LOMIN )/60.00+0BLE( LOSEG )/3600.D0 
LON=-LON 

DLAM=( LON-CM }*3600.000 


А=6578155.000 
R=298. 2800 


KO=0.9996D0 
8 = A¥(R-1.00)/R 


N = (A-B)/C A+B) 


АР с Аж((1.00-М)%5.00/4.00%(Мжж2-МжҰ3 )%81.00/60.00ж(Мжж0-Мжж5)) 

ӨР с 3.00/2.00жАж( (М-Мжж2)%7.00/8.00ж МжҰ3-Мжжд4 )%55.00/64.00жМ4жж5) 
CP - 15.00/16.00*A*(N»**2-N**3423.00/4.00*( Nxx& -Nxx5 ) ) 

DP = 35.D0/48.00*A*( Nx«3 -N»x4*11.D0/16 . DOXNXX5 ) 

EP 3 315.D0/512.D0*A*( Nx*x4-N»*5 ) 


PHIMIN - PHI*60.00*2.90888208666D-4 

PI-OARCOS( 1.00) 

PHI-PHI/180.00*xPI 

S * APXPHIMIN-BP*DSIN( 2. DOxPHI )*CP*DSIN( 6 .DOXPHI) 
5 - DP*DSIN(L6.00xXPHI )*EPX*OSIN(8.00XPHI ) 

R1 = KO*S 


SINSEC = (1.00/3600.00)/180.00xPI 
SINSEC=OSIN(SINSEC } 


ESQ = (A%ww2-B8x%2 )/Aw%2 
ESQP = ESQ/(1.00-ESQ) 
RM = AX¥(1.D0-ESQ)/(DSQRT(1.D0-ESQXOSIN( PHI )w%%2 ) ) 33 


RP = RM¥(1.00+ESQP*OCOS( PHI )**2 ) 
m2 - RPXDSIN( PHI )*OCOS( PHI )¥SINSEC **2/2.00*KOx1.008 
RS = SINSEC##G*RPXOSIN( PHI )*DCOS( PHI )**3/24¢.D0%(5.D0-DTAN( PHI )жє2 


$ + 9, DO*ESQP*DCOS( PHI )€2+4 .00*ESQP¥ESQPxDCOS( PHI }**4 )*KO*1.016 
RG = RP*DCOS( PHI J¥SINSEC*KO*1 .DS 

= SINSEC¥*¥3*RPXOCOS( PHI )**3/6.00%(1.00-DTAN( PHI )**2 
$ * ESQPXDCOS( PHI )**2)xK0x*1.D12 


P = .000100*0LAM 


P2 = Pxx2 

PZ = PxxZ 

PG = PxxG 

P5 = Pxx5 

рө с ржжб 

A6 3 P6XSINSECX*6XRPXDSIN( PHI JX*0COS(PHI )x5/720.00 


$ * (61.00-58.DOXOTAN( PHI )**2«DTAN( PHI )»x& 
$ + 270.D0*ESQP*DCOS( PHI )**2-330 . DOXESQPXDSIN( PHI )*e*2 )xKO*1. 026 
B5 = PS*xSINSECX*5xRPX*DCOSt( PHI ):**5/120.00*(5.00-18.00xX0TAN( PHI )*x2 
$ * OTAN(PHI )33*€6* 16. DOXESQPX*0COS( PHZ ) 2 
$ - 58.00x*ESQPX*OSIN(PHI )**2)xK0x*x1.D20 
NORTH-R1*R2XxP24R3*P6*A6 
EAST-z( RGXP«R5xP3«4B5)*500000.00 
RETURN 
END 


SUBROUTINE UTMGP( NORTH ,EAST »LADEG,LAMIN,LASEG,LOOEG ,LOMIN,LOSEG) 
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DOUBLE PRECISION A,R,N>B, BP,CP,D0P,EP,S,R1,ESQ,ESQP,RM,RP ,KO 
DOUBLE PRECISION R7,R8,R9,E5,Q,Q2,Q3,Q6,Q5,Q6,D6, AP,SINSEC 
DOUBLE PRECISION R10,DLAM>,NORTH,EAST,RPHI,DPHI ,PHIMIN 
DOUBLE PRECISION EPRIME,DELTA,D6NUM,CM,PHI,LON,PI 

REAL LASEG,LOSEG 


THIS SUBROUTINE COMPUTS THE UTM COORDINATES OF GP IN WGS 72 
IN ZONE 10 CENTRAL MERIDIAN 123 00 00 W 


G OoOo 


CM--123.0D0 


O 


A=6378135.000 
R=298.26D0 


KO=0.9996D0 


PI=DARCOS(-1.00) 
B = A¥(R-1.00)/R 


М = (А-В)/(А+В) 
A*(C1.D0-N)*5.D0/4.DO*( N3e2 -N9ee3 ) t31.D0/64 . DOx ( Мзееф-Мвеб )) 
$.D0/2.DO0*A*( (N-Nx*2 ) *7.D0/8. DO*( NX9X3 -N9eeá ) +55.00/64.ООж Веб ) 
15.00/16.О0ждж( Ме? -Мжж3 +5.00/4.00#( Мжзф-Мекб )) 
$35.D0/48.D0*A*X( Nx*Z -N»x4*11.00/16.DO0*N»**5) 
315.D0/512.D0*A*( N»9& -N»*x5 ) 


Ω 
79 
пин ии 


FIRST APROXIMATION OF PHI 


ооо 


PHI-NORTH/30.8D0/3600.D0 


COMPUT TRUE MERIDIONAL DIST AND APROXIMATE PHI 


ооо 


DO 100 I-1,8 
PHIMIN-PHI*60.D0*2.90888208666D-4 
RPHI=PHI/180 .DO*PI 
5 с: APXPHIMIN-BP*DSIN(2.DO*RPHI ) *CP*XDSIN(4.DOX*RPHI) 
$ -DPX*DSIN(6.DO*RPHI ) *EP*DSIN(8.DO*XRPHI) 
Rl = KOX¥S 
DELTA=NORTH-R1 
PHI=(DELTA/30.8/3600. )+PHI 
100 CONTINUE 
RPHI-PHI/180.DO0*PI 


SINSEC 
SINSEC 


(1.D0/3600D0)/180.DO0*PI 
DSIN( SINSEC ) 


ESQ = (Aww2-B%%2 )/ Aw% 2 
ESQP - ESQ/(1.00-ESQ) 
RM * A*(1.D0 -ESQ)/( DSQRT( 1.00 -ESQXDSIN(RPHI )92 ) )x3 


RP * RMX*(1.DO*ESQP*DCOS(RPHI P962) 

R7*UDTANCRPHI )/( 2, DO*RP3s2x*SINSEC ) )*( 1. DO *ESQP*DCOS( RPHI )»2) 
$ x1.D12/KO*«2 

R8z( DTAN(LRPHI )/( 26 . DO«RPXQG*SINSEC ) )«( 5. DO* Z. DO*DTAN(RPHI )жж2 
$ *6.DOXESQP*DCOS( RPHI )**2-6.DO0X*ESQPX*DSIN(RPHI )92-3.DO*ESQPsex2 
3 *XDCOS(RPHI 96€4-9,DOX*ESQP329C2*DCOS ( RPHI )*2*DSIN(RPHI )*x2 ) 
$ #1.024/К03е4 

R9z1.D0/DCOS( RPHI )/( RP*SINSEC )*1.D6/KO 

R1071.D0/DCOS( RPHI )/( 6. DO*XRP»X*S3S*SINSEC )*( 1.D0*2.D0*DTAN( RPHI )*x2 
$ *ESQPXDCOS(RPHI )**2 )*1.D18/KO0s*3 


EPRIME-EAST-500000.D0 
9 с .00000100жЕРБІМЕ 


Q2 = Q**2 
Q3 = Q**3 
9% = Q**4 
Q5 = Qxx5 
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/* 


//GO. 
7/60. 
//60. 
//GO. 
//GO. 
//GO. 


/” 


Qó = Q%%6 
DeNUM-Q6*DT ANtRPHI ) 
D6z(D6NUM)/(720.DOXRPX*6*SINSEC )*(61.00«90.D00X0TAN( RPHI )**2 
$ +45.D0*DTAN( RPHI )*¥4+107.DO*ESQP*OCOS( RPHI )**2 
$ -162 .DOXESQP*DSIN( RPHI )**2-45. DOXESQPXDTAN( RPHI )**2 
5 XOSIN(RPHI ) 2 )#1.056/КОЖж»6 
E5s(Q5*1.00/0COS(t RPHI ))/(120.00X*RP»*5x*SINSEC )*(5.00*28.00 
$ XDTAN(RPHI )*2*24.00xDT ANCRPHI )*4*6. 00*ESQPXDCOSt RPHI )*x*2 
$ *8.DOXESQPXDSIN(RPHI 2 )*1.030/K036&*5 


DPHI=( -R7*Q2Z +R8*¥Q4-D6 )/3600.D0 
OLAM=(R9*Q-R10*Q3+E5 }/3600.00 
PHI=PHI+DPHI 

LON=CM+DLAM 

CALL OMS( PHI ,LADEG,LAMIN,LASEG ) 
CALL OMS(LON,LOOEG,LOMIN,LOSEG) 
RETURN 


ENO 


SUBROUTINE OMS(DEC,LOEG,MIN,SEC ) 


DOUBLE PRECISION DEC,XNUM,XMIN 


XNUM=DABS( DEC ) 
LOEG=DINT ( XNUM ) 
XMIN=( XNUM-DFLOAT(LDEG ) )*60.00 
MIN-DINT t XMIN) 
XNUM- t XMIN-OFLOATtMIN))*60.00 
SEC-SNGL ( XNUM ) 
IF (LDEG.GE.360) LDEG-LDEG-360 
RETURN 


END 


SUBROUTINE SUB(TT,DAY,HOUR »,MIN,SEC ) 
INTEGER DAY, HOUR, MIN 
DOUBLE PRECISION TIME,TT 


REAL SEC 


ТТМЕ <ТТ 
DINT(TIME / (24.00x3600.00)) 


DAY - 
TIME 
HOUR 
TIME 
MIN = 
TIME 
SEC = 
RETUR 
END 


ΕΤΟΕΕΟΟΙ 
ЕТ07Ғ001 
ЕТ5ОЕО01 
FT50F001 
FT51F001 
FT52F001 


TIME - DFLOAT(DAY) w 24.D0 ж 3600.00 
DINT(TIME / 3600.00) 
TIME - DFLOAT(HOUR) ж 5600.00 


DINT(TIME / 60.00) 


TIME - DFLOAT(MIN) ж 60.00 


SNGL ( TIME ) 


N 


SYSOUT =* 

DSN=MSS .S$0812.GPS.TRANSDUC. POS ,DISP=SHR 
DSN=MSS ..S0812.CVDATA,DISP=SHR 

DSN=MSS .S0812 .GPS. ANTENNA .POS ,DISP=SHR 
DSNzMSS.S0812.PITROLL ; OISP-SHR 

OSN=MSS .S0812.S0AS .COURSE ,OISP=SHR 
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APPENDIX G 
PROGRAM COMPARE POSITION. SOURCE LISTING 


// JOBCOMP A JOB (0812,9999), ‘EZEQUIEL' ,CLASS=C 
//ЖМАІМ ORG=NPGVM1.0812P ,LINES=( 99) ,CARDS=( 99 } 
//*FORMAT PR,DDNAME -GO.FTO6FO001,; 

//*DESTzNPGBACH 


// 


EXEC FORTVCLG 


//FORT.SYSIN DD * 


GO GO O G O G O O O O O O O O O O O GO O OG O G OO GO G GO O оО 


ооо 


ооо 


ооо 


РКОСКАМ СОМРАКЕ РОЗТТТОМ 


AUTHOR: AUGUSTO EZEQUIEL 
DATE: MARCH 27, 1987 


DESCRIPTION: 


THIS PROGRAMS TAKES EACH GPS POSITION , FINDS THE 
CORRESPONDING POSITION USING THE MR FALCON, COMPUTES THE 
4» Y AND RADIAL DIFFERENCES AND PRINTS THE DATE TIME TAG 
GEOGRAPHIC POSITIONS AND DIFFERENCES. THIS PROGRAM READS 
THE PERIODS TO COMPARE THE DATA FROM A SEPARATE FILE. 


THE PROGRAM RUNS IN MVS 
I/O SPECIFICATIONS: SEE END OF THIS JOB. 


THIS PROGRAM IS RESTRICT AS IS TO THE SEAFLOOR BENCHMARK 
EXPERIMENT 


THE PROGRAM WILL RUN WITH ANY AMOUNT OF DATA ONLY LIMITED TO DISK 
SPACE 


ANY BLANK LINES WILL TERMINATE THE PROGRAM IN ERROR 


DOUBLE PRECISION XPOS(3),YPOS(3),TIME(5),SECD,SECH,RATE 
DOUBLE PRECISION DIFX,DIFY ,DIST,XNEW,YNEW 

REAL SEC(5),LASEG(3),LOSEG(3) 

INTEGER MONTH(5),DAY(5)  YEARC5 ) ;HOUR( 5 ) ;MIN( 5 ) ;NSAT 

INTEGER LADEG( 3) ;LAMIN( 3 ) ; LODEG( 3 ) | LOMIN(C3);LINE ;LOOP ; INDEX 


INITIALIZATION OF CONSTANTS 
SECH=3600.D0 
SECD=24¢.00*SECH 
OIFX=0.0 
DIFY=0.0 
OIST=0.0 
INPUTS THE NUMBER OF PERIODS TO COMPARE THE DATA 
КЕА0(25,%,ЕМ02100) LOOP 
LOOP OF THE PROGRAM TO THE NUMBER OF TIMES DESIRED 
DO 1000 INDEX = 1, LOOP 
LINE =60 
REWIND( 30 ) 
REWIND( 50 } 


READ(25,%,END=100) MONTH(4),DAY(4),YEAR(4) »HOUR(4) »MIN(G) ,SEC(G) 
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$ ,MONTH(5),0AY(L5) , YEARC 5) ,HOUR( 5 ), MIN(C 5), SEC( 5) 


C 
C STARTING TIME 
C 
TIME(4)-zDFLOAT(DAY( 4) )*SECD*DFLOAT ( HOUR( 4 ) JXSECH 
$ + DFLOAT(MIN(4))0*6D.DD*DBLE(SEC(4)) 
C 
C ENDING TIME 
C 
TIME( 5)zDFLOAT(DAY( 5) )J*SECD*DFLOAT ( HOUR(S ) )XSECH 
$ + DFLOAT(MIN( 5) )*60.DD+DBLE(SEC(5)) 
C 
C READS TWO MR FALCON POSITIONS 
C 
10 READ(30,END=1DD) MONTH(1),DAY(1),YEAR(1),HOUR(1),MIN(1),SEC(1), 
$ LADEG(1),LAMIN(1),LASEG(1),LODEG(1),LOMIN( 1),LOSEG(1) 
TIME (1 J=DFLOAT( DAY( 1) )*SECD+DFLOAT( HOUR( 1) )*SECH 
$ + DFLOAT(MIN( 1) )¥60.D0+DBLE(SEC(1)) 
C 
C POSITIONS THE FIRST POSITION INSIDE THE PERIOD. THE REST, OF THE 
C PROGRAM WILL ADJUST 8Y IT SELF 
C 
IF( TIME(1).LT.TIME(4) ) GO TO 10 
C 
C 
READ( 3D,ENDz100) MONTH(2),DAY(2) , YEARC2 ) ;HOUR(C2 ) ,MIN( 2), SEC( 2), 
$ LADEGt( 2 ) ; LAMIN( 2 ) ; LASEG( 2 ) , LODEG( 2 ) ; LOMIN(C2 ) ,LOSEG( 2) 
TIME (2 J=DFLOAT( DAY (2 ) J* SECD*DFLOAT ( HOUR( 2 ) J*SECH 
5 + DFLOAT(MIN( 2) )*60.D0+DBLE(SEC(2)) 
C 
C COMPUTES THE UTM COORDINATES 
C 
CALL GPUTM(LADEG(1),LAMIN(1),LASEG(1),LODEG(1),LOMIN(1),LOSEG(1), 
$ YPOS(1),XPOS(1)) 
CALL GPUTM( LADEG( 2) ; LAMIN( 2 ) ; LASEG( 2 ); LODEG( 2) ; LOMIN(C2 );LOSEG( 2) 
$ YPOS(2),XPOS(2)) 
C 
C READS GPS POSITION OF THE SHIP 
C 
2D READ(5D,END=1DD) MONTH(3),DAY(3),YEAR(3),HOUR(3),MIN(3),SEC(3); 
$ LADEG(3),LAMIN(3),LASEG(3),LODEG(3),LOMIN( 3 ),LOSEG(3),NSAT 
CALL GPUTM( LADEG(3),LAMIN(3),LASEG(3),LODEG(3),LOMIN(3),LOSEG(3), 
$ YPOS(3),XPOS(3)) 
С 
TIME(3)=DFLOAT(DAY(3))%SECD+DFLOAT(HOUR( 3) )XSECH 
$ + DFLOAT(MIN(3))%60.DD+DBLEI(SEC(3)) 
C 
C POSITION OUTSIDE OF THE PERIOD 
C PROGRAM MILL PICK ANOTHER PERIOD 
C 
ІР( TIME(3).GT.TIME(5) ) GO TO 1DDD 
C 
C POSITION IN TIME BETWEEN THE TWO MR FALCON POSITIONS 
C 
25 CONTINUE 
C 
C POSITION TO EARLY FOR MR FALCON DATA 
C 


IF(TIME(3).GE.TIME(1)0) GO TO 40 
Cx WRITE(6,30) MONTH(3),DAY(3),YEAR(3),HOUR(3),MIN(3),SEC(3) 
Сж 30 FORMAT(' GPS PSN dð ',12,1X,12,1X,I16,1X,I2,1X, 12, 1X,F4.1, 
Сж $ ' REJECTED. NO MR FALCON DATA') 


GOTO 2D 
C 
C POSITION LATER THAN THE SECOND MR FALCON POSITION 
C 
40 CONTINUE 
IF(TIME(3).LE.TIME(2)) GOTO 70 
С 
С MOVES THE SECOND MR FALCON POSITION TO THE FIRST SET 
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MONTH( 1 )z MONTH(2 ) 
DAY(C1)20AY(2) 
YEARC1)ZYEAR(C2) 
HOUR( 1 )ZHOUR( 2 ) 
МІМ(13-МІМ(2) 
SEC( 2 J=SEC( 2) 
TIME(1)-TIME(2) 
LADEG(1)-LADEG(2) 
LAMIN(1)-LAMIN(2) 
LASEG(1 J=LASEG( 2) 
YPOS(1)=YPOS( 2) 
LODEG( 1 J=LODEG( 2) 
LOMIN( 1 J=LOMIN( 2 ) 
LOSEG(1 J=LOSEG( 2) 
XPOS(1 }=XPOS( 2 ) 


READS A SECOND SET OF MR FALCON DATA 


READ(30,END=100) MONTH(2),DAY(2),YEAR( 2),HOUR(2),MIN(2),SEC(2), 
$ LADEG( 2},LAMIN( 2),LASEG(2),LODEG(2),LOMIN( 2 ),LOSEG( 2} 
TIME(2)zDFLOAT(DAY( 2) )*SECD«*DFLOAT ( HOUR( 2 ) )JXSECH 
$ * DFLOAT(MIN(2))*X60.D0*DBLE(SEC(2)) 
CALL GPUTM( LADEG( 2) ; LAMINC2) ; LASEG( 2 ); LODEG( 2) ; LOMINC2) , LOSEG(2); 
$ YPOS(2) ,XPOSU2)) 
GOTO 25 


GPS POSITION WITHIN THE TWO MR FALCON POSITIONS 
70 RATE-(TIME(3)-TIME(1))/(TIMEC2)-TIME(C1)) 


XNEW=XPOS( 1 ) +( ХРОЅ( 2 )-XPOS(1) )*RATE 
ҮМЕЙ=ҮРОЅ(1 )+(ҮРОЅ( 2 )-ҮРОЅ(1 ) )жВАТЕ 


ОІРХ=ХРОЅ(3 )-ХМЕМ 
DIFY=YPOS( 3 )-YNEW 
DIST=DSQRT( DIFX**¥2+DIFY *¥2 ) 


OUTPUTS THE RESULT 


IF (LINE.LT.60) GO TO 85 
WRITE( 6,80 ) 
80 FORMAT(1H1,/,17X%,'DIFFERENCES BETWEEN GPS AND MR FALCON POSITIONS' 
»//,29X%,'AVOIDING SV11' , 
//э9Х» 'DATE' , 6X, 'TIME',09X, 'LATITUDE' , 7X, 'LONGITUDE' ; 
12X,'(METERS)',9X, '$8' ,/,2X, 'M' 2X, 'D' 1X, 'YEAR' ; 1X, 'HH' ; 
1X,'MM',1X,'SS.S',2X, 'DD', 1X, 'MM', 1X, ''SS.SSSS' , 3X, 'DDD' ; 1X, 
'"MM',1X,'SS.SSSS', 5X, DX', 5X, DY' ,B5X, DIST. ', 1X,'SAT. ' ,7) 
LINE=0 
85 LINE=LINE+1 
WRITE(6,90) MONTH(3),DAY(3),YEAR(3 ),HOUR(3),MIN(3),SEC(3)>5 
$ LADEG(3 ),LAMIN(3 ),LASEG(3 }),LODEG(3 ),LOMIN(3 ),LOSEG(3); 
$ DIFX,DIFY »DIST »NSAT 
90 FORMAT(2(1X,12),1X,I4,1X, 12, 1X, 12; 1X; F6.1, 1X, I3, 1X, I2, 1X5 F7.6, 
$ ' N',1X,I3,1X,12,1X,F7.4,' W',1X,F6.25,1X,;F6.2,1X, F7.25,1X5 12) 


ww «4» «4» 4» 


READS ONE MORE POSITION 
GOTO 20 


END OF PROGRAM 
1000 CONTINUE 
100 CONTINUE 


STOP 
END 


SUBROUTINE GPUTM(LADEG,LAMIN,LASEG,LODEG,LOMIN, LOSEG,NORTH>EAST ) 
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DOUBLE PRECISION A,R,N,AP,BP,CP,DP,EP,S,RI,ESQ,ESQP , RM, RP,KO 
DOUBLE PRECISION R2,R3,;,R6,R5,P,;P2,P2,P4,P5,P6,A6,B5,SINSEC 
DOUBLE PRECISION PHI,DLAM,NORTH,EAST ,B,PHIMIN,PI,LON 

INTEGER LADEG,LAMIN,LODEG,LOMIN 

REAL LASEG,LOSEG 


THIS SUBROUTINE COMPUTS THE UTM COORDINATES OF GP IN WGS 72 
IN ZONE 10 CENTRAL MERIDIAN 123 00 00 W 


CM=-123.0D0 

PHI=DFLOAT( LADEG )+DFLOAT( LAMIN )/60.00+DBLE( LASEG )/3600.D0 
LON=DFLOAT(LODEG )+DFLOAT( LOMIN)/60.D0+DBLE( LOSEG )/3600.D0 
LON=-LON 

DLAM=( LON-CM )*3600.0D0 


A=6378135.0D0 
R=298.26D0 


ΚΟΞΟ. 999600 
В = A¥(R-1.D0)/R 


N = (А-В)/(А+В) 


АР = AX((1.D0-N)*5.D0/4.DO*( NX*2-Nx*3)«81.D0/64 . DOX( N3XPEG -NDEXEB ) ) 

BP = 3.D0/2.D0*A*( ( N-Nxx2)*7.D0/8.DO0x( NX -NXPI$ 1 +55 .D0/64. DOXN**5 ) 
СР = 15.D0/16.D0*A*( Nx 2-N»9 3 43. D0/6 . DOX( N»PG -NXEXE5 ) ) 

ОР = 35.D0/658.D0X*A*( Nxx*3-Nx*4*11.D0/16.D0xN395) 

EP = 315.D0/512.DO0*A*( NxsGá -N3ex 5 ) 


РНІМІМ г РНІ%60.00%2.908882086660-% 

PI=DARCOS(-1.D0) 

PHI-PHI/180.DO0*PI 

S = APX*PHIMIN-BP*DSIN(2.DO*PHI )+CP¥DSIN(4.DOXPHTI ) 
$ - DPXDSIN(6.DOXPHI)*EPXDSIN(8.DOXPHI) 

Rl » KO*S 


SINSEC - (1.D0/3600.D0)/180.D0x*PI 
SINSEC-DSIN( SINSEC ) 


ESQ - (Axx2-B»e*x2 )/Axxe2 
ESQP = ESQZ(1.D0-ESQ ) 
RM = Aw(1.D0-ESQ )Z/( DSQRT( 1. DO-ESQXDSIN( PHI ):e2 ) )x>*3 


RP = RM¥(1.D0+ESQP*DCOS( PHI )*2 ) 
R2 - RPX*DSIN(PHI )*xDCOS( PHI JS*SINSEC*2/2.D0*KO*1.0D8 
R3 » SINSEC*4X*RPXDSIN(PHI )*DCOS( PHI )x*3/24.D0*(5.D0-DTAN(PHI ) 2 


5 * 9.DOX*ESQPXDCOS ( PHI ex2*«4 . DO*XESQPX*ESQP*DCOS ( PHI )3eà )xK0*1.D16 
RG = RP*DCOS(PHI )*SINSECXKO*1.D4 
R5 - SINSEC*3x*RPx*DCOS( PHI )»9*3/6.D0*( 1. DO-DTAN(PHI )xxe2 

$ * ESQPX*DCOS( PHI )**2 )*KOx*x1.D12 


P = .0001DOX*DLAM 


P2 * P*x2 
рз = Ржжз 
PG = PREG 
p5 = P%x5 
P6 - P3*6 
Аб * P6X*SINSEC3?o»*6*RPX*DSIN( PHI )*DCOS ( PHI )*3€5/720.D0 


$ κ (61.00-58.ΡΟΧΡΤΑΝί ΡΗΙ ) 692 +ОТАМ(РНТ ) 4 
$ + 270.DO0X*ESQPX*DCOS ( PHI )*2- 330. DO*ESQPX*DSIN(PHI )9*2 )xKO*1.D24 
В5 = PB5XSINSECX*5x*RP*DCOS( PHI )::€*5/120.D0*( 5. D0-18.DOXDTAN( PHI ) 2 
$ * DTAN(PHI )*6*16. DOXESQPXDCOS( PHI x2 
$ - 58.D0X*ESQPXDSIN(PHI )*2 )*KO*1.D20 
МОКТНЕВ 1 + К2жр2 + ЗжР4 + д6 
EAST=(R4*P+R5*P3+B5 )+500000. 00 
RETURN 
END 
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ЕМО 
/* 
//GO.FTO6FOOL DD SYSOUT=% 
/760.ҒТ25Ғ001 00 DSNzMSS.S0812.PERIODS.COMPARE , DISP=SHR 
7/60.ЕТ20Ғ001 00 DSN=MSS.S0812.FALCON.TRANSDUC .POS,DISP=SHR 
7/60.ЕТ50Ғ001 00 DSNzMSS.S0812.GPS.TRANSDUC .POS ,DISP-SHR 
/7 
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APPENDIX H 
PROGRAM COMPARE PLOT. SOURCE LISTING 


//EZEQUIEL JOB (0812,9999), ‘EZEQUIEL' ,CLASS=C 


// 


ЕХЕС FRTVCLGP 


//FORT.SYSIN DD * 


G GO O GO GO O G G G O G oO G 


aan 


ооо ооо ооо ооо ооо 


ооо 


PROGRAM COMPARE PLOT 

RUNS IN FORTRAN VS 

AUTHOR: AUGUSTO EZEQUIEL 

DATE : 31 MARCH 1987 

THIS PROGRAM MAKES A PLOT OF THE POSITIONS OF THE SHIP 
USING GPS AND MR FALCON POSITIONS OF THE TRANSDUCER 


DOUBLE PRECISION XPOS >YPOS,XLEFT >YLEFT , ТТМЕ( 4 ) „ЗЕСО „ЗЕСН 
REAL SEC(3),XPLT,YPLT,BLXH,BLYH,SCALE ,VALUE ,LASEG, LOSEG 
REAL YORG,XORG 
INTEGER MONTH(3),YEAR(3),DAY(3),HOUR(3),MIN(3) 
INTEGER IPEN,LADEG,LAMIN,LODEG,LOMIN 
INTEGER I0,LOOP,INDEX 
DIMENSIONS OF SHEET, LEFT CORNER AND SCALES 
BLXH=15. 
BLYH=20.5 
SCALE-1./5000. 
AORG-0.0 
YORG=-25.0 
INITIALIZATION OF CONSTANTS 
ТІМЕ(4):0.00 
SECH=3600.D0 
SECD=24¢.D0*SECH 
PLOTTER INITIALIZATION 
CALL PLOTS (0,0,0) 
READS THE AMOUNT OF PERIODS TO COMPARE 
READ(30,%) LOOP 
DO 100 INDEX=1,LOOP 
READS THE PERIODS AND INFERIOR LEFT CORNER OF EACH AREA 


READ(30,*) MONTH(1),DAY(1),YEAR(1),HOUR(1),MIN(1),SEC(1) 
5 ,MONTH(2),DAY(2),YEAR(2),HOUR(2),MIN(2),SEC(2), XLEFT,YLEFT 


STARTING TIME 


TIME( 1 )=DFLOAT(DAY( 1) )¥SECD+DFLOAT(HOUR( 1 ) )¥SECH 
$ + OFLOAT(MIN( 1) )¥60.DO+O0BLE(SEC(1)) 


ENDING TIME 


TIME( 2 )=DFLOAT( DAY( 2) )¥SECD+DFLOAT( HOUR( 2 ) )ж5ЕСН 
$ + OFLOAT(MIN( 2) )*60.D0+DBLE(SEC(2)) 


109 


oan ооо 


ооо 


ооо 


ООО ооо ооо ооо ооо оо о ооо 


ооо 


POSITIONS THE PEN ІМ ΤΗΕ ORIGIN 


YORG=-YORG 
XORG225.0*(1-(-1)x*x*xINDEX)/2 
IFC INDEX.EQ.2) XORG=0 

CALL PLOT ( XORG,YORG,-3) 


GRIO 
CALL PLOT (0.1,0.1,3) 
CALL PLOT (0.1,BLYH-0.1,2) 
CALL PLOT (BLXH-0.1,BLYH-0.1,2) 
CALL PLOT (BLXH-0.1,0.1,2) 
SCALE 
CALL METER( SCALE } 


TITLE 


CALL SYMBOL(08.1,19.2,0.20,29HSEAFLOOR BENCHMARK EXPERIMENT,O. 


CALL SYMBOL(08.1,18.8,0.20,29H PHASE II 
CALL SYMBOL(08.1,18.4,0.20,29H 16 AUGUST 1986 
CALL SYMBOL(08.1,18.0,0.20,29HAVOIDING SV11. NO CONSTRAIN 


CALL SYMBOL(O1.0,03.375,0.15,3,0.,-1) 


CALL SYMBOL(1.20,03.3,0.15,29HMR FALCON TRANDUCER POSITIONS;0. 


CALL SYMBOL(01.0,03.075,0.15,11,0.,-1) 
CALL SYMBOL(1.20,03.0,0.15,29HGPS TRANOUCER POSITIONS 


SETS SYMBOL FOR MINI RANGER POSITIONS 
ISN=3 
REMINDS THE FILES 


REWINO( 31) 
REWIND( 32 ) 


LOOP FOR THE TYPES OF POSITION 
00 80 10:51,22 

SETS PEN UP 
IPEN=3 


LOOP WITHIN ONE TYPE OF POSITION 


10 CONTINUE 


READ THE DATA 


READ( I0,END=50 JMONTH( 3) ,DAY(3) ,YEAR(3),HOUR(3); 
$ MIN(3),SEC(3),LADEG,LAMIN,LASEG, LODEG,LOMIN,LOSEG 


COMPUTES THE UTM COORDINATES 
CALL GPUTM( LADEG,LAMIN,LASEG ,LODEG,LOMIN,LOSEG,YPOS »XPOS ) 
COMPUTES THE TIME IN SECS 


TIME (3 J=DFLOAT(ODAY(3 ) JX¥SECD+DFLOAT( HOUR(S ) )¥SECH 
$ + DFLOAT(MIN( 3 ) )*60.00+DBLE( SEC(3)) 
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229) 
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TESTS IF NITHIN THE PERIOD 


IF(TIME(3).GT.TIME(2)) GOTO 50 
IF(TIME(3).LT.TIME(1)) GOTO 10 


COMPUTES THE PLOTTER COORDINATES 


XPLTz(XPOS-XLEFT )*SCALE*100. 
YPLTz(YPOS-YLEFT )*SCALE*100. 


TESTE IF INSIDE AREA 


IF(XPLT.LT.O.) GOTO 30 
IF(YPLT.LT.O.) GOTO 30 
IF(XPLT.GT.BLXH) GOTO 30 
IF(YPLT.GT.BLYH) GOTO 30 


TESTS IF POSITIONS ARE AWAY MORE THEN 10 SECS IN TIME 


IF(TIME(3)-TIME(%).GT.10.00) IPEN=3 
TIME (4 )=TIME(3 } 


PLOTS THE POSITION 


CALL PLOT(XPLT ,YPLT ; IPEN) 
IPEN=2 


PLOTS THE TIME EVERY MINUTE 


IF(IFIX(SEC(3)).GT.0) GO TO 10 

VALUE zFLOAT ( HOUR( 3 ) )* 100 *F LOAT( MIN(3)) 

CALL NUMBER(IXPLT*0.15,YPLT*0.15,0.15, VALUE 0. ,-1) 
CALL SYMBOL(XPLT,YPLT,0.15,ISN>0O. »-1) 

GOTO 10 


PEN UP WHILE THE POSITIONS ARE OUT OF THE SHEET 


30 IPEN=3 
GOTO 10 


END OF PLOT 
50 CONTINUE 


ISN=11 
IPEN=3 
80 CONTINUE 
100 CONTINUE 
CALL PLOT(O. »0. »+999) 
STOP 
END 


SUBROUTINE METER( SCALE ) 


Х0:2.0 
YO=1.0 
CALL PLOT(X0,Y0,-3) 


DO 10 J=1,12 

XP= FLOAT(J)-1.0 

CALL PLOT(XP,0.0,3) 

CALL PLOT(XP,0.25,2) 
10 CONTINUE 


DO 20 J=1,9 
XP=FLOAT(J)*0.1 
CALL PLOT(XP,0.0,3) 
CALL PLOT(XP,0.2,2) 


lli 


оосо с) 


20 CONTINUE 


CALL PLOT(0.,0.0,3) 
CALL PLOT(11.,0.,2) 
CALL PLOT(0.,0.2,3) 
CALL PLOT(11.,.2,2) 


VALUE=1./SCALE/100. 

CALL NUMBER(-0.25,0.27,0.15,VALUE ,0.0,-1) 
CALL NUMBER(0.95,0.27,0.15,0.0,0.0,-1) 
VALUE=5.0/SCALE/100. 

CALL NUMBER(5.775,0.27,0.15, VALUE 0.0, -1) 
VALUE=10.0/SCALE/100. 

CALL NUMBER(10.7,0.27,0.15,VALUE ,0.0,-1) 
CALL SYMBOL(5.5,-0.25,.15, 6HMETERS,0.0,6) 


CALL SYMBOL(10.0,12.0,1.5,62,0.0,-1) 
CALL SYMBOL(10.1,12.4,0.2,85,0.0,-1) 
ΧΟΞ-ΧΟ 

ΥΟΞ-ΥΟ 

CALL PLOT(X0,Y0,-3) 


RETURN 
END 


SUBROUTINE GPUTM( LADEG,LAMIN,LASEG ,LODEG,LOMIN,LOSEG »NORTH,EAST ) 


DOUBLE PRECISION A>R,N,AP>BP>CP,DP,EP>S,R1,ESQ,ESQP ,RM,RP KO 
DOUBLE PRECISION R2,R3,R4,R5,P,P2,P23,PG,P5,P6,A6,B5,SINSEC 
DOUBLE PRECISION PHI,DLAM,NORTH,EAST »,B,PHIMIN,PI,LON 

REAL LASEG,LOSEG 


THIS SUBROUTINE COMPUTS THE UTM COORDINATES OF GP IN WGS 72 
IN ZONE 10 CENTRAL MERIDIAN 123 00 00 W 


CM=-123.0D0 

PHI=DFLOAT(LADEG )+DFLOAT(LAMIN)/60.D0+DBLE(LASEG )/3600.D0 
LON=DFLOAT( LODEG )+DFLOAT( LOMIN)/60.00+DBLE( LOSEG )/3600.D0 
LON- -LON 

DLAM={( LON-CM }¥3600.0D0 


A=6378135.000 
R=298.26D0 


KO=0.9996D0 


В = A¥(R-1.D0)/R 


М = (A-B)/(A*B) 

АР = А#(1(11.00-М)+5.00/4.00#( Мед -Мж3 )+81.00/64.00#( Кеф-МжзеБ )) 

BP = 3.D0/2.D0*A*(t (N-Nxx2)*7.D0/8.DO*( Nxx3-N3xx&4 ) *55.D0/66 , DOXNXPX5 ) 
СР = 15.D0/16.D0*A*( Nxx2-Nxx3«3.D0/6 . DOx( NxX&4 -N3x5 ) ) 

DP = 35.D0/68.D0*A*( Nxx3-N»x4 *11.D0/16.DOxN3x*5) 

EP 7 315.D0/512.D0*A*( Nx*&-N»ex5 ) 


PHIMIN - PHI*60.D0*2.90888208666D-4 

PI-DARCOS(-1.D0) 

PHI-PHI/180.DO0XPI 

5 - APXPHIMIN-BPXDSIN(2.DOxPHI )*CPX*DSIN(4.DOXxPHI) 
$ - DPXDSIN(6.DOXPHI )*EPXDSIN(8.DOXPHI) 

Rl = KOxS 


SINSEC г (1.00/3600.002/180.00жРІ 
SINSEC-DSIN(SINSEC ) 


ESQ - (Ax*2-Bs*2 )/Axx2 
ESQP - ESQ/(1.D0-ESQ) 


112 


RM 


Ax(1.D0-ESQ)/L( DSQRT(C1.DO-ESQXDSIN(PHI )x*2 ) )*3 


ë 
RP = RM*(1.D0+ESQP*DCOS(PHI ) 92) 
R2 = RP*DSIN(PHI)*DCOS(PHI)%SINSECx*2/2.D0%wKO*1.0D8 
АЗ = SINSECX*4xRPXDSIN( PHI )XDCOS( PHI )x*2/26 .D0x( 5. DO-DTAN(PHI )2 
$ * 9.DO0XESQPXDCOS( PHI 32*4.DOXESQPXESQGPXDCOS( PHI )**4 )xKOx1.D16 
R& s RPX*DCOS(PHI)XSINSECXKOx1.D4 
В5 = SINSECXx*zxRPXxDCOS(PHI 392/6.D0x( 1. D0-DTAN(PHI «2 
$ + ESQPXDCOS( PHI )%*2)*KO*¥1.D12 
5 
P = .0001D0XxDLAM 
P2 = Pxx2 
рз = PxxZ 
ра = ржжё 
P5 = ржж5 
P6 = Pxx6 
С 
AG = P6XSINSECXX*6XRPXDSIN( PHI )*DCOS( PHI )**5/720.00 
$ зж (61.00-58.DOXDTAN(PHI 12 *DTANCPHI 966 
$ + 270.00XESQPXDCOS(PHI 9«2-320.DOXESQPXDSIN(PHI x2 )*KOXx1.D26 
B5 = P5X*SINSECXx*5x*RPXDCOS(PHI )*5/120.D0x*(5.D0-18.D0XDTAN( PHI 1x2 
$ + DTAN( PHI )**4+14.00%ESQP*DCOS( PHI )**2 
$ = 58.DO*ESQP*DSIN( PHI )**2 )*KO¥1.020 
МОВТН-В1%В2жР24В3жР44А6 
ЕАбТ-( ВежраВ5жр3%85)%500000.00 
RETURN 
END 
/ж 


//GO . PLOTPARM 


DD ж 


&PLOT XMIN=0. »XMAX=300. »YMIN=0. »YMAX=51. ,SCALE=1. ,UNITS=2.540 &END 


//GO.FTO6FOO1 
//GO.FT30FO001 
//GO.FT31F001 
//GO.FT32F001 
/7 


DD SYSOUT-* 

DD DSN-MSS.S0812.PERIODS.COMPARE ;DISPzSHR 

DD DSN-MSS.S0812.FALCON. TRANSDUC . POS ;DISPzSHR 
DD DSNzMSS.S0812.GPS.TRANSDUC.POS , DISP -SHR 
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APPENDIX I 


CONTROL FILE FOR KALMN2 PROGRAM. AS DEFINED BY 
BROWN(1986) 


985590 700005. 1 6 10 5 1 1 ο. 0 
0 -2740.506 -2321 156 3772.0620 «διαδτή ο ο το ο RO 
6378.135 ο. 5 
1.01.0 1.0 1 .Е-ТА 2.1 τ ΕΞ ΤΕ 1: νι. 7’; 
1 0 1 0 0 0 0 1 0 0 
0.000005 0.0 
1 1 0 0 0 0 1 0 1 
ος ο. ου 
0.75Е-4 5н-4 ΤΕ. 1Ε-5 1.0 1.0.1.0 2.07 ΕΞ: 0: 
1.0 1.0 1.0 1.0 1.0 1.0 ЕРОЛ БЫ ЕТЕ 
-. 031 
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APPENDIX J 
CONTROL FILE FOR CVFICA РКОСКАМ. 


The first two lines may have a title up to 80 columns. The third line has the 
month dav and vear of the starting dav of the GPS week. The fourth line has the 
number of satellites to be used, and other two constants are kept at zero. The last line 
has the estimate position in WGS 72 coordinates. 


The next is the one used for this thesis: 


SEAFLOOR BENCH MARK PHASE II 
R/V POINT SUR 16 AUGUST 1986 
8 10 1986 
- 0 Q. 
-2739.58 -4340.77 3773.14 
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APPENDIX K 
CONTROL FILE FOR COMPARE POSITION AND PLOT PROGRAMS 


The first lines has the number of periods to compare the data. This line must be 
followed by the same number of lines as periods are to compare. Each line the date 
time tags from the start and of the periods has the left inferior corner of the area 
corresponding to each period. 


The next 1s the one used for this thesis: 


ll 

8 16 1986 19 140.0 8 16 1986 19 20 3.0 567320.0 4038400.0 
8 16 1986 19 20 0.0 8 16 1986 19 25 3.0 566920.0 4039000.0 
8 16 1986 19 25 0.0. 8 16 1986 19 30 3.0 566670.0 4039600.0 
16 1986 19 30 0.0 8 16 1986 19 35 3.0 566270.0 4040200.0 
10 1986 19 35 0.0 8 16 1986 19 40 3.0 566020.0 4040800.0 
16 1956 19 47 0.0 8 16 1986 19 53 3.0 564920.0 4040950.0 
16 1986 19 570.0 8 16 1986 20 05 3.0 564970.0 4040200.0 
16 1986 20 070.0 8 16 1986 20 13 3.0 565570.0 4039500.0 
16 1986 20 13 0.0 8 16 1986 20 20 3.0 566020.0 4039000.0 
16 1986 20 200.0 8 16 1986 20 25 3.0 566420.0 4038500.0 
8 16 1986 20 300.0 8 16 1986 20 37 3.0 565985.0 4038400.0 


ο co co co Οὐ со OO 
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